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ABSTRACT 


This  paper  documents  the  development  o-f  a  non— para- 
metric statistics  package  -for  the  TRS— 80  microcomputer. 
The  package  is  comprised  of  ten  programs  v*aith  the  major 
emphasis  on  non— parametr i c  hypothesis  testing. 

The  programming  language  is  TRS— SO  Level  II  Disk 
Basic.   The  package  is  compatible  with  any  disk  based 
TRS— SO  Model  I/I  1 1  compatible  microcomputer-   With 
modification  of  the  screen  display  and  the  disk  access 
commands,  the  package  is  transportable  to  microcomputers 
produced  by  other  manufacturers. 

The  statistical  analysis  capability  implemented  on  a 
relatively  inexpensive  system  provides  a  useful  tool  to  the 
student  or  the  trained  analvst- 
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INTRODUCTION 


Statistical  Analysis  o-f  data  is  utilized  as  a  decision 
aid  in  most  every  -field.   Personnel  in  technical  -fields 
have  been  exposed  to  and  require  the  use  o-f  many 
statistical  analysis  procedures  during  the  course  o-f  their 
career.   In  a  career  in  which  statistical  analysis  o-f  data 
is  required  ,  but  not  on  a  daily  basis,  the  knowledge  o-f 
the  proper  statistical  procedure  is  o-ften  retained  where 
the  ioechanics  o-f     implementing  the  procedure  requires 
extensive  review.   A  so-ftware  package  which  eliminates  this 
review  and  handles  the  mechanics  o-f  implementing  the 
procedure  is  most  valuable. 

The  perscnal  computer  has  evolved  to  the  degree  that 
it  is  now  well  suited  -for  the  implementation  o-f  statistical 
analysis  procedures  -for  small  to  moderately  sized  data 
bases.  The  statistical  analysis  so-ftware  that  is 
commercial!'/  available  is  generally  unsophisticated  and 
unsuited  for  the  serious  user. 

An  excellent  package  which  models  parametric 
distributions  and  implements  parametric  hypothesis  testing 
was  developed  by  Mr.  R.  P.  Isbell  and  was  the  subject  o-f 
his  thesis  work  at  the  Naval  Postoraduate  School.   The 


purpose  Q-f  this  thesis  is  to  provide  a  so-ftware  package 
which  io^plements  the  more    commonly  used  non— parametric 
procedures,   accurately  approximates  their  distributions, 
and  provides  a  convenient  data  base  manager  use-ful  -for  both 
this  package  and  that  developed  by  Mr.  Isbell. 


S 


II.   OVERVIEW 

A.   PACKAGE  CONTENTS 

The  package  may  be  logically  divided  into  three 
partitions;  the  data  base  manager,  ncn— parametr i c 
distributions,  and  ncn— parametr i c  hypothesis  tests.   The 
data  base  management  programs  permit  data  entry  and  storage 
to  be  in  a  columnar  -format  or  in  a  tabular  -format.   The 
programs  provide  a  convenient  -file  editing  capability  as 
well  as  a  -flexible  -file  printing  capability. 

The  more  commonly  used  non— parametric  distributions 
and  their  inverse  distributions  ars    approximated.   The 
distributions  3.rez 

1.  Wilcoxson  signed  rank  distribution 

2.  Mann— Whitney  distribution 

3.  Smirnov  distribution 

4.  Kolmogorov  distribution 

5.  Lilliefors  modification  to  the  K— S  test  -for 
Normal  Distribution 

6.  Lillie-fors  modification  to  the  K— S  test  -for 
Exponential  Distribution 

7.  Chi— Square  distribution 

The  package  contains  procedures  -for  implementing  the 
-following  hypothesis  tests: 

1-  Wilcoxson  signed  rank  test 

2.  Mann— Whitney  test 

3.  Smirnov  test 

4.  Test  -for  Normal  Distribution 

a.  Mean  and  variance  unknown 

b.  Mean  known,  variance  unknown 

c.  Mean  and  vairi  3.nc:e    known 

5.  Test  for  Exponential  Distribution 


a.  Parameter  unknown 

b.  Parameter  known 

6.    Contingency  Table  Chi -Square  test 


B.  MICROCOMPUTER  SYSTEM 

The  package  was  developed  on  an  LNW  Model  I 
microcomputer  system  with  48K  bytes  o-f  memory  and  one  5—1/4 
inch  -floppy  disk  drive.   The  package  is  totally  compatible 
with  any  TRS-80  Model  I  or  Model  III  disk  based 
microcomputer  system  and  is  not  dependent  on  any  unique 
-features  o-f  the  advanced  TRS— 80  Model  I/I  1 1  operating 
systems. 

The  programming  language  is  TRS— 80  Disk  Basic 
implemented  by  Microso-ft.   Compatibility  with  all  TRS— SO 
Model  I/I  1 1  systems  is  maintained  by  using  only  TRSDOS 
compatible  commands.   The  package  is  not  directly 
transportable  to  ncn— TRS— SO  compatible  computers;  however, 
the  Microso-ft  interpreter  is  commonly  used  by  other 
manufacturers.   The  display  commands  would  have  to  be 
altered  to  be  compatible  with  microcomputers  produced  by 
other  manu-f acturers. 

C.  LIMITATIONS 

With  any  program  development  on  a  microcomputer,  the 
programmer  must  be  cognizant  of  its  speed,  memory,  and 
precision  limitations.   To  achieve  maximum  utilization  of 
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the  computer  memcry,  the  -floppy  disk  is  used  as  a  virtual 
memory  device  pulling  in  routines  on  an  as  needed  basis. 
This  method  permits  a  minimum  of  3000  data  observations  to 
be  resident  in  the  computer  memory. 

The  speed  of  execution  is  a  major  concern  when  using 
an  interpreted  language.   To  achieve  an  acceptable 
execution  speed,  asymptotic  distributions  were  used  in  lieu 
of  exact  distributions  and  faster  but  more  memory  intensive 
sort  routines  were  routines  were  used. 

Precision  is  certainly  a  concern  in  this  package; 
however,  single  precision  variables  proved  to  be  adequate 
in  most  cases.   Double  precision  variables  were  used  only 
in  the  approximations  to  the  Wilcox son  and  the  Mann— Whitney 
di  str i  but i  ons. 
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III.   DISCUSSION  OF  PROGRAMS  AND  ALGORITHMS 

A.  STRUCTURE 

The  package  consists  o-f  ten  programs,  i.e.  a  control 
module,  two  data  I/O  programs,  a  program  calculating 
distributions  and  six  programs  performing  hypothesis  tests. 
The  choice  o-F  the  program  to  be  executed  is  made  through 
the  menu  driven  control  module.   The  user  is  -free  to 
execute  any  program  in  the  package  through  the  control 
module  with  all  variables  being  initialized  as  a  new 
program  is  executed. 

All  data  I/O  is  conducted  through  the  two  data  I/O 
programs.   The  data  is  entered  to  the  hypothesis  test 
programs  through  a  menu  o-f  established  data  -files  which 
were  created  by  the  data  I/O  programs. 

B.  PROGRAM  1,  MENU/BAS 

MENU/BAS  is  the  control  module  -from  which  all  programs 
Sirs    executed.   The  program  consists  o-f  the  master  menu  and 
two  sub— menus.   The  functional  flow  of  the  program  is  shown 
in  figure  1. 
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BLOCK  DIAGRAM  OF  MENU /B AS 


pATAJORHATJENy 

1.  Coluaar  Data 

2.  Tabular  Data 


Execute 
DBHGR/BAS 


MASTERJENU 

1.  Data  Base  Manager 

2.  Probability  Distributions 

3.  Hypothesis  Tests 


Execute 
CTABLE/BAS 
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Execute 
aiLCOX/BAS 


Execute 
I1ANN/3AS 


Execute 
NP3TAT/BAS 


Execute 
:«IRNQV/9AS 


Execute 
NLILL/3AS 


3 

HYPOTHESIS  TEST  HEHU 

I.  yilcoxscn  Sianed  Ranks  Test 

2.  Hann-yhitnev  Test 

3.  Sflirnov  Test 

4.  Norsal  Test  Distribution 

5.  Exponential  Test  Distribution 

6.  Contincency  Table  Test 

Execute 
ELILL/BAS 


Execute 
CHISSU/BAS 


FIGURE  1. 
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C.   PROGRAM  2,   DBMGR/BAS 

DBMGR/BAS  provides  the  data  I/O  for  the  hypothesis 
test  programs.   The  program  is  menu  driven  and  o-f-fers  the 
-following  selections: 

1.  Create  Mew  Data  Files 

2.  Edit  Data  Files 

3.  Purge  Data  Files 

4.  Print  Data  Files 

5.  Return  to  Master  Menu 

A  special  purpose  data  -file  named  EDITDATA/DIR  is  used 
by  this  program  to  contain  the  -filenames  o-f  the  data  -files 
created.   It  is  used  as  a  data  -file  directory  and  is  called 
by  the  hypothesis  test  programs  to  display  the  filenames  of 
all  current  data  files  and  to  permit  a  convenient  menu 
selection  of  the  data  files. 

C.l.   NEW  FILE  CREATION  ROUTINE 

The  new  file  creation  routine  prompts  the  user  for  the 
number  of  files  to  be  created.   The  user  is  then  prompted 
for  the  disk  drive  number,  i.e  O  —  3,  to  contain  the  files. 
The  data  files  on  the  selected  drive  a.rs    displayed  to 
prevent  any  undesired  duplicate  file  naming.   The  user  is 
then  prompted  for  the  new  filenames.   If  a  previously  used 
filename  is  chosen,  the  user  is  alerted  that  a  duplicate 
•filename  has  been  used  and  the  user  is  given  the  option  of 
renaming  the  file  or  overv-^ri  t  i  ng  the  existing  file. 
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A-fter  the  -files  3ir&    named,  the  input  screen  is 
displayed.   The  screen  will  display  -four  files  (columns) 
and  -five  observations  (rows).   The  rows  and  columns  may  be 
scrolled  to  display  other  observations  as  needed. 

The  input  area  is  located  in  the  top— le-ft  section  o-f 
the  screen.  Depressing  ENTER  transfers  the  input  a^rea.  to 
the  bracketed  destination  area. 


There  are    several  special  -function  keys  which  must  be 
de-f  i  ned: 


ENTER 


Arrows 

Shi -ft  arrow 

D 
I 

CLEAR 
3 


trans-fers  contents  o-f  the  input  arsa. 
to  the  destination  area.   If  the  input 
area  is  blank,  the  destination  Sirsa.    is 
moved  to  the  right. 

moves  the  destination  area 

moves  the  destination  area  to  its 
most  extreme  position 

deletes  the  current  row 

inserts  a  blank  row  immediately  before 
the  current  row 

Returns  to  menu 

writes  the  files  to  disk  and  returns  to 
the  menu 


The  data  is  written  to  disk  in  a  sequential  file  with 
the  first  field  containing  the  number  of  observations  and 
each  field  thereafter  containing  the  observations.   The 
routine  was  designed  for  maximum  user  convenience;  hcv^ever. 


these  convenience  features  will  permit  undesired  "zero" 
valued  observations  i -f  the  user  is  not  aware  o-f  the  design 
c-f  the  program.   Some  precautions  and  guidelines  to  -follow 
are: 


1.   All  displayed  -files  will  be  assigned  a  common 

value  -for  the  number  o-f  observations.   Therefore,  only 
files  which  will  have  the  same  number  of  observations 
should  be  entered  si mutaneousl y . 


2.   It  is  sometimes  tempting  to  "play"  with  the  arrow  keys 
to  move  the  destination  area  around  the  screen.   The 
number  of  observations  is  determined  by  the  maximum  row 
reached  by  the  destination  Sirea.^     therefore,  "playing" 
with  the  arrow  keys  may  result  in  an  unwanted  value  for 
the  number  of  observati ons- 


C.2.   EDIT  EXISTING  FILES  ROUTINE 

The  file  editing  routine  operates  similiarly  to  the 
nevw  file  creation  routine  with  cnly  minor  variation.   The 
user  is  prompted  for  the  disk  drive  containing  the  files 
and  then  prompted  for  the  selection  of  the  files  to  be 
edited.   The  files  a.re    displayed  on  the  screen  and  editing 
is  performed  exactly  as  described  in  the  new  file  creation 
routine.   Prior  to  writing  the  files  to  disk,  the  user  is 
given  the  opportunity  to  rename  the  file. 

C.3.   PURGE  FILES 

The  data  files  ars    displayed  one  by  one  and  the  user 
is  asked  if  it  is  desired  to  erase  the  file  from  the  disk. 


16 


All  erasing  of  data  -files  should  be  performed  from  the 
purge  routine  and  not  from  the  operating  system  routines  ai 
the  purge  routine  also  removes  the  filename  from  the  data 
file  directory,  EDITDATA/DIR. 

C.4.   PRINT  DATA  FILES 

This  routine  permits  formatting  the  data  files  for 
printing.   The  user  is  prompted  for  the  files  to  be 
printed,  the  printer  line  width,  column  headers,  and  given 
the  option  of  numbering  the  observations.   The  maximum 
number  of  spaces  between  columns  is  calculated  and 
displayed.   The  user  is  prompted  for  the  desired  number  of 
spaces  between  columns.   Finally  the  user  is  asked  for  the 
reoort  title  and  the  files  ar^    orinted. 


D.   PROGRAM  3,  CTABLE/BAS 

CTABLE/3AS  is  used  to  create  tables  used  in  the 
contingency  table  test  program.   The  program  operates 
similiarly  to  the  DBMGR/BAS  program  with  the  difference 
listed  below: 

1.  The  data  file  names  are    stored  in  CTABLE/DIR  versus 
EDITDATA/DIR. 

2.  In  the  new  file  creation  routine,  the  user  is 
prompted  for  the  number  of  rows  and  columns  in  the 
table. 

3.  The  "D"  and  "I"  kevs  have  no  function. 
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E.   PROGRAM  4,  WILCGX/BAS 

WILCOX/BAS  performs  the  WilcoKBon  Signed  Rank  test 
using  data  stored  in  data  -files  created  by  DBMGR/BAS.   The 
user  is  prompted  -for  the  two  data  files  for  the  test.   The 
test  statistic  is  computed  and  the  user  is  prompted  for  the 
hypothesis  testing  problem,  i.e.  one  sided  or  two  sided- 
The  test  statistic,  sample  size  and  p— value  a.re    displayed 
and  the  user  is  prompted  for  the  significance  level  of  the 
test.   The  significance  level  is  compared  with  the  p— value 
and  acceptance  cr  rejection  of  the  null  hypothesis  is 
determined. 

E.l.   I'JILCOXSQN  SIGNED  RANK  TEST 

The  data  consists  of  two  paired  samples  of  n 
cbservat  i  ons,  (  x  ^  ,  y-?)  ,  (xo,  y->)  ,  .  .  ,  ^'^n  '  ''^n '  * 
The  difference  d^=  >'i~yi  is  computed  for  each  of 
the  n  pairs.   Pairs  with  a  difference  of  zero  3irs    omitted. 
The  absolute  differences,  Idj,  !,  of  the  remaining  pairs 
are  ranked  with  rank  1  assigned  to  the  pair  ()t-,y^) 
with  the  smallest  absolute  difference-   If  one  or  more 
pairs  have  the  same  absolute  difference,  the  average  rank 
is  assigned  to  the  pairs.   The  following  test  statistics 
are    computed: 

T   =  Sum  of  ranks  with  positive  d^  ""  s 

T   =  Sum  of  ranks  with  negative  d^ ' s 


IS 


T   =  MIN(T'*".T  ) 


The  assumptions  -for  the  test  are: 

1.   The  distribution  o-f  each  d^  is  symmetric 
2-   The  pairs  (>:^,y^)  constitute  a  bivariate  random 
samp 1 e . 

The  hypothesis  may  be  stated  in  one  o-f  three  ways: 


A. 

^o 

: 

E(X) 

= 

E(Y) 

•^1 

m 

E'X) 

7^ 

E(Y) 

B. 

Ho 

m 
m 

E(X) 

< 

E(Y) 

Hi 

E(X  ) 

E(Y) 

C. 

'  'O 

m 

E'X) 

•-. 

E<Y) 

^O 

■ 

E  ( X ) 

< 

E  ( Y ) 

The  decision  rule  for  the  three  hypotheses  is  to  accept 

H^  i-F  : 

A.   P''Test  Statistic  <  T)       Significance  Level /2 
S.   PCTest  Statistic  ■•!  T)   !'■■   Sinni-ficance  L.evei 
w«   P'CTest  Statistic  'C    T  )       Si'^ni-ficance  Level 

E.2.   DISTRIBUTION  THEORY  AND  ALGORITHM 

The  distribution  o-f  the  Wilcoxscn  statistic.  T.  m3.y    bt 
■found  by  randomizing  the  signed  ranks.   Since  there  are  twt 
possible  signs,  +  and  — ,  and  n  ranks,  the  number  o-f 
combinations  in  the  randomization  is  2*^.   The  exact 
cumulative  distribution  may  be  computed  by  examining  the 
2   combinations  and  counting  those  with  the  sum  c-f  ranks 
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less  than  or  equal  to  the  test  statistic  and  dividing  by 
2*^.   Such  a  method  is  cbvioLisly  slow  for  all  but  very 
small  sample  sizes. 

A  number  o-f  appro>;i  mati  ens  to  the  exact  distribution 

have  been  proposed.   The  most  accurate  o-F  the  proposed 

approximations  is  the  Edgewcrth  approximation  which  is 

documented  in  Re-ferences  1.  2,  and  3.   This  approximation 

is  based  upon  the  fact  that  the  distribution  of  T  is 

asymptotically  normal  with  mean  and  variance: 

u-f-  =   n  (n-^1 )  /4 

0:^2=   n  (n+1)  <2n  +  l) /24 

The    normalized    I'iilcoxscn    Statistic,     z,     is    defined: 

z    =     (T    —    Uj)  /  <r^ 

The  Edgewcrth  expansion  to  terms  of  size  l/n'-  is 

P(T  <  t)  =  F(z>  +    L^f  *^'*  (z) /4!  + 

I  ,  i  ^  5  ^  ,'  -  'k  /  A  !  +^,-.  t'  L  ,  )  -^-c  '"  ^  ^  /  '^  ' 

where  F(z)  is  the  cumulative  standard  normal  distribution 

and  : 

f  *     is  the  i''th  deri^/ative  of  the  normal  densitv 

f unct i  on 

L4  =  -12  (3n'^+3n-l ) /5n  (n  +  l)  ■:2n-i-i  ) 

L.  =  576(3n^+6n^-3n+l) /7<n (n+1)  <2n+l) )^ 
o 
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The  algorithm  impiementsd  in  this  package  uses  a  speed 
optimized  e^act  distribution  when  the  execution  time  is 
small  and  the  Edgeworth  approximation  otherwise.   The 
decision  rule  is  : 

Approximation       i -f     T  >  27  or  (N  >  9  and  T  >  N) 

Exact  i -f     otherv-iise 

The  algorithm  must  use  the  normal  cumulative 
distribution-   The  algorithm  used  in  this  package  -for  the 
normal  cumulative  distribution  was  taken  -from  the  thesis 
work  a-f    Mr.  R.P.  Isbell  and  is  documented  in  Re-ference  4. 


F.   PROGRAM  5,  MANN/BAS 

MANN/EAS  per-forms  the  Mann-Whitney  test  using  data 
stored  in  data  -files  created  by  DBMGR/BAS.   The  user  is 
prompted  for  the  two  data  -files  for  the  test.   The  test 
statistic  is  computed  and  the  user  is  prompted  for  the 
hypothesis  testing  problem,  i.e.  one-sided  ar    tv^jo— si  ded. 
The  test  statistic,  sample  sizes  and  p— value  are  displayed 
and  the  user  is  oromoted  for  the  sionif icance  level  of  the 
test.   The  significance  level  is  compared  with  the  p— value 
and  acceotance  or  rejection  of  the  null  hvDothesis  is 
deter mi  ned . 


F.i.     mann-iajhitney  test 

The  data  consists  o-f  two  randcm  samples, 

X  j^,;:2?  •  •  ■  ?"n  ^'^^  V'l  ?  ¥2?  •  •  *  ?  y^-   ^^^  ^'^° 

samples  3.rs    combined  and  ranked  -from  1  to  n+m  with  the 

smallest  observation  being  assigned  the  rank  o-f  1.   The 

test  statistics  are  : 

T,.  =  nm  +  n(n+l)/2  -  (Sum  o-f  ranks  -for  saiiiple  a) 

T   =  nm  -»-  m<m+l)/2  —  (Sum  o-f  ranks  for  sample  y) 

T   =  MIN(Tj,  ,T,^) 

The  assumptions  for  the  test  a.rs    : 

i.   Both  samples  Ars    random  samples- 

2.   The  two  samples  3.r&    independent. 


T-C  there  is  a  dif-ference  bet  ween  the  two 
population  distributions,  that  difference  is 

V 


in  location  only,  i.e.  F,.  (u)  =  F,  m-i-c! 


where  c  is  a  constant. 
The  hypothesis  may  be  stated  in  one  of  three  ways 

A .  H^  :  E ( X )  =  E ( Y ) 
H 1  :  E  (  X  )  ?^  E  ( Y ) 

B.  Hq  :  E(X)  i;  E(Y) 
H  ^  :  E  (  X  )  >  E  ( Y ) 

C.  hi^     :  E(X)  ^  E'Y) 
H.  :  E(X)  <  E(Y) 


The  decision  rule  for  the  three  hypotheses  is  to  accep' 


A.  P(Test  Statistic  <  T)   >    Significance  Level/ 

B.  P(Test  Statistic  <  T„ )       Significance  Level 


C.   P(Test  Statistic  <  T_)       Siqnificance  Lev 


y 


■::}• 


.1 


F.2.   DISTRIBUTION  THEDRY  AhiD  ALGORITHM 

The  distribution  o-f  the  Mann— Whitney  statistic  may  be 
■found  by  randomizing  the  ranks  of  the  combined  data.   There 
3.rs     (n+m)  !  /  (n  !  )  (m!)  "ways  in  which  the  ranks  may  be  grouped. 
The  e;<act  distribution  may  be  computed  by  examining  all 
possible  combinations  o-f  ranks  and  counting  those  with  the 
rank  sum  less  than  or  equal  to  nm+n(n+l)/2  minus  the  test 
statistic  and  dividing  by  (n+m)  I /  (n i  )  (m!).   Such  a  method 
is  slow  -for  all  but  very  small  sample  sizes. 

A  num.ber  o-f  approximations  to  the  exact  distribution 
have  been  proposed.   The  most  accurate  o-f  those  proposed  is 
the  Edgeworth  approximation  documented  in  Re-ference  5.   The 
approximation  is  based  on  the  fact  that  the  test  statistic 
is  asvmtct i cal 1 V  normal  with  mean  and  variance  : 


Ua^  =  nm/2 

OX  =  nm  (n-t-m-^1 )  /  12 


The  normalized  Mann— Whitney  Statistic,  z,  is  defined  : 


■^t '  '  ^ 
The  Edgeworth  expansion  is  : 

P<t<t)  =  FCz)  -i-L^f  '^Mz) /4'  +  L^f  ''"^  <z) /6!  -<- 

35(L4)-f ,7j (z)/S! 

where  F(z)  is  the  standard  normal  cumulative  distribution 


■f'''^   is  the  i'th  derivative  c-f  the  standard  normal 

density 
L^  =  -  (jr,'^+n'^+n+!Tj++nm)  /  (20nm  (n+m+1 )  ) 
L^  =  NUM  /  DENOM 
where  : 
NUM  =  2  (n'*+m'*)H-4nm  (n^+fn^) +6n'^fTi'^ 

-f-4  (n~'+fT«"^)  -i-7nm  <n+m)  -^n'^+fn'^H-2nm— n— m 
DENOM  =  210n^fn-(n+m-f-l)^ 

The  algorithm  that  is  implemented  in  this  package  use? 
a  speed  optimized  exact  distribution  -for  small  sample  size-; 
and  the  Edgewcrth  approK i mat i en  for  larger  sample  sizes. 
The  decision  rule  is  : 

appro;:  i  mat  i  on      if       m  >  9  or  T  >  27  —  n(n+l)/2 
exact  i-f       otherwise 

where  n  and  m  ars    the  smaller  and 

M    tZl  1   '~^  v_  1      ^3  w^k  1  a  t  Laf  ^  i^    —e  ^.    ^—   i.^  —j  i   >— .  —^  *—f  \—,  ^-^    ^  i>   «  «^  ^~      w 

G.   PRCGRAM  6,  3MIRN0V/EAS 

3MIRNGV/BA3  per-forms  the  Smirnov  test  on  data  stored 
in  data  "files  bv  IDSMGFc/'SAS.   The  user  is  or  emoted  -for  the 
two  data  -files  -for  the  t'^s*^.   The  test  statistic  is 
computed  and  the  user  is  prompted  for  the  hypothesis 
testing  problem,  i.e.  one— sided  or  two-sided.   The  test 
statistic,  sample  sizes,  and  p— value  3.r^    displayed.   The 
user  is  prompted  for  the  significance  level  of  the  test. 
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The  significance  level  o-f  the  test  is  compared  to  the 

p— value  and  acceptance  or  rejection  o-f  the  null  hypothesi 

is  determined. 

G. 1.   SMIRNOV  TEST 

The  data  consists  o-f  two  independent  random  samples 
with  unknown  cumulative  distribution  -functions  F(x)  a.nd, 
G(::)-   Sj(x)  and  S--v(x)  are  their  empirical  cumulative 
distribution  -functions.   The  test  statistics  a.rs    : 


„  + 


=  MAX (Si (x)  -  So^k) 


D   =  MAX (S2(-)  -  3^  (x)  ) 


D   =  MAX(D"^,D  ) 


The  assumptions  -for  the  test  3.rs    : 

1.  The  samples  airs    random  samples 

2.  The  tv^o  samples  are  independent 

3.  The  random  samcies  ars    continuous 


The  hypotheses  may  be  stated 


^. 


B. 


H^  :  F(x)  =  G(x 
H^  :  F(x)  ^  G'x 


hT  :  F(x)   >  G(x 


^■7 


F<x)  1    G(x 
F  ( X  )  <  G  C  X 


The  decision  rul 
H„  if  : 

A.  PCTest  Statistic  <   D) 

B.  PCTest  Statistic  <  D^) 


for  the  three  nvoothese  is  to  acceot 


Sianificance  Level /2 


Significance  Level 


C.   PCTest  Statistic  <  D  )   >   Significance  Level 

G.2.   DISTRIBUTION  THEORY 

The  exact  distribution  of  the  test  statistic  is 
di-f-ficult  to  analyze  -for  the  general  case.   E.F.  Dricn,  in 
Reference  6,  derived  the  distribution  for  the  special  case 
of  equal  sample  sizes.   It  is  known  that  the  test  statistic 
has  the  asymptotic  distribution  : 

G<S)  =  1  -  ^   {-l)^~'^sy.pK-2i-s'^) 
where  G(S)  is  known  as  the  Smirnov  distribution 

S  =  D   SQR ( (n+m) /nm)   where  n  and  m  are  sample  sizes 

It  is  generally  agreed  that  sample  sizes  of  100  or 
more  3.rs    required  for  the  overestimating  effects  of  the 
asymptotic  distribution  to  become  negligible.   The  problem 
is  then  to  improve  upon  the  Smirnov  distribution  fcr  sampls 
sizes  of  less  than  100.   P.J.  Kim,  in  Reference  7, 
presented  an  improvement  which  reduced  the  overestimating 
effects  for  the  smaller  sample  sizes.   The  improvements 
presented  by  P.J.  Kim  and  implemented  in  this  package  ^rs     : 
S  =  D  ♦  (SDRvnm/ (n-*-m)  )  ) -s-2/ (3*  (SDR(m)  )  )  for  m=kn,  k  =  l,2...... 

S  =  <D-l/2m) • (SQR(n) )  for  n/m  <  0.1 

S  =  D -(SQR  (nm/ (n-^m)  )  )  -*-  2/ (5  •  (SQR  (n )  )  )  otherwise 
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H.   PROGRAM  7,  NLILL/BAS 

NLILL/EAS  performs  a  test  -for  normal  dietributicn  on  a 
data  set  stored  in  a  data  -file  created  by  DBMGR/BAS.   The 
test  may  be  per-fcrmed  for  the  three  cases  : 

1.  mean  and  variance  unknown 

2.  mean  known  and  variance  unknown 
3-   mean  and  variance  known 

The  user  is  prompted  for  the  data  file  to  be  tested 
and  prompted  for  the  specific  case  to  be  tested.   In  the 
cases  where  a  parameter  is  known,  the  user  is  prompted  for 
the  known  parameter.   The  test  statistic,  sample  size  and 
p— value  3.rs    displayed.   The  user  is  prompted  -far    the 
significance  level  a-f    the  test.   The  significance  level  is 
compared  with  the  p— value  and  acceptance  or  rejection  of 
the  null  hypothesis  is  determined. 

H.l.   TEST  FOR  NORMAL  DISTRIBUTION 

The  data  consists  of  a  random  sample  >:■;,....  x^ 
from  an  unknown  cumulative  distribution  F(;;).   G<>;>  is  the 
normal  cumulative  distribution  with  parameters  either 
specified  or  estimated  from  the  samel e.  deoending  uoon  the 
specific  case.   The  test  statistics  s.rs    : 

D"^  =  MAX  (i/n  -  Giy.^  )  ) 

D~  =  MAX<G(x.)  -  <i-l)/n)) 

D   =  MAX  (D'*".D~) 


Z/ 


The  assumption  -for  the  test  is  that  the  sample  is  a 

random  sample-   The  hypotheses  sirs    : 

Case  1.    Hj_,  :  F(m)  is  normal  with  unspecified 

mean  and  variance 
Hf     :  F(x)  is  not  normal 

Case  2.    H   :  F-'x)  is  normal  with  specified  mean 
H<  :  F<k)  is  not  normal  with  specified 
mean 

Case  3.    H   :  F<>:)  is  normal  with  specified  mean 

and  variance 
H-t  :  F''>:)  is  not  normal  with  specified 
mean  a^nd    ^y 3.r i  ainas 

The  decision  rule  is  accept  H„  if  : 

P(Test  Statistic  <  D)  >  Significance  level 

H.2.   DISTRIBUTION  THEORY  AND  ALGORITHM 

The  distribution  of  the  test  statistic  differs  for 
each  of  the  three  cases.   For  Case  1,  where  the  mean  and 
variance  3.rs    unknown,  the  distribution  was  approximated 
using  Monte  Carlo  techniques  by  H.W.  Lilliefors  and 
published  in  Reference  8.   The  distribution  is  presentd  in 
tabular  form  and  is  not  described  in  an  analytical  sense. 
In  Reference  S,  Lilliefors  recognized  that,  for  a 
particular  significance  level,  the  distribution  of  the  tesi 
statistic  differed  by  a  multiplicative  constant  from  that 
of  the  Kolmogorov  distribution.   This  package  capitalizes 
on  this  concept;  however,  the  multiplicative  value  is  a 
function  of  the  value  of  the  cumulative  distribution.  The 
solution  must  therefore  be  solved  iterativelv.   The 
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algorithm  -follows: 


T  < Test  Statistic 

Denominator  < 1 

Do  Until  ABS<Tdencm  -  Denominator)  <  .0001 

D  < T/Denominator 

S  < D:*(-.01  +  SDR(N)  -i-  O.S5/SQR(N)) 

Tdenom  < Denominator 

CDF  < 1  -  (. 135)2 

I-F  CDF  <  .85 

Then  Denominator  < 1  —  .24>i CDF 

If  .35  <  CDF  <  .90 

Then  Denominator  < 1.337  —  .c6xCDF 

I-F  .90  <  CDF  <  .  95 

Then  Denominator  < 1.339  —  .64;;CDF 

I-f  .95  <  CDF  <  .975 

Then    Denominator    < 1.775    —    1.12i<CDF 

I-f    CDF    >    .975 

Then  Denominator  < 2.133  —  1.467xCDF 

End  Do 


For  Case  2,  where  the  mean  is  known  and  the  variance 

is  unknown,  the  distribution  was  approximated  by  M.A. 

Stephens  and  published  in  Re-ference  9-   A  si  mi  liar  solution 

was  -found  to  exist  -for  this  case  as  existed  -for  Case  1, 

i.e.  the  cumulative  distribution  must  be  calculated 

iteratively.   The  algorithm  -for  the  iterative  solution  is  : 

T< Test  Stastistic 

Denominator  < 1 

Do  Until  AES (Tdenom  —  Denominator)  <  .0001 

D  < T/Denominator 

S  < D:t(.35  +  SQR(N)  +  .2S3/SQR(N)) 

Tdenom  < Denominator 

CDF  < 1  -  (. 135)" 

I-f  CDF  <  .9 

Then  Denominator  < 1.33  —  .3xCDF 

I-f  .9  <  CDF  <  .95 

Then  Denominator  < 1.47  —  .4xCDF 

I-F  .95  <  CDF  <  .975 

Then  Denominator  < 1 . 35  —  .3xCDF 

I-f  CDF  >  .975 

Then  Denominator  < 1.07 

End  Do 
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For  Case  3,  where  the  mean  and  variance  3.rs    completely 
specified,  the  distribution  is  the  Kclmogorov  distribution. 
The  Kolmogorov  distribution  is  a  special  case  o-f  the 
Smirnov  distribution  with  one  sample  size  approaching 
in-finity  and  has  the  asymptotic  distribution  : 


(S)  =  1  -  2  ^-l)^ 


li^v 


EXP<-2i^S^) 
with  S  =  D-3QR(N) 


M.A.  Stephens  proposed  improvements  to  the  asymptotic 
distribution  in  Re-ference  9  which  improve  the  accuracy  of 
the  asymptotic  distribution  for  small  sample  sizes.   This 
package  implements  the  improvements  proposed  by  Stephens  in 
Reference  9.  i.e.  : 

S  =  D  •  (  .  12  ^    3QR  (n  )  h-  .11  /SGR  (n  )  ) 


1.   PROGRAM  S,   ELILL/BAS 


ELILL/BAS  performs  a  test  for  exoonential  distribution 
on  data  stored  in  data  files  by  DBMGR/BAS.   The  test  may  be 
perf or-Tied  for  the  cases  where  the  oarameter  is  known  and 
where  the  parameter  is  unknown. 

The  user  is  prompted  for  the  data  file  to  be  tested 
and  for  the  specific  case,  i.e.  parameter  known  cr    unknown. 
In  the  case  where  the  parameter  is  known,  the  user  is 
prompted  for  the  parameter  value.   The  test  statistic. 


sample  size  and  p— value  &rs    displayed-   The  user  is 
prompted  -For  the  si  gni -f  i  cance  level  o-f  the  test.   The 
signi-f icance  level  is  compared  with  the  p— value  and 
acceptance  or  rejection  o-f  the  null  hypothesis  is 
determi  ned. 
I.l.   TEST  FOR  EXPONENTIAL  DISTRIBUTION 

The  data  consists  c-f  a  random  sample  >;•,,....  x^, 
-from  some  unknown  cumulative  distribution  F(>i).   G(;c)  is 
the  exponential  cumulative  distribution  -function  with  the 
parameter  either  estimated  +rom  the  data  or  specified, 
depending  upon  the  speci-fiic  case.   The  test  statistics 


3.rsz 


D"^  =  MAX  (i/n  -  G(>;)  ) 

D~  =  MAX  ( G  ■:  >: )  -  (  i  - 1 )  /n  ) 

D   =  MAX  ■D'^,D~) 

The  assumption  for  the  test  is  that  the  sample  is  a 

random  sample.   The  hypotheses  are: 

Case  1.  Hj_j  :  FCx)  is  exponential  —  parameter  specified 
H-,     :  F(x)  is  not  exponential 

Case  2.  H^    :  F<x)  is  exponential  —  parameter  unspecified 
H^  :  F(x)  is  not  exponential 

The  decision  rule  is  to  accept  H   if  : 

P(Test  Statistic  <  D)  >  Significance  Level 

1.2-   DISTRIBUTION  THEORY  AND  ALGORITHM 

The  distribution  of  the  test  statistic  differs  for  thi 
two  cases.   For  the  case  where  the  parameter  is  unknown. 


the  distribation  was  appro;: ijua ted  using  Monte  Carlo 
techniques  by  H.W.  Liliiefors  and  published  in  Re-ference 
iO.   M-A.  Stephens,  in  Reference  9,  expanded  the  work  c-f 
Lif-fie-fors  by  providing  improvements  for  the  small  sample 
cases.   The  algorithm  is  an  iterative  method  which  assumes 
a  solution  and  computes  the  test  statistic  vjhich  would 
yield  the  assumed  value  for  the  cumulative  distribution 
function  (CDF).   The  assumed  CDF  is  updated  until  the 
difference  between  the  calculated  test  statistic  and  the 
actual  test  statistic  is  sufficiently  small.   The  algorithm 


T  < TEST  STATISTIC 

DENOMINATDR  < .26  -^    SQR(n)  +  .5/SQR(N) 

F  i  < o 

CTT*  .•■■ 1 

DO  UNTIL  ERROR  <  .0001 

CDF  < (F2  -  Fl)/2 

IF  CDF  >  .95 

THEN  TS  < 5.  125;<CDF  -  3.7S0S 

IF  .90  <  CDF  <  .95 

THEN  TS  < 2.  02:; CDF  -  .8367 

IF  .30  <  CDF  <  . 90 

THEN  TS  < 1.124xCDF  -  .0249 

IF  .70  <  CDF  <  .30 

THEN  TS  < 0.743:<CDF  +  .2799 

IF  .50  <  CDF  <  .  70 

THEN  TS  < 0.546;<CDF  +  .4137 

IF  .30  <  CDF  <  .50 

THEN  TS  < 0.464xCDF  +  .4538 

IF  .20  <  CDF  <  .30 

THEN    TS    < 0.435>;CDF    -f-    .4525 

IF    .10    <     CDF    <     .20 

THEN  TS  < 0.537): CDF  ^     .4321 

IF  CDF  <  . 10 

THEN  TS  < 0.346;: CDF  +  .4062 

TS  < TS/DENOMINATDR  +  . 2/N 

ERROR  < T-TS 

IF  ERROR  >  O  THEN  Fl  < CDF  ELSE  F2  < CDF 

END  DO 


For  Case  2,  with  the  parameter  known,  the  distribution 
is  the  Kolmogorov  distribution.   The  distribution  and  the 
algorithm  used  in  this  package  is  described  in  the  -fully 
specified  case  of  the  test  for  normal  distribution. 

J.   PROGRAM  9,   CHIBQU/BAS 

CHISDU/BAS  performs  the  r  ;<  c  contingency  table  test 
using  data  stared  in  data  files  created  by  CTABLE/SAS.   The 
user  is  prompted  for  the  data  file  to  be  tested.   The  test 
statistic,  degrees  of  freedom,  fraction  of  cells  with 
expected  frequency  less  than  5,  number  of  cells  with 
expected  frequency  less  than  1,  and  the  p— value  are 
printed.   The  user  is  prompted  for  the  significance  level 
of  the  test.   The  significance  level  is  compared  with  the 
p— value  and  acceptance  \Dr    rejection  of  the  null  hypothesis 
is  determined. 
J.l.   CONTINGENCY  TABLE  TEST 

The  contingency  table  test  may  be  used  to  present  a 
tabulation  of  data  contained  in  several  samples  to  test  the 
hypothesis  that  the  probabilities  do  not  differ  from  sample 
to  sample.   Another  use  of  the  r  ;<  c  contingency  table  is 
with  a  single  sample  where  each  element  in  the  sample  may 
be  classified  into  one  of  r  different  categories  according 
to  one  category  and  into  one  of  c  different  categories 
according  to  a  second  category.   The  two  applications  Eiirs 


treated  the  same  in  the  statistical  analysis,  but 
basic  di -f -f  erences  between  the  two  applications  justify 
separate  discussions- 
J.1.1   TEST  FOR  DIFFERENCES  IN  CELL  PROBABILITIES 

The  data  consists  o-f  random  samples  drawn  -from  r 
populations.   Each  observation  in  each  sample  may  be 
classified  into  one  of  c  different  categories.   Let  Oj  ^ 
be  the  number  of  observations  from  the  i ' th  sample  that 
fall  into  category  j.   Let  n^  be  the  number  of 
observations  from  the  i ' th  population,  c^  be  the  number 
of  observations  in  the  j ' th  category,  and  N  be  the  total 
number  of  observations.   The  test  statistic  is  : 
T    =      E    E(0,  ^  -  E-  ,)2/E.  , 
where  E^  ^  =  n^c^/N 

The  assumptions  ars    that  each  sample  is  a  random 

sample,  the  samples  3.r&    mutually  independent,  and  each 

observation  may  be  categorized  into  exactly  one  of  c 

categories.   The  hypotheses  may  be  stated  : 

H   :  All  of  the  probabilities  in  the  same  column 
a.rB    equal 

H,  :  At  least  two  of  the  probabilities  in  the  same 
column  a.r&    not  equal 

The  decision  rule  is  to  accept  the  null  hypothesis  if  : 

P''Test  Statistic  <  T)   >  significance  level 


J.  1-2   TEST  FOR  INDEPENDENCE 

The  data  consists  o-f  one  random  sample  o-f  size  N-   The 
data  may  be  classi-fied  by  one  o-f  two  criteria.   Using  the 
-first  criteria,  each  observation  is  associated  with  one  c-f 
r    rows  and  using  the  second  criteria,  each  cbservatin  is 
associatd  with  one  of  the  c  columns.   Let  O-  ^  be  the 
number  o-f  observations  in  row  i  and  column  j.  n-  be  the 
number  o-f  observations  in  the  i  '  th  row.  and  c^  be  the 
number  c-f  observations  in  the  j '"  th  column.   The  test 
statistic  is  as  de-fined  in  the  test  -for  di -f -f  erences  in  cell 
probabi 1 i ties. 

The  hypothesis  may  be  stated  : 

H   :  P(row  i,  column  j)  =  P(row  i)  x  P (column  j) 
H-t  :  P(row  i-  column  j)  =  P(row  i)  k  P  (column  j) 

The  decision  rule  is  to  accept  the  null  hypothesis  i -f  : 
PvTest  Statistic  <  T)   >   Significance  Level 

J. 2   DISTRIBUTION  THEORY  AND  ALGORITHM 

The  exact  distribution  o-f  T  is  di-f-ficult  to  calculate 
-for  all  but  the  2  x  2  case.   It  is  well  known  that  T  has 
asymptotically  a  Chi— Square  distribution.   The  Chi— Square 
approximation  is  valid  -for  large  exected  cell 
probabilities,  i.e.  E^^  ^ .   The  approximation  ^s 


considered  to  be  satisfactory  i -f  no  Ej^  ^  is  less  than  1 
and  not  more  than  2.07.    of  the  Ej^  ^'s  3.rs    less  than  5. 


The  algorithm  used  in  this  package  for  the  Chi— Square 
distribution  was  extracted  from  the  thesis  of  Mr.  R.P. 
Isbeli.   The  details  of  his  appro;:  i  mat  ion  to  the 
distribution  ars    detailed  in  Reference  1.   For  the  2  x  2 
case,  a  correction  factor  proposed  by  F.  Yates  in  Reference 
11  is  used  to  improve  the  Chi— Square  appro;-:  i  mati  on  to  the 
e;:act  distribution. 

K.   PROGRAM  10,  NPSTAT/BAS 


NPSTAT/BAS  contains  seven  non— parametr i c  probability 
distributions,  the  normal  distribution,  and  their  inverse 
probability  distributions.   The  distributions  appro;;  i  mated 

Normal  Distribution 

Wilco;:son  Signed  Rank  Distribution 

Mann— Whitney  Distribution 

Smirnov  Distribution 
5.   Kolmogorov  Distribution 
■^■-        Li  1  lief  or  s  Test  for  Normal  Distribution 

Lilliefcrs  Test  for  E;;oonential  Distribution 


The  user  is  orovided  a  menu  to  select  the  distributior 
desired.   The  user  is  prompted  for  the  required  inputs  and 
the  cumulative  distribution  is  displayed. 

The  approximations  to  the  distributions  3ir&    described 
in  detail  in  the  hypothesis  test  prcgram.s.   The  3.cc\j.r 3.cy 


c=o 


ccmpariscns  o-f  this  packages  approK  i  mati  ons  to  the  exact 
distributions  aire    provided  in  Tables  1  through  6. 
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TABLE  1 
WILCOXSON  DISTRIBUTION  ACCURACY  COMPARISON 


Sample 
Size 

Pcsi0.05 

Pc^rO.  10 

P^0.20 

Exact 

Approx 

Exact 

Approx 

Exact 

Approx 

5 

.0313 

.0313 

.  0938 

.  0938 

.  1 563 

.  1563 

10 

.  0420 

.  0420 

.  0962 

.0970 

.  1875 

.  1875 

15 

.  0473 

.0473 

.0938 

.0938 

.  1947 

.  1947 

20 

.0487 

.  0487 

.0947 

.0947 

.  1942 

.  1942 

25 

.0479 

.  0479 

.  0954 

.0954 

.  1942 

.  1942 

30 

.0481 

.0481 

.  0990 

.  0990 

.  1965 

.  1965 

35 

.0484 

.  0484 

.0977 

.  0977 

.  1972 

.  1970 

40 

.0486 

.  0486 

.0984 

.  0984 

.  1987 

.  1987 

45 

.0493 

.  0498 

-0981 

.0981 

.  1973 

.  1973 

50 

.0495 

.0495 

.  0988 

.  0988 

.  1994 

.  1994 

38 


TABLE  2 
MANN-WHITNEY  DISTRIBUTION  ACCURACY  COMPARISON 


Sample 

Sampl e 

Test 

Size 

Size 

Statistic 

E;<ajit 

ApproK 

6 

6 

-7 

.  0465 

.  0463 

6 

6 

10 

.  1 20 1 

.  1200 

6 

6 

12 

-  1970 

.  1965 

6 

8 

11 

.  0539 

.  0542 

6 

3 

14 

.  1142 

.1143 

6 

3 

16 

.  1725 

.  1724 

e 

Q 

16 

.  0524 

.  0524 

S 

3 

20 

.1172 

.  1173 

8 

3 

-~>'T* 

.  1641 

.  1641 

S 

1 0 

21 

.  0506 

.  0506 

3 

10 

25 

.  1015 

.  1016 

S 

10 

23 

.  1577 

.  1 577 

10 

10 

23 

.  0526 

.  0526 

10 

10 

•Jy-j> 

.  loss 

-  1033 

10 

10 

36 

.1575 

.  1575 

TABLE  3 
SMIRNOV  DISTRIBUTION  ACCURACY  COMPARISON 


Probabl ity 

=  0.95 

n/m 

m  = 

25 

m  = 

50 

m  = 

100 

Exact 

Approx 

Exact 

Approx 

Exact 

Approx 

.  1 

.5209 

.5185 

.3797 

.3737 

_  2 

.5321 

.  5339 

.  3903 

.3910 

.  2830 

.  4£l3oO> 

.  o 

.  3393 

.3434 

.2443 

.2463 

.4 

.4188 

.4277 

.  3044 

3086 

.2192 

.  2213 

.5 

.2769 

.2765 

.  2003 

.  2003 

.6 

.  3704 

.  3733 

.2680 

.2694 

.  1925 

.  1932 

.  7 

.2579 

.2571 

.  1848 

.  1844 

.8 

.3423 

.  3429 

.2474 

.2475 

.  1774 

.  1775 

.9 

.  2409 

.  2397 

.  1726 

,  1719 

1.0 

.  3087 

.3082 

.2256 

.2258 

.  1 634 

.  1635 

Probabl ity 

=  0 . 975 

n/m 

m  = 

25 

fn  = 

50 

1 

fn  =  10< 

J 

Exact 

.Approx 

Exact 

Approx 

Exact    Af 

3prox 

.  1 

.5775 

.5743 

.4224 

.4142 

'n 

.5910 

.5996 

.4348 

.  4375 

.3152 

.3161 

■  -Jj 

.3777 

.  3829 

.2719 

.2742 

.4 

.4764 

.4779 

.  3390 

.3441 

.  2440 

.2464 

.  5 

.  3096 

.  3094 

.2236 

.  6 

.4126 

.4172 

.2985 

.  3004 

.2143 

.2151 

.7 

.2871 

.2867 

.  2056 

.  2053 

.8 

.3818 

.3832 

.  2754 

.2759 

19744 

.  1976 

.  9 

.2682 

.2673 

.  1920 

.  1914 

1 

,0 

.3464 

.3462 

.2524 

.2526 

.  1822 

.  1325 

40 


TABLE  4 
KOLMOGOROV  DISTRIBUTION  ACCURACY  COMPARISON 


Sampl e 
Size 

P  = 

0.90 

P  = 

0.95 

P  = 

0 .  975 

Exact 

ApproM 

E>:act 

Appro:: 

Exact 

Approx 

5 

.447 

.446 

.509 

.509 

.  563 

.564 

10 

.323 

•  Oj£.C> 

.369 

.369 

.409 

.409 

15 

.266 

.267 

-304 

.  304 

.333 

•  ooS 

20 

.232 

.265 

.265 

.294 

.294 

.208 

.209 

.238 

.233 

.264 

.264 

30 

.  194 

.  191 

.218 

.213 

.  242 

.  242 

35 

.  177 

-  177 

.  202 

.202 

.  "^24 

.  224 

40 

.  165 

.  166 

.  189 

.  189 

.210 

.210 

41 


TABLE  5 

LILLIEFORS  NORMAL I TY  TEST  DISTRIBUTION 
ACCURACY  COMPARISON 


Sampl e 

P  =  0.90 

P  =  C 

.95 

P  =  ■: 

-     m        r       f 

Size 

E;;act 

Appro; 

Exact 

Appro;; 

Exact 

Appro; 

5 

.315 

.314 

.337 

.  343 

-  405 

.396 

i  0 

.239 

.239 

.258 

.261 

.  294 

.  302 

15 

.201 

.  200 

.220 

.219 

.257 

.253 

20 

.  174 

.  176 

.  1 90 

.  192 

.231 

TT?'? 

25 

.  15S 

.  159 

.  173 

.  173 

.  200 

.  200 

30 

.  144 

.  146 

.  161 

.  159 

.  137 

.  134 

42 


TABLE  6 

LILLIEFORS  EXPONENT I AL I TY  TEST  DISTRIBUTION 
ACCURACY  COMPARISaN 


P  =  0 

.80 

P  =  0 

.90 

P  =  0 

.95 

Sample 

Sire 

E>;act 

Approx 

Exact 

Apprcx 

Exact 

Approx 

5 

.  3603 

.  3615 

.  4045 

.  4028 

.  4420 

.  4400 

10 

.  2626 

.2642 

.2955 

.2956 

.3244 

.3239 

15 

.2174 

.2135 

.2448 

.2448 

.  2690 

.  2636 

20 

.  1S93 

.  1  905 

.2132 

.2137 

.2345 

.  2346 

25 

.  1703 

.  1711 

.  1918 

1921 

.2110 

.211 0 

30 

.  1559 

.  1567 

.  1756 

.  1760 

.  1932 

.  1933 

35 

.  1447 

.  1454 

.  1  630 

.  1633 

.  1:^93 

1795 

40 

.  1356 

.  1362 

.  1528 

.1531 

.  1681 

.  1683 

45 

.  12S1 

.  1286 

.  1443 

.  1445 

.  1583 

.  1589 

50 

.1217 

.  1221 

.1371 

.  i  o  /  o 

.  1509 

.  1510 

60 

.1113 

.1117 

.  1255 

.  1256 

.1381 

.1381 

70 

.  1033 

.  1035 

.1164 

.1164 

1281 

.  1281 

SO 

.  096S 

.  0969 

.  1090 

.  1  09  1 

.  1 200 

.  1200 

90 

.0914 

.0914 

.  1029 

-  1029 

.1132 

.1132 

1 00 

.  0868 

.  0868 

.  0977 

.  0'=?77 

.  1075 

.  1075 

10  CLEAR  350 

20  DEFINT  I-K  :  DEFS7R  Z 

30  LB$=CHR*<123)  :  RB$=CHR* ( 125)  :  MB$=CHR* ( 176)  : 

BL$=CHR*(252) 

40  CLS: PRINT  TAB (10)  "NON-PARAMETRIC  STATISTICAL  ANALYSIS 

SOFTWARE" 

50  PRINT  TAB (31)  "by" 

60  PRINT  TAB (21)  "Robert  Lee  Zangmeister" 

70  PRINT  :  PRINT  TAB(13)  "Submitted  in  partial  -ful-fillment 

D-f  the"  :  PRINT  TAB  (17)  "requirements  -for  the  degree  o-f " 

80  PRINT  :  PRINT  TAB (19)  "MASTER  OF  APPLIED  SCIENCE" 

90  PRINT  :  PRINT  TAB  (28)  "-from  the" 

100  PRINT  :  PRINT  TAB (19)  "NAVAL  POSTGRADUATE  SCHOOL" 

110  PRINT  TAB (25)  "December  1982" 

120  PRINT: PRINT  TAB (17)  "Press  SPACE  BAR  to  continue" 

130  ZI=INKEY$  :  IF  ZI=""  THEN  130 

140  CLS  :  PRINT  TAB (20) ; "PROGRAM  SELECTION  MENU" 

150  PRINT  5)128,  "(1)    Data  Base  Manager" 

160  PRINT  "(2)    Probability  Distributions" 

170  PRINT  "(3)    Hypothesis  Tests" 

180  KP=896  :  KS=6  :  GOSUB  500 

1 90  ON  I O  GOTO  200 , 250 , 1 000 

200  CLS: PRINT  "DATA  BASE  OPTION" : PRINT 

210  PRINT  "(1)    Column  Format  Data  Entry" 

220  PRINT  "(2)    Table  Format  Data  Entry" 

230  KP=896:KS=2: GOSUB  500 

240  IF  10  =  1  THEN  RUN  "DBMGR/BAS"  ELSE  RUN  "CTABLE/3AS" 

250  RUN  "NPSTAT/BAS" 

500  REM   *  *  *  KEYBOARD  ENTRY  ROUTINE  *  * 

510  PRINT  3KP+5, "OPTION  DESIRED  " ; LB* ; MB$; RB*: 

520  ZI=INKEY*  :  IF  ZI=""  THEN  520ELSE  IO=VAL(ZI) 

530  PRINT  SKP+21,  ZI;:  FOR  1=1  TO  100  :  NEXT 

540  IF  I0<1  OR  IO>KS  THEN  PRINT  S)KP,  BL$;:  PRINT  3KP,"t* 

NOT  A  VALID  OPTION    ttt" ; z     FOR  1=1  TO  1000  :  NEXT  :  PRINT 

S)KP,  BL$;  :  GOTO  510 

550  RETURN 

1000  CLS: PRINT  TAB (20) ; "HYPOTHESIS  TEST  MENU": PRINT 

1010  PRINT  3128, "(1)   Wilcox  son  Signed  Ranks  Test" 

1020  PRINT  "(2)   Mann-Whitney  Test" 

1030  PRINT  "(3)   Chi-Square  Contingency  Table  Test" 

1040  PRINT  "(4)   Smirnov  Test" 

1050  PRINT  "(5)   Test  -For  Normality" 

1060  PRINT  "(6)   Test  -for  E;:ponenti  al  i  ty " 

1070  KP=896:KS=6: GOSUB  500 

1080  ON  10  GOTO  1100,1110,1120,1130,1140,1150 

1100  RUN  "WILCOX/BAS" 

1110  RUN  "MANN/BAS" 

1120  RUN  "CHISQU/3AS" 

1130  RUN  "SMIRNGV/BAS" 

1140  RUN  "NLILL/BAS" 

1150  RUN  "ELILL/BAS" 
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i^-FtOCBFt*^^!^   S     —     lI>BMC3F^/:B<=i^ 


10  CLEAR  1000 
20  ON  ERROR  GOTO  19000 
30  DEFSTR  Z 

40  DIM  F$(95) ,F1*(11) ,M*(11) 

50  LB*=CHR«  < 1 23 ) : RB*=CHR* ( 1 25 ) : MB$=CHR* ( 1 76 ) : BL*=CHR$ ( 252 ) 
60  CLS 

70  PRINT  TAB ( 20 >; "DATA  FILE  MANAGEMENT  PROGRAM" 
80  PRINT:PRINT:PRINT  "MENU:" 
90  PRINT  "(1)   Create  New  Data  Files" 
100  PRINT  "(2)   Edit  Existing  Files" 
110  PRINT  "(3)   Purge  Existing  Files" 
120  PRINT  "(4)   Print  Data  Files" 
130  PRINT  "(5)   Return  to  Master  Menu" 
140  KP=S96:KS=5:G0SUB  ISO 
150    ON  10  GOSUB  260,310,330,340,170 
160     GOTO  60 
170  RUN  "MENU/BAS" 

ISO  REM  t    t       KEYBOARD  INPUT  SUBROUTINE  * 
190  PRINT  5)KP+5,  "OPTION  DESIRED  "  ;  LB$;  MB*;  RB$; 
200  ZI=INKEY«:IF  ZI=""  THEN  200  ELSE  IO=VAL(ZI) 
210  PRINT  3KP+21,ZI; :FOR  1=1  TO  100: NEXT 

220  IF  I0<1  OR  IO>KS  THEN  PRINT  3KP, BL$; : PRINT  3KP,"«* 
NOT  A  VALID  OPTION   **«";:  FOR  1=1  TO  1000:  NEXT:  PRINT 
3KP,  BL$;:GOTO  190 
230  RETURN 

240  DN*=INKEY$: IF  DN*=""  THEN  240  ELSE  I=VAL(DN5) 
250  IF  KO  OR  I>3  THEN  240  ELSE  RETURN 

260  ^SUB  TO   PERFORM  NEW  DATA  FILE  CREATION  ROUTINES 
270  CLS  :  INPUT  "How  Many  Variables  Sire    to  be  Used 
(1-11) ";NF:CT=NF:RT=1 
280  IF  NF<1  OR  NF>1 1  THEN  270 

290  GOSUB  14000  :  GOSUB  12000  :  N=l  :  CR=1  :  GOSUB  6020 
300  CLS: GOSUB  15050: RUN 
310  GOSUB  350  :  CR=1  :  GOSUB  6020 
320  GOSUB  15000  :  RUN 
330  GOSUB  5000  :  RUN 
340  GOSUB  350  :  GOSUB  4000  :  RUN 
350  REM  t    t    t       READ  MULTIPLE  ARRAYS  t    t 
360  K=0  :  DL$="/" 

370  CLS: PRINT  "Disk  Drive  Containing  Files  (0-3)?" 
380  GOSUB  240 

390  SL=17 15: GOSUB  1200: IF  SL=0  THEN  GOTO  370 
400  GOSUB  1400  :  PRINT  "" 

410  PRINT: IF  10  =  2  PRINT  "Enter  the  Files  To  Be  Edited 
(e.g.  1/3/4)";  ELSE  PRINT  "Enter  the  Files  To  Be  Printed"; 
420  INPUT  B$  :P1=0 
430  P2=INSTR(1,BS,DL$) 
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440  IF  P2=0  THEN  NV=VAL(B*):  Pl=l  :  GOTO  460 

450     NV=VAL(LEFT$(B$,P2-1) )  :  B*=RIGHT$ (B$, LEN (B$) -P2) 

460  IF  NV>FO  THEN  PRINT: PRINT"  *  «  «  INVALID  FILE  #  « 

«":GOTO  410 

470     K=K+1  :  IF  K=12  THEN  NF=11  :  GOTO  530 

480     F1*(K)=F*(NV) 

490     IF  P1=0  THEN  430 

500  B*="": PRINT: PRINT  "Do  You  Wish  To  Read  Files  From 

Another  Di5k(Y/N)?  "; 

510  ZI=INKEY$:IF  11=""    GOTO  510  ELSE  IF  2I="Y"  PRINT 

ZI:GaTO  370 

520  NF=K:CT=K: PRINT  "N" 

530  MM=-1 

540  FOR  1=1   TO  NF 

550  OPEN  "I",1,F1$(I)  :  INPUT  #1,M  :  CLOSE  1 

560  IF  M  >  MM  THEN  MM=M 

570  NEXT  I 

580  RT=MM 

590  GOSUB  12000  :  PRINT 

600  FOR  J=l  TO  NF 

610      PRINT  "READING  IN  FILE   — >  ";F1$(J) 

620      OPEN  "I",1,F1$(J)  :  INPUT  #1,N 

630      FOR  1=1  TO  N  :  INPUT  #1,X(I,J)  :  NEXT  I 

640      CLOSE  1 

650  NEXT  J 

660  N=MM  :  RETURN 

1000  REM   DIRECTORY  CHECK   *  * 

1010  IF  F0=0  THEN  RETURN 

1020  P2=0 

1030  FOR  Pl=l  TO  FO 

1040      OPEN  "I", 1,F$(P1-P2) 

1050      CLOSE  1 

1060  NEXT  PI 

1070  IF  P2=0  THEN  RETURN  ELSE  GOSUB  3400  :  RETURN 

1100  REM  *  «  *   CORRECT  DIRECTORY   *  t 

1110  IF  P1-P2=F0  THEN  1130 

1120     FOR  I=P1-P2  TO  FO-1  :  F$ (I ) =F$ (I+l )  :  NEXT  I 

1130  F0=F0-1  :  P2=P2+1  :  RETURN 

1200  REM   *  *  «   READ  DIRECTORY   *  t 

1210  CLS 

1220  F$="EDITDATA/DIR: "+DN$ 

1230  OPEN  "I", 1,F* 

1240  INPUT  #1,F0 

1250  IF  F0=0  THEN  1290 

1260  FOR  1=1  TO  FO 

1270       INPUT  #1,F$(I) 

1280  NEXT  I 

1290  CLOSE  1  :  GOSUB  1000 

1300  CLOSE  1 

1310  RETURN 

1400  REM  t    t    t       DISPLAY  DIRECTORY  t    t 
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1410  CLS 

1420  IF  F0=0  THEN  1510 

1430  FOR  1=1  TO  FO  STEP  4 

1440      FOR  J=l  TO  4 

1450  IF  I+J-1>F0  THEN  1510 

1460  IF  J=l  AND  LEN(F*<I+J-1) )<12  AND  I>1  THEN 

PRINT  "" 

1470         PRINT  TAB(16*(J-1) ) ; " (";STR$(I+J-1) ; ") 

";F*(I+J-1) ; 

1480      NEXT  J 

1490  NEXT  I 

1500  PRINT 

1510  RETURN 

2300  REM   *  *  *   CREATE  DIRECTORY   *  « 

2305  CLOSE  1 

2320  F$="EDITDATA/DIR: "+DN* 

2340  OPEN  "O", 1,F* 

2345  PRINT  #1,0 

2360  CLOSE  1 

2370  RETURN 

3400  REM  t    t    t       ENTER  NEW  FILES   *  * 

3410  F*="EDITDATA/DIR: "+DN* 

3420  OPEN  "O", 1,F$ 

3430  PRINT  #1,F0 

3440  FOR  1=1  TO  FO 

3450      PRINT  #1,F*(I) ; ", "; 

3460  NEXT  I 

3470  CLOSE  1 

3480  RETURN 

3500  REM  *    t    t       FILE  NAME  CHECK   *  * 

3510  EX=0 

3520  FOR  1=1  TO  FO 

3530       IF  F$=F$(I)  THEN  EX=1  :  RETURN 

3540  NEXT  I 

3550  RETURN 

3600  REM  t    t    t       DIRECTORY  DISPLAY   *  « 

3610  CLS 

3620  INPUT  "DISK  DRIVE  NUMBER  (0,1,2,3)  TO  FIND  DIRECTORY 

OF";DN$ 

3630  IF  VAL(DN$)<0  OR  VAL(DN$)>3  THEN  3610 

3640  GOSUB  1200  :  SOSUB  1400  :  PRINT  ""  :  INPUT  "PRESS 

ENTER  TO  CONTINUE" ;J1 

3650  RETURN 

4000  REM  *    t     t       PRINT  DATA  FILES   *  * 

4005  CLS  :  INPUT  "Input  Printer  Width";PW 

4010  PF=0  :  ML=-1  :  bL$=" ! "  :  TL=-1 

4020  FOR  1=1  TO  N 

4030      FOR  J=l  TO  NF 

4040  IF  LEN(STR$(X (I, J) ) )  >  ML  THEN 

ML=LEN ( STR$ ( X  < I , J ) ) ) 

4050      NEXT  J 
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4060  NEXT  I 

4061  ML=ML+1 

4070  CLS: PRINT: PRINT  "Enter  Column  Headers.   Separating 

lines  with" 

4080  PRINT"the  delimiter  ' ! '  to  produce  a  two  column 

header" 

4100  FOR  1=1  TO  NF 

4110  PRINT:PRINT  "Enter  the  Column  Header  -for  Column";  I 

4120  INPUT  M$(I) 

4130  NEXT  I 

4140  MS=-1 

4150  FOR  1=1  TO  NF 

4160      P1=0 

4170      P2=INSTR(P1+1,M*(I) ,DL$) 

4130      L=P2-P1-1 

4190       IF  P2=0  THEN  L=LEN (M$ ( I ) ) -PI 

4200       IF  L>MS  THEN  MS=L 

4210       IF  P2=U  THEN  4240 

4220  P1=P2 

4230  GOTO  4170 

4240  NEXT  I 

4250  IF  ML>MS  THEN  MX=ML  ELSE  MX=MS 

4260  FOR  1=1  TO  NF 

4270      P1=0  :  NL=0  :  T$="" 

4290      P2=INSTR(P1+1,M*(I) ,DL$) 

4300      L=P2-P1-1 

4310       IF  P2=0  THEN  L=LEN (M$ ( I ) ) -PI 

4320      SF=INT( (MX-L) /2)  :  SB=MX-L-SF 

4330      T$=T*+STRING$(SF, " 

")-HP1ID*(M*(I)  ,P1  +  1,L)+STRING$(SB,  "  ") 

4335      NL=NL+1  :  IF  NL>TL  THEN  TL=NL 

4340      IF  P2=0  THEN  M$ ( I ) =T*  :  GOTO  4370 

4350         P1=P2 

4360  GOTO  4290 

4370  NEXT  I 

4380  PRINT: PRINT  "Would  you  like  to  number  the  observations 

(Y/N)?  "; 

4385  A$=INKEY$:IF  A$=""  THEN  4385  ELSE  PRINT  A* 

4390  IF  A*="Y"  THEN  PI =5  ELSE  PI=0 

4400  MS=INT( (PW-MX«NF-PI) /NF) 

4410  PRINT:  PRINT"The  maximum  #  o-f  spaces  that  can  be  used 

between  columns  is";MS;"." 

4420  PRINT:  INPUT"Enter  the  number  o-f  spaces  between 

columns" ; SP 

4445  IF  SP<0  THEN  INPUT  "TO  RENAME  COLUMNS  TYPE  Y  ELSE 

N";T$  :  IF  T$="Y"  THEN  4010  ELSE  RETURN 

4450  IF  SP<0  OR  SP>MS  THEN  RETURN 

4455  PRINT: INPUT  "Enter  Report  Title" ;T$ 

4472  P1=(NF-1)  *  (MX+SP)-t-PI+MX-l 

4474  Pl=INT(Pl/2)  :  P2=PW-P1  :  IF  LEN(T$)>PW  THEN 

T*=LEFT$(T*,PW) 
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4476  IF  LEN(T*)>P2  THEN  P1=INT ( (PW-LEN (T*) ) /2) 

4480  LPRINT  TAB(P1);T* 

4490  LPRINT  TAB(Pl);  :  FOR  1=1  TO  LEN(T$)  :  LPRINT 

"-";:NEXT  I  :  LPRINT 

4500  LPRINT  ""  :  LPRINT  "" 

4505  FOR  J=l  TO  TL 

4510      FOR  1=1  TO  NF 

4520  LPRINT 

TAB( (I-1)*(MX+SP)+PI) ;MID*(M$(I) , ( J-1 ) X+1 , MX) ; 

4530      NEXT  I 

4540      LPRINT  "" 

4550  NEXT  J 

4560  A$=""  :  FOR  1=1  TO  MX  :  A$=A$+"-"  :  NEXT  I 

4570  FOR  1=1  TO  NF 

4580      LPRINT  TAB ( (I-l ) « (MX+SP) +PI ) ; A*; 

4590  NEXT  I 

4600  LPRINT  ""  :  K=PI+MX-1-INT ( (MX-ML) /2) 

4610  FOR  1=1  TO  N 

4620       IF  PI=5  THEN  LPRINT  STR* ( I ) ; " : " ; 

4630      FOR  J=l  TO  NF 

4640  LPRINT 

TAB( (J-1) «(MX+SP)-LEN(STR$(X (I, J) ) )+K)  ; X  (I, J) ; 

4650      NEXT  J 

4660      LPRINT  "" 

4670  NEXT  I 

4680  RETURN 

5000  REM  t    t    t       PURGE  DATA  FILES   *  * 

5010  CLS 

5020  PRINT  "Disk  Drive  Containing  Files  to  Be  Purged 

(0-3)?" 

5030  GOSUB  240 

5040  SL=5010  :  GOSUB  1200  :  IF  SL=0  THEN  GOTO  5010 

5045  K=0 

5046  PRINT  TAB(20);"  tttttt    PURGE  DATA  FILES  **:**" 
5050  FOR  1=1  TO  FO 

5060  A*="": PRINT  F*(I-K);"      KILL  IT?  (Y/N/Q)  "; 

5061  ZI=INKEY$:IF  ZI=""  THEN  5061 

5062  IF  ZI="Q"  THEN  PRINT  Z I : GOTO  5095 

5064  IF  ZI="Y"  THEN  PRINT  Z I : GOSUB  5100  ELSE  PRINT  "N" 
5090  NEXT  I 
5095  GOSUB  3400 

5099  RETURN 

5100  KILL  F$(I-K) 

5110  IF  I-K=F0  THEN  5130 

5120     FOR  J=I-K  TO  FO-1  :  F$ ( J ) =F$ ( J+1 )  :  NEXT  J: 

5130  F0=F0-1  :  K=K+1 

5140  RETURN 

6000  REM  ttt    EXAMINE  DATA  tt 

6005  INPUT  CR 

6010  IF  CR  <  1  THEN  CR=1 

6015  IF  CR  >  N  THEN  CR=N 


49 


6020  E*=CHR*(31)  :  N$=""  :P1=0:  CLS 

6025  C1=1:R1=1:B1=1 

6030  CR=R1:CC=C1 

6035  GOSUB  9500 

6040  GOSUB  6300  :  GOSUB  6350 

6045  P1=0 

6050  A*=INKEY$  :  IF  A$=""  THEN  6050 

6055  IF  A*=CHR*<91)  AND  P1=0  THEN  N*="":  GOSUB  6500  :  GOTO 

6050 

6060  IF  A*=CHR*(27)  AND  P1=0  THEN  N$=" ": CR=1 : GOSUB 

6520: GOTO  6050 

6065  IF  A$=CHR$(10)  AND  P1=0  THEN  N*="":  GOSUB  7000  :  GOTO 

6050 

6070  IF  A$=CHR*(26)  AND  P1=0  THEN  N$=" ": CR=RT: GOSUB 

7020:G0T0  6050 

6075  IF  A$=CHR$(3)  AND  P1=0   THEN  N$="":  GOSUB  7500  :  GOTO 

6050 

6080  IF  A*=CHR*(9)  AND  P1=0  THEN  N$="":  GOSUB  SOOO  :  GOTO 

6050 

6085  IF  A*=CHR*(24)  AND  P1=0  THEN  N$=" ": CC=1 : GOSUB 

7530: GOTO  6050 

6090  IF  A$=CHR$(25)  AND  P1=0  THEN  N$=" ": CC=CT: GOSUB 

8030:G0T0  6050 

6095  IF  (A$>"/"  AND  A$< " : " )   OR  A*="."  OR  A*="-"  OR  A*="+" 

OR   A$="E"  THEN  Pl=l  :  GOSUB  9000  :  GOTO  6050 

6100  IF  A$=CHR$':i3)  THEN  P1=0:  GOSUB  S500  :  GOTO  6050 

6105  IF  A*=CHR$(a)  AND  Pl=l  THEN  GOSUB  9020  :  GOTO  6050 

6110  IF  A$="D"  AND  P1=0  THEN  N$="": GOSUB  6320: GOSUB 

10000: GOTO  6050 

6115  IF  A$="I"  AND  P1=0  THEN  N$=" ": GOSUB  6320: GOSUB 

10500: GOTO  6050 

6120  IF  A$=CHR«(31)  AND  P1=0  THEN  RUN 

6125  IF  A$="3"  AND  P1=0  THEN  N=RT  :  RETURN 

6130  GOTO  6050 

6300  REM  t    t     t       CLEAR  ENTRY  AREA  t    t 

6310  CC=1 

6320  PRINT  ^2, " 

6330  PRINT  5)  2-CHR$(95> 

6340  RETURN 

6350  REM   *  *  *   BRACKET  ENTRY  *  :* 

6360  K=324+12S«(CR-Rl)-*-15:<(CC-Cl) 

6370  TY=15360+K 

6380  POKE  TY, 183  :  POKE  TY+14, 107 

6390  RETURN 

6400  POKE  TY, 123  :  POKE  TY+14, 128 

6410  RETURN 

6500  REM  t    t    t       PROCESS  UP  ARROW   tc  * 

6510  CR=CR-1:IF  CR< 1  THEN  CR=1:G0T0  6540 

6520  IF  CR<R1  THEN  R1=CR: RL=R 1+4: GOSUB  9500 

6530  GOSUB  6400: GOSUB  6350 

6540  RETURN 
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7000  REM   *  «  *   PROCESS  DOWN  ARROW   «  * 

7010  CR=CR+1 

7020  IF  CR>RL  THEN  RL=CR: Rl=RL-4: GOSUB  6400:G0SUB  9500 

7030  IF  CR  >  RT  THEN  RT=RT+1 

7040  GOSUB  6400: GOSUB  6350 

7050  RETURN 

7500  REM   *  *  *   PROCESS  LEFT  ARROW  «  « 

7510  CC=CC-1 

7520  IF  CC<1  THEN  CC=1:G0T0  7550 

7530  IF  CC<C1  THEN  C1=CC: CL=C 1+3: GOSUB  9500 

7540  GOSUB  6400: GOSUB  6350 

7550  RETURN 

8000  REM  *    t    *       PROCESS  RIGHT  ARROW   *  * 

8010  CC=CC+1 

8020  IF  CC  >  CT  THEN  GOSUB  9200 

8030  IF  COCL  THEN  CL=CC:  CI  =CL-3:  GOSUB  9500 

8040  GOSUB  6400: GOSUB  6350 

S050  RETURN 

8500  REM  t    t    t       CHANGE  VALUE  IN  FIELD   *  * 

8510  GOSUB  6320 

8520  IF  N$='*"  THEN  GOSUB  8000  :  RETURN 

8530  IF  CR>RT  THEN  RT=CR 

8540  X ( CR , CC ) =VAL ( N* ) 

8550  N*="" 

8560  K=325+12S* (CR-Rl ) +15* (CC-Cl ) 

8570  PRINT  5)K+1,X(CR,CC)  ; 

8580  RETURN 

8590  K=K-2:PRINT  SK, CHR* ( 128) ; 

8600  PRINT  USING  ZZ;X(CR,CC); 

8610  RETURN 

9000  REM  t    t     $■  PROCESS  NUMBER  ENTRY   *  :* 

9010  N*=N$+A*  :  PRINT  S)2,  "  "  :  PRINT  3  2,N$;:PRINT 

CHR«(95)  :  RETURN 

9020  IF  LEN(N*)>0  THEN  N*=LEFT$  <N$,  LEN  (N*) -1 )  :  PRINT  :D  2," 

":  PRINT  a  2,N* 
9030  RETURN 

9200  REM  t     t    t       LAST  COLUMN  PROCESSING   *  t 
9210  CC=1  :  CI  =  CC  :  CL=Cl+3 
9220  CR=CR+1 

9230  IF  CR>RL  THEN  RL=CR  :  Rl=RL-4 
9240  GOSUB  6400 

9250  IF  CT  >  4  OR  CR  =  RL  THEN  GOSUB  9500 
9260  IF  CR>RT  THEN  RT=RT+1 
9270  RETURN 

9500  REM  t    *    t       SCREEN  DISPLAY   «  « 

9510  PRINT  5)196,  CHR*(  143)  :  IF  CT<Cl+4  THEN  CL=CT  ELSE 
CL=Cl+3 
9520  RL=Rl+4 
9530  FOR  I=C1-1  TO  CL-1 

9540  PRINT  3197+15* (I+l-Cl) ,STRING$ (6, CHR$( 143) ); :PRINT 
CHR* < 128) ;: PRINT  CHR$ (65+1 );: PRINT  CHR$ < 128) ;: PRINT 
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STRING* (A, CHR*( 143) ) 

9550  NEXT  I 

9560  PRINT  ;D195,CHR$(170)  :PRINT  3)259,  CHR$  (  170) 

9570  FOR  I=R1-1  TO  RL-1 

9580  IF  I  >  3  THEN  K1=0  ELSE  Kl=l 

9590  PRINT  3320+128* ( I+l-Rl ) +K1, STR$ < I+l ); zPRINT  CHR*(170) 

9600  FOR  J=C1-1  TO  CL-1 

9610  PRINT  5)326+15«(J  +  l-Cl)+12S*(I  +  l-Rl)  ,  X(I  +  1,  J  +  1) 

9620  NEXT  J 

9630  PRINT  3387+128* ( I+l-Rl ) ,CHR*( 170) 

9640  NEXT  I 

9650  RETURN 

10000  REM  t    t    t       DELETE  A  LINE   *  * 

10010  FOR  I=CR  TO  RT-1:F0R  J=l  TO  CT: X ( I , J ) =X ( I+l , J )  :  NEXT 

J  :  NEXT  I 

10020  FOR  J=l  TO  CT:  X(RT,J)=0  :  NEXT  J 

10030  RT=RT-1  :  IF  Rl  >  1  AND  RL  >  RT  THEN  R1=R1-1  : 

CR=CR-1 

10040  GOSUB  9500  :  GOSL'B  6350 

10050  RETURN 

10500  REM  t    t     t        INSERT  A  LINE   *  * 

10510  FOR  I=RT  TO  CT  STEP  -1  :  FOR  J=l  TO  CT  : 

X  (I  +  l, J)=X (I, J)  :  NEXT  J  :  NEXT  I 

10520  FOR  J=l  TO  CT  :  X(CR,J)=0  :  NEXT  J 

10530  RT=RT+1  :  GOSUB  9500  :  GOSUB  6350 

10540  RETURN 

11000  REM   *  *  *   VERIFY  FILENAME   *  # 

11010  FW=0 

11020  P1=INSTR(1,S$, "/") 

11030  IF  PI  <>  O  THEN  IF  LEN(S$)-P1>3  THEN  FW=1  :  RETURN 

11040  IF  PI  >8  THEN  FW=1  :  RETURN 

11050  IF  P1=0  THEN  IF  LEN(S*)>8  THEN  FW=1  :  RETURN 

11060  IF  PI  <>  0  THEN  IF  ASC (MID$ (S$, Pl+1 , 1 > )<65  OR 

ASC(MID*(S$,P1+1, 1) ) >90  THEN  FW=1  :  RETURN 

11070  IF  ASC(L£FT$(S$,  1)  )  ::65  OR  ASC  (LEFT$  (S*,  1 )  )  >90  THEN 

FW=1  :  RETURN 

llOSO  IF  P1=0  THEN  P2=LEN(3$)  ELSE  P2=P1-1 

11090  P3=l 

11100  FOR  I=P3  TO  P2 

11110        A$=MID$(S$, I, 1) 

11120        IF  <ASC(A$)>47  AND  ASC(A$)<58)  OR  (ASC(A$)>64 

AND  ASC(A$)<91  )  THEN  11140 

11130  FW=1  :  RETURN 

11140    NEXT  I 

11150  IF  P3=i  AND  PI  <>  O  THEN  P3=P1+1  :  P2=LEN<S$)  :  GOTO 

11100 

11160  RETURN 

12000  REM  t    t    t       DIMENSION  ARRAYS  «  * 

12010  NV=MEM-500  :  NV=INT (NV/ (4* (NF+1 ) ) ) -1 

12020  IF  MM  >  NV  THEN  PRINT  "INSUFFICIENT  SPACE  TO  READ  IN. 

": INPUT  "PRESS  ENTER  TO  CONTINUE"; A* 


12030  DIM  X (NV,NF) 

12040  MN=NV 

12050  RETURN 

13000  REM  t    t    t  RENAME  FILES  ROUTINE   *  * 

13010  CLS 

13020  PRINT  "Disk  Drive  To  Contain  Files  (0-3)?" 

13030  GOSUB  240 

13040  SL=8240: GOSUB  1200  :  GOSUB  1400  :   PRINT  "" 

13050  FOR  J=l  TO  NF 

13060      PRINT  "Variable  #" ; J; "Currentl y  Uses  Filename 

— >";F1$<J) 

13070      INPUT  "New  Name";S$  :  GOSUB  11000 

13080      IF  FW=1  THEN  PRINT  "INVALID  FILENAME"  :  GOTO 

13060 

13090         F$=S$  :  GOSUB  3500 

13100         IF  EX=1  THEN  PRINT: PRINT  "A  File  Currently  Uses 

the  New  Fi 1 ename. ": PRINT  "Do  You  Want  to  Use  it  Anyway 

•:y/n)?"  else  goto  13120 

13110  A$=INKEY*: IF  A$=""  THEN  13110  ELSE  PRINT  A*:  IF  A$<> 

"Y"  THEN  13060 

13120  F1$(J)=F* 

13130  NEXT  J 

13140  RETURN 

14000  REM   *  «  *   NAME  FILES  ROUTINE   *  * 

14010  CLS 

14020  PRINT  "Disk  Drive  To  Contain  Files  (0-3) ?": GOSUB  240 

14030  SL=8330  :  GOSUB  1200  :  GOSUB  1400  :  PRINT  "" 

14040  PRINT 

14050  FOR  J=l  TO  NF 

14060  PRINT  "FILENAME  FOR  VARIABLE  #";J  :  INPUT  S$  :  GOSUB 

11000 

14070  IF  FW=1  THEN  PRINT  "INVALID  FILENAME"  :  GOTO  14060 

14080     F$=S$  :  GOSUB  3500 

14090     IF  EX=1  THEN  PRINT  "FILENAME  IS  ALREADY  USED.  TYPE 

Y  TO  USE  NAME  ANYWAY  ELSE  N";:  INPUT  A*  :  IF  A*  <>  "Y"  THEN 

1 4060 

14100     F1$<J)=F$ 

14110  NEXT  J 

14120  RETURN 

15000  REM   *  *  «   EXITING  ROUTINE  FOR  OPTIONS  1  ?<  2  t    t 

15010  CLS 

15020  PRINT  "Do  You  Wish  To  Rename  Files  (Y/N)?" 

15030  ZI=INKEY$:IF  ZI=""  GOTO  15030 

15040  IF  2I="Y"  THEN  GOSUB  13000 

15050  SL=8550  :  GOSUB   1200  :  GOSUB  1400  :  PRINT  "" 

15060  PR I NT: PR I NT 

15070  FOR   K=l  TO  NF 

15080         PRINT  "WRITING   TO  DISK  — >  ";F1$(K) 

15090       OPEN  "O",  1,F1S(K)+":  "-HDN$ 

15100       PRINT  #1,N  :  FOR  1=1  TO  N  :  PRINT  #1,X<I,K);  : 

NEXT  I 
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15110       CLOSE  1  :  F*=F1*(K)  :  GOSUB  3500  :  IF  EX=1  THEN 

15130 

15120       F0=F0+1  :   F* (FO) =F14 (K)  :  GDSUB  3400 

15130  NEXT  K 

15140  RETURN 

19000  'ERROR  HANDLING  ROUTINE 

19015  IF  ERR=10  AND  ERL=8520  THEN  N*=""  :  RESUME  8530 

19020  IF  ERR  <>  106  THEN  191 OO 

19021  IF  ERL=5100  THEN  RESUME  5110 

19022  IF  ERL=1040  THEN  GOSUB  1100  :  RESUME  1050 
19025  IF  ERL  <>  1230  THEN  19100 

19030  IF  SL=1715  OR  SL=5010  OR  SL=5410   THEN  PRINT  "NO 

FILES  ON  THIS  DRIVE.  PRESS  ENTER  TO  CONTINUE":  INPUT  AS  : 

SL=0  :  RESUME  1300 

19045  IF  SLIsA8240  OR  SL=8330   OR  SL=S550  THEN  GOSUB  2300  : 

RESUME  1210 

19100  PRINT  "ERR=";ERR; "ERL=";ERL  :  ON  ERROR  GOTO  O  : 

RESUME 
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i=-Fi:oc3Ft<=»i^  :3  _   c:"rrf=%:BL_E:/]Bi!=i^ 


10  CLEAR  lOOO 

20  ON  ERROR  GOTO  19000 

30  DEFSTR  Z 

40  DIM  F*(95) ,M*(11) ,N*(11) 

50  LB*=CHR« ( 1 23 ) : RB$=CHR« (125): MB*=CHR$ (176): BL*=CHR* ( 252 ) 

60  Z$="######"  :  ZZ  =  "##.  ##i:  ZLL" 

70  CLS 

80  PRINT  TAB (20) ; "DATA  FILE  MANAGEMENT  PROGRAM" 

90  PR I NT : PR I NT : PR I NT  " MENU : " 

100  PRINT  "(1)   Create  New  Data  Files" 

110  PRINT  "(2)   Edit  Existing  Files" 

120  PRINT  "(3)   Purge  Existing  Files" 

130  PRINT  "(4)   Print  Data  Files" 

140  PRINT  "(5)   Return  to  Master  Menu" 

150  KP=S96:KS=5:G0SUB  190 

160    ON  10  GOSUB  1000,1100,1200,1300,130 

170     GOTO  70 

180  RUN  "MENU/BAS" 

190  REM  *  t       KEYBOARD  INPUT  SUBROUTINE  * 

200  PRINT  3KP+5, "OPTION  DESIRED  " ; LB$; MB$; RB$: 

210  ZI=INKEY*:IF  ZI=""  THEN  210  ELSE  IO=VAL(ZI) 

220  PRINT  3KP+21, ZI; :FOR  1=1  TO  100: NEXT 

230  IF  I0<1  OR  IO>KS  THEN  PRINT  S)KP,  BL$;  :  PRINT  5)KP,"*S 

NOT  A  VALID  OPTION   ***";:  FDR  1=1  TO  1000:  NEXT:  PRINT 

3KP,  EL$;:GOTO  200 

240  RETURN 

250  DN$=INKEY*: IF  DN$=""  THEN  250  ELSE  I=VAL(DN$) 

260  IF  I<0  OR  I>3  THEN  250  ELSE  RETURN 

lOOO  REM  :*  *  *   NEW  FILE  ROUTINE  *  « 

1010  CLS  :  INPUT  "How  oiany  rows";RT 

1020  INPUT  "How  many  columns" ;CT 

1030  BT=1 

1040  GOSUB  14000  :  GOSUB  12000  :  N=l  :  CR=1  :  GOSUB  6040 

1050  CLS: GOSUB  15050: RUN 

1100  GOSUB  1500  :  CR=1  :  GOSUB  6040 

1110  GOSUB  15000  :  RUN 

1200  GOSUB  5000  :  RUN 

1300  GOSUB  1500  :  GOSUB  4000  :  RUN 

1500  REM  t    t    t    READ  DATA  ARRAY  t    t 

1510   DL$="/" 

1520  CLS: PRINT  "Disk  Drive  Containing  Files  (0-3)?" 

1530  GOSUB  250 

1540  SL= 17 15: GOSUB  2400: IF  SL=0  THEN  GOTO  1520 

1550  GOSUB  2800  :  PRINT  "" 

1560  PRINT: IF  10  =  2  PRINT  "Enter  the  File  To  Be  Edited"; 

ELSE  PRINT  "Enter  the  File  To  Be  Printed"; 


1570  INPUT  NV 

1580  IF  NV>FO  THEN  PRINT: PRINT"  t    t    t     INVALID  FILE  #  « 

*":GOTO  1560 

1590     F1$=F$(NV) 

1600  OPEN  "I", 1,F1$ 

1610  INPUT  #1,RT 

1620  INPUT  #1,CT 

1630  INPUT  #1,BT 

1640  CLOSE  1 

1650  GOSUB  12000  :  PRINT 

1660  OPEN  "I", 1,F1$ 

1670  INPUT  #1,RT: INPUT  #1,CT: INPUT  #1,BT 

1680  FOR  1=1  TO  RT 

1690  FOR  J=l  TO  CT 

1700  INPUT  ttl,X(I,J) 

1710  NEXT  J 

1720  NEXT  I 

1730  CLOSE  1 

1740  RETURN 

2000  REM  *  *  *  CHECK  FOR  FILE  DIRECTORY  *  * 

2010  IF  F0=0  THEN  RETURN 

2020  P2=0 

2030  FOR  Pl=l  TO  FO 

2040      OPEN  "I", 1,F$(P1-P2) 

205O      CLOSE  1 

2060  NEXT  PI 

2070  IF  P2=0  THEN  RETURN  ELSE  GOSUB  3400  :  RETURN 

2080  'SUB  TO  CORRECT  DIRECTORY 

2090  IF  P1-P2=F0  THEN  2110 

2100     FOR  I=P1-P2  TO  FO-1  :  F$ (I ) =F* ( I+l >  :  NEXT  I 

2110  F0=F0-1  :  P2=P2+1  :  RETURN 

2400  REM  t     t     t    READ  DIRECTORY  OF  FILES  «  t 

2410  CLS 

2420  F$="CTABLE/DIR: "+DN* 

2430  OPEN  " I " , 1 , F$ 

2440  INPUT  #1,F0 

2450  IF  F0=0  THEN  2490 

2460  FOR  1=1  TO  FO 

2470       INPUT  #1,F*(I) 

2480  NEXT  I 

2490  CLOSE  1  :  GOSUB  2000 

2500  CLOSE  1 

2510  RETURN 

2SOO  REM   *  *  *   DISPLAY  DIRECTORY  t     * 

2810  CLS 

2820  IF  F0=0  THEN  2910 

2830  FOR  1=1  TO  FO  STEP  4 

2840      FOR  J=l  TO  4 

2850  IF  I+J-1>F0  THEN  2910 

2860  IF  J=l  AND  LEN (F$ ( I+J-1 ) )  C 12  AND  I>i  THEN 

PRINT  "" 
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2870         PRINT  TAB(16*(J-1) ) ; " (";STR$(I+J-1) ; ") 

";F*(I+J-1) ; 

2830      NEXT  J 

2890  NEXT  I 

2900  PRINT 

2910  RETURN 

3200  REM   *  *  *  CREATE  "CTABLE/DIR"  *  * 

3210  DN$="0" 

3220  CLOSE  1 

3230  F*= "CTABLE/DIR: "+DN« 

3240  OPEN  "O", 1,F$ 

3250  PRINT  #1,0 

3260  CLOSE  1 

3270  RETURN 

3400  REM   *  *  *   WRITE  NEW  FILES  IN  DIRECTORY   *  * 

3410  F*= "CTABLE/DIR: "+DN$ 

3420  OPEN  "O", 1,F$ 

3430  PRINT  #1,F0 

3440  FDR  1=1  TO  FO 

3450      PRINT  #1,F$(I)  ;  ",  ••; 

3460  NEXT  I 

3470  CLOSE  1 

3480  RETURN 

3600  REM  t    t     t    CHECK  IF  FILENAME  WAS  ALREADY  USED  *  * 

3610  EX=0 

3620  FOR  1=1  TO  FO 

3630      IF  F*=F*(I)  THEN  EX=1  :  RETURN 

3640  NEXT  I 

3650  RETURN 

4000  REM  *  *  «  PRINT  ROUTINE  *  t 

4010  CLS  :  INPUT  "Input  Printer  Width ";PW 

4020  PF=0  :  ML=-1  :  DL$=" ! "  :  TL=-1 

4030  FOR  1=1  TO  N 

4040      FOR  J=l  TO  CT 

4050  I F  LEN ( STRS ( X ( I , J ) ) )  >  ML  THEN 

ML=LEN ( STR$ ( X ( I , J ) ) ) 

4060      NEXT  J 

4070  NEXT  I 

4080  ML=ML+1 

4090  CLS: PRINT: PRINT  "Enter  Column  Headers.   Separating 

lines  with" 

4100  PRINT"the  delimiter  ' !'  to  produce  a  two  column 

header" 

4110  FDR  1=1  TO  CT 

4120  PRINT:PRINT  "Enter  the  Column  Header  -for  Column";  I 

4130  INPUT  M$(I) 

4140  NEXT  I 

4150  MS=-1 

4160  FOR  1=1  TO  CT 

4170      P1=0 

4180      P2=INSTR(P1+1,M$(I) ,DL$) 
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4190      L=P2-P1-1 

4200       IF  P2=0  THEN  L=LEN (M$ ( I ) ) -PI 

4210       IF  L>MS  THEN  MS=L 

4220       IF  P2=0  THEN  4250 

4230         P1=P2 

4240  GOTO  4180 

4250  NEXT  I 

4260  IF  ML>MS  THEN  MX=ML  ELSE  MX=MS 

4270  FOR  1=1  TO  CT 

4280      P1=0  :  NL=0  :  T$="" 

4290      P2=INSTR(P1+1,M*(I) ,DL$) 

4300      L=P2-P1-1 

4310       IF  P2=0  THEN  L=LEN (M$ ( I ) ) -PI 

4320      SF=INT( (MX-L) /2)  :  SB=MX-L-BF 

4330      T$=T$+STRING$(SF, " 

")+MID$(M$(I) ,P1+1,L)+STRING$(SB, "  ") 

4340      NL=NL+1  :  IF  NL>TL  THEN  TL=NL 

4350      IF  P2=0  THEN  M$ ( I ) =T$  :  GOTO  4380 

4360  P1=P2 

4370  GOTO  4290 

4380  NEXT  I 

4390  PR I NT: PR I NT  "Would  you  like  to  number  the  rows  (Y/N)7 

I*  ■ 

4400  A*=INKEY$:IF  A*=""  THEN  4400  ELSE  PRINT  A$ 

4410  IF  A*="Y"  THEN  PI=5  ELSE  PI=0 

4420  MS=INT( <PW-MX  *  CT-PI)/CT) 

4430  PRINT:  PRINT"The  maicimum  #  o-f  spaces  that  can  be  used 

between  columns  is";MS;"." 

4440  PRINT:  INPUT"Enter  the  number  o-f  spaces  between 

columns" ; SP 

4450  IF  SP<0  THEN  INPUT  "TO  RENAME  COLUMNS  TYPE  Y  ELSE 

N";T$  :  IF  T$="Y"  THEN  4020  ELSE  RETURN 

4460  IF  SP<0  OR  SP>MS  THEN  RETURN 

4470  PRINT: INPUT  "Enter  Report  Title" ;T$ 

4480  Pl=(CT-l)*(MX-»-SP)+PI+MX-l 

4490  Pl=INT(Pl/2)  :  P2=PW-P1  :  IF  LEN(T*)>PW  THEN 

T$=LEFT$(T$,PW) 

4500  I F  LEN ( T$ ) >P2  THEN  P 1 = I NT ( ( PW-LEN ( T$ ) ) /2 ) 

4510  LPRINT  TAB(P1);T$ 

4520  LPRINT  TAB(Pl);  :  FOR  1=1  TO  LEN(T$)  :  LPRINT 

"-";:NEXT  I  :  LPRINT 

4530  LPRINT  ""  :  LPRINT  "" 

4540  FOR  J=l  TO  TL 

4550      FOR  1=1  TO  CT 

4560  LPRINT 

TAB( (I-l) « (MX+SP)+PI) ;MID*(M*(I) , (J-1) X+1,MX) ; 

4570      NEXT  I 

4580      LPRINT  "" 

4590  NEXT  J 

4600  A$=""  :  FOR  1=1  TO  MX  :  A*=A$+"-"  :  NEXT  I 

4610  FOR  1=1  TO  CT 
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4620      LPRINT  TAB ( ( I-l ) * (MX+SP) +PI ) ; A*; 

4630  NEXT  I 

4640  LPRINT  ""  :  K=PI+MX-1-INT ( (MX-ML) /2) 

4650  FOR  1=1  TO  RT 

4660      IF  PI=5  THEN  LPRINT  STR* (I ) ; " : " ; 

4670      FOR  J=l  TO  CT 

4680  LPRINT 

TAB( (J-1)*(MX+SP)-LEN(STR$<X(I, J) ) )+K) ;X(I, J) ; 

4690      NEXT  J 

4700      LPRINT  "" 

4710  NEXT  I 

4720  RETURN 

5000  REM   *  :*  *   PURGE  DATA  FILE  DIRECTORY  t    t 

5010  CLS 

5020  PRINT  "Disk  Drive  Containing  Files  to  Be  Purged 

(0-3)?" 

5030  GOSUB  250 

5040  SL=5010  :  GOSUB  2400  :  IF  SL=0  THEN  GOTO  5010 

5050  K=0 

5060  PRINT  TAB  (20);"  ««**;*:«  PURGE  DATA  FILES  ***«" 

5070  FOR  1=1  TO  FO 

5080  A*="":PRINT  F$(I-K);"      KILL  IT?  (Y/N/Q)  "; 

5090  ZI=INKEY*:IF  ZI=""  THEN  5090 

5100  IF  ZI="Q"  THEN  PRINT  Z I : GOTO  5130 

5110  IF  ZI="Y"  THEN  PRINT  Z I : GOSUB  5150  ELSE  PRINT  "N" 

5120  NEXT  I 

5130  GOSUB  3400 

5140  RETURN 

5150  KILL  F$(I-K) 

5160  IF  I-K=FO  THEN  5180 

5170     FOR  J=I-K  TO  FO-1  :  F* ( J) =F$ ( J+1 )  :  NEXT  J: 

5180  F0=F0-1  :  K=K+1 

5190  RETURN 

6000  REM  ttt    EXAMINE  DATA  *t 

6010  INPUT  CR 

6020  IF  CR  <  1  THEN  CR=1 

6030  IF  CR  >  N  THEN  CR=N 

6040  E$=CHR*(31)  :  N*=""  :P1=0:  CLS 

6050  C1=1:R1=1:B1=1 

6060  CR=R1:CC=C1 

6070  GOSUB  6760 

6080  GOSUB  6240  :  GOSUB  6290 

6090  P1=0 

6100  A$=INKEY$  :  IF  A*=""  THEN  6100 

6110  IF  A$=CHR$<91)  AND  P1=0  THEN  N*="":  GOSUB  6360  :  GOTO 

6100 

6120  IF  A$=CHR$(27)  AND  P1=0  THEN  N$=" ": CR=1 : GOSUB 

6380:G0T0  6100 

6130  IF  A*=CHR*(10>  AND  P1=0  THEN  N*="":  GOSUB  6410  :  GOTO 

6100 

6140  IF  A$=CHR$(26)  AND  P1=0  THEN  N$=" ": CR=RT: GOSUB 
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6440: 

6150 

6100 

6160 

6100 

6170 

6500: 

6180 

6560: 

6190 

6720 

6200 

6210 

6220 

6225 

6230 

6240 

6250 

6260 

6270 

6280 

6290 

6300 

6310 

6320 

6330 

6340 

6350 

6360 

6370 

6380 

6390 

6400 

64  lO 

6420 

6430 

6440 

6450 

6460 

6470 

6480 

6490 

6500 

6510 

6520 

6530 

6540 

6550 

6560 

6570 

6580 


GOTO  6100 

IF  A$=CHR*<8)  AND  P1=0   THEN  N*= 

IF  A$=CHR*(9)  AND  P1=0  THEN  N*=" 


•":  GOSUB  6470 
':  GOSUB  6530  : 


GOTO 
GOTO 


•":CC=1: GOSUB 
•":CC=CT: GOSUB 


IF  A$=CHR*(24)  AND  P1=0  THEN  N*=' 

GOTO  6100 

IF  A$=CHR*(25)  AND  P1=0  THEN  N*=' 

GOTO  6100 

IF  (A$>"/"  AND  A*<":")   OR  A*="E"  THEN  Pi=l  :  GOSUB 

:  GOTO  6100 

IF  A*=CHR$(13)  THEN  P1=0:  GOSUB  6590  :  GOTO  6100 

IF  A$=CHR$(8)  AND  Pl=l   THEN  GOSUB  6740  :  GOTO  6100 

IF  A*="Q"  AND  P1=0  THEN  RETURN 

IF  A«=CHR$(31)  AND  P1=0  THEN  RUN 

GOTO  6100 

REM   *  *  *   CLEAR  ENTRY  AREA   *  * 

CC=1 

PRINT  32, " 

PRINT  3  2,CHR*(95) 

RETURN 

REM   *  «  *   BRACKET  ENTRY  *  « 

K=320+12S* (CR-R1)+11* (CC-Cl)+5 

TY=15360+K 

POKE  TY,183  :  POKE  TY+1 1,187 

RETURN 

POKE  TY, 128  :  POKE  TY+1 1,128 

RETURN 

REM  t    t    t       PROCESS  UP  ARROW  t     t 

CR=CR-1:IF  CR<1  THEN  CR=1:G0T0  6400 

IF  CR<R1  THEN  R1=CR: RL=R 1+4: GOSUB  6760 

GOSUB  6340: GOSUB  6290 

RETURN 

REM  t    t    t       PROCESS  DOWN  ARROW  t     t 

CR=CR+1 

IF  CR>RT  THEN  CR=RT:GOTO  6460 

IF  CR>RL  THEN  RL=CR: Rl=RL-4: GOSUB  6340: GOSUB  6760 

GOSUB  6340: GOSUB  6290 

RETURN 

REM   «  *  «   PROCESS  LEFT  ARROW  *  * 

CC=CC-1 

IF  CC<1  THEN  CC=1:G0T0  6520 

IF  CC<C1  THEN  C1=CC:CL=C 1+4: GOSUB  6760 

GOSUB  6340: GOSUB  6290 

RETURN 

REM   *  *  *   PROCESS  RIGHT  ARROW   *  * 

CC=CC+1 

IF  CC  >  CT  THEN  GOSUB  7000 

IF  COCL  THEN  CL=CC:  Cl=CL-4:  GOSUB  6760 

GOSUB  6340: GOSUB  6290 

RETURN 
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CHANGE  VALUE  IN  FIELD   *  * 


RETURN 


t    t 
PRINT  3  2,N$; 


PRINT 


PRINT  3  2,  ■• 


ELSE 


6590  REM   *  *  * 

6600  GOSUB  6260 

6610  IF  N*=""  THEN  GOSUB  6530 

6620  X(CR,CC)=VAL(N*) 

6630  N*="" 

6640  K=326+12S*<CR-R1)+11*(CC-C1) 

6650  IF  X(CR,CC) >999999  THEN  GOTO  6690 

6660  PR  I  NT  5)K+ 1 ,  CHR*  (128); 

6670  PRINT  USING  Z$;X(CR,CC); 

6680  RETURN 

6690  K=K-2:PRINT  5)K,  CHR$  ( 128)  ; 

6700  PRINT  USING  ZZ;X(CR,CC); 

6710  RETURN 

6720  REM  *    t     t  PROCESS  NUMBER  ENTRY 

6730  N$=N$-»-A$  :  PRINT  92,"  "  ; 

CHR*(95)  :  RETURN 

6740  IF  LEN(N$)>0  THEN  N$=LEFT$ (N$, LEN (N*) -1) 

":  PRINT  3  2,N$ 
6750  RETURN 

6760  PRINT  3196, CHR$( 143)  :  IF  CT<Cl+5  THEN  CL=CT 
CL=Cl+4 

6770  IF  RT<Rl+4  THEN  RL  =  RT  ELSE  RL=Rl+4 
6780  FOR  I=C1-1  TO  CL-1 

6790  PRINT  3197+1 1* ( I+l-Cl ) ,STRING$ (4, CHR$( 143) ); :PRINT 
CHR* ( 123) ;: PRINT  CHR* (65+1 );: PRINT  CHR* ( 128) ;: PRINT 
STRINGS (4, CHR* ( 143) ) 
6800  NEXT  I 

6810  PRINT  3195, CHR*( 170) : PRINT  3259, CHR* ( 170) 
6820  FOR  I=R1-1  TO  RL-1 
6830  IF  I  >  8  THEN  K1=0  ELSE  Kl=l 

6840  PRINT  3320+128* (I+1-R1)+K1,STR«(I+1) : rPRINT  CHR$(170) 
6850  FOR  J=C1-1  TO  CL-1 

6860  PRINT  3327+11* (J+1-C1)+128* (I+l-Rl) ,CHR*(128) ; :PRINT 
USINGZ$; X (I+l, J+1) 
6870  NEXT  J 

6880  PRINT  3387+128* (I+l-Rl) ,CHR$( 170) 
6890  NEXT  I 
6900  RETURN 

7000  REM  *  *  *   PROCESS  LAST  COLUMN  ENTRY   *  * 
7010  CC=1  :  C1=CC  :  CL=Cl+4 
7020  CR=CR+1  :  IF  CR  >  RT  THEN  CR=RT 
7030  IF  CR  >  RL  THEN  RL=CR  :  Rl=RL-4 
7040  GOSUB  6340 

7050  IF  CT  >  5  OR  CR  =  RL  THEN  GOSUB  6760 
7060  RETURN 

11000  REM  *  *  *  VERIFY  FILENAME  -  FW=1  IF  FAILS  *  * 
11010  FW=0 

11020  P1=INSTR(1,S$, "/") 

11030  IF  PI  <>  O  THEN  IF  LEN(S*)-P1>3  THEN  FW=1  :  RETURN 
11040  IF  PI  >8  THEN  FW=1  :  RETURN 
11050  IF  P1=0  THEN  IF  LEN(S$)>8  THEN  FW=1  :  RETURN 
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11060  IF  PI  <>  O  THEN  IF  ASC (MID$ (S$, Pl+1 , 1 ) )<65  CR 
ASC(MID$<S*,P1+1, 1) ) >90  THEN  FW=1  :  RETURN 

11070  IF  ASC(LEFT*(S$, 1) )<65  OR  ASC (LEFT* <S$, 1 ) ) >90  THEN 
FW=1  :  RETURN 

llOSO  IF  P1=0  THEN  P2=LEN(S$)  ELSE  P2=P1-1 
11090  P3=l 

11100  FOR  I=P3  TO  P2 
11110       A*=MID*(S$, I, 1) 

11120       IF  (ASC(A$)>47  AND  ASC(A$)<58)  OR  (ASC<A*)>64 
AND  ASC(A*)<91  )  THEN  11140 
11130  FW=1  :  RETURN 

11140    NEXT  I 

11150  IF  P3=l  AND  PI  <>  O  THEN  P3=P1+1  :  P2=LEN(S$)  :  GOTO 
11100 

11160  RETURN 

12000  REM  *  *  *  DIMENSION  ARRAYS  *  t 

12010  NV=MEM-500  :  NV=INT (NV/ (4* (RT  «  CT  «  BT+1)))-1 
12020  IF  MM  >  NV  THEN  PRINT  "INSUFFICIENT  SPACE  TO  READ  IN. 
": INPUT  "PRESS  ENTER  TO  CONTINUE" ;A$ 
12030  DIM  X(RT,CT) 
12040  MN=NV 
12050  RETURN 

13000  REM  *  *  *    RENAME  FILES    :*  * 
13010  CLS 

13020  PRINT  "Disk  Drive  To  Contain  Files  (0-3)?" 
13030  GOSUB  250 

13040  SL=S240: GOSUB  2400  :  GOSUB  2800  :   PRINT  "" 
13050      INPUT  "New  Name";S*  :  GOSUB  11000 
13060       IF  FW=1  THEN  PRINT  "INVALID  FILENAME"  :  GOTO 
13050 

13070         F$=S$  :  GOSUB  3600 

130SO         IF  EX=1  THEN  PRINT: PRINT  "A  File  Currently  Uses 
the  New  Fi lename. ": PRINT  "Do  You  Want  to  Use  it  Anyway 
(Y/N)?"  ELSE  GOTO  13100 

13090  A*=INKEY$:IF  A$=""  THEN  13090  ELSE  PRINT  A$:  IF  A$<> 
"Y"  THEN  13050 
13100  F1$=F$ 

13110  RETURN 

14000  REM  *  *  «    NAME  FILE     t    t 
14010  CLS 

14020  PRINT  "Disk  Drive  To  Contain  File  (0-3) ?": GOSUB  250 
14030  SL=8330  :  GOSUB  2400  :  GOSUB  2800  :  PRINT  "" 
14040  PRINT 

14050  PRINT  "FILENAME  FOR  TABLE";  :  INPUT  S$  :  GOSUB  11000 
14060  IF  FW=1  THEN  PRINT  "INVALID  FILENAME"  :  GOTO  14050 
14070     F$=S$  :  GOSUB  3600 

14080  IF  EX=1  THEN  PRINT  "FILENAME  IS  ALREADY  USED.  TYPE 
Y  TO  USE  NAME  ANYWAY  ELSE  N";:  INPUT  A$  :  IF  A$  <>  "Y"  THEN 
14050 

14090     F1$=F$ 
14100  RETURN 
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15000 

15010 

15020 

15030 

15040 

15050 

15060 

15070 

15080 

15090 

15100 

15110 

15120 

15130 

15140 

15150 

15170 

15160 

15170 

1 9000 

19010 

19020 

19030 

19040 

19050 

19060 

FILES 

SL=0  : 

1 9070 


*  *  *  EXITING  ROUTINE  FOR  OPTIONS  1  & 


REM 
CLS 

PRINT  "Do 
ZI=INKEYi 
IF  ZI="Y" 
SL=a550  : 


«  * 


TO  DISK  — > 


1,F1*+":  '■+DN* 


You  Wish  To  Rename  File 
IF  21=""  GOTO  15030 
THEN  GOSUB  13000 
GOSUB   2400  :  GOSUB  2S00 
PRINT: PRINT 

PRINT  "WRITING 
OPEN  "O" 

PRINT  #i,RT: PRINT  #1 , CT: PRINT#1 , BT 
FOR  1=1  TO  RT 
FOR  J=l  TO  CT 
PRINT  #1,X(I, J) ; 
NEXT  J 
NEXT  I 

CLOSE  1  :  F*=F1$  :  GOSUB  3600 


(Y/N) 


PRINT 


Fl* 


IF  EX=1  THEN 


F*(FO)=Fi* 


GOSUB  3400 


F0=F0+1 
RETURN 

REM  t    t    *    ERROR  HANDLING  ROUTINE  *  * 
IF  ERR=10  AND  ERL=6590  THEN  N$=' 

ERR  <>  106  THEN  19080 

ERL=5150  THEN  RESUME  5160 

ERL=2040  THEN  GOSUB  2080  :  RESUME  2050 

ERL  <>  2430  THEN  19080 


RESUME  6630 


IF 
IF 
IF 
IF 
IF 
ON 


SL=5010  OR  SL=5410   THEN  PRINT  "NO 
PRESS  ENTER  TO  CONTINUE":  INPUT  A$ 


SL=1715  OR 
THIS  DRIVE. 
RESUME  2500 
IF  SL=S240  OR  SL=8330   OR  SL=8550  THEN  GOSUB  3200 
RESUME  2410 
19080  PRINT  "ERR=":ERR 
RESUME 


"ERL=";ERL  :  ON  ERROR  GOTO  O 
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i=-ROGF<^f^ 


iaJ  I  i_c::cDX  /Brf=*^ 


10  CLEAR  1000 

20  ON  ERROR  GOTO  19000 

30  DEFINT  P:DEFSTR  2:DEFDBL 


F,L 


40  DIM  S(1500, 1) ,F*(95) ,P(1500) ,X(1500) ,B<50) 

50  LB$=CHR* ( 1 23 ) : RB$=CHR$ ( 1 25 ) : MB*=CHR$ ( 1 76 ) : BL$=CHR$  <  252 ) 

60  REM  «  *  *   WILC0X30N  SIGNED  RANKS  TEST   *  * 

70  CLSzPRINT  5>20,  "WILCOXSON  SIGNED  RANKS  TEST" 

75  PRINT: PRINT: PRINT  "Disk  Drive  Containing  Data  Files 

(0-3)?" 

SO  DN$=INKEY$: IF  DN$=""  THEN  SO  ELSE  I=VAL(DN$) 

85  IF  KO  OR  I>3  THEN  75 

90  SL=100  :  GOSUB  11000  :  IF  SL=0  THEN  70  ELSE  GOSUB  12000 

:  PRINT  "" 

100  PRINT: INPUT  "Enter  the  Two  Data  Files  (e.g.  3,5)";F1,F2 

IF  F1>F0  OR  F2>F0  OR  FKl  OR  F2^ 

OPEN  "I", 1,F$(F1)+": "+DN$ 

INPUT  #1,N 

N  :  INPUT  #1,S(I 


110 
120 
130 
140 
150 
160 
170 
ISO 
190 
200 
210 
You 


TO 


O) 


Fi les  (e. g, 
:i  THEN  100 


NEXT  I 


F$(F2)+": "+DN$ 


FOR  1=1 

CLOSE  1 

OPEN  " I " , 1 . 

INPUT  #1,M 

FOR  1=1  TO  M  :  INPUT  #1,S(I,1)  :  NEXT  I 

CLOSE  1 

IF  N=M  GOTO  250 

PRINT:  PRINT  "The  Data  Sets  are    NOT  o-f  Equal  Size.   Do 

Want  to  Continue  using  the  Smaller  #  o-f  data  points. 
Enter  other  Data  Files,  or  Return  to  the  Master  Menu 
(C/E/R)?"; 

220 

GOTO  240: ELSE  IF  2I="E"  THEN 

THEN  PRINT  "R" : RUN"MENU/BAS" 


220  ZI=INKEY*:IF  ZI=""  THEN 
230  IF  2I="C"  THEN  PRINT  21 
PRINT  "E":RUN  ELSE  IF  21 "R" 
ELSE  GOTO  220 
240  IF  N  >  M  THEN  N=M 

N 

=S(I,0)-S(I,  1) 
1)=0  THEN  J=J-1 


250 

J=l 

260 

FOR 

1  =  1 

TO 

270 

S(J. 

■  1)^ 

280 

IF  S(J 

290 

J=J  +  1 

300 

NEXT  I 

310 

N=J- 

-1 

320 

FOR 

1  =  1 

TO 

330 

M=l 

:  N; 

1=N 

10000 

340 

JS=] 

L 

350 

FOR 

K=2 

TO 

>60 


N  :  X(I)=ABS(S(I, 1) )  :  NEXT  I 
:  FOR  1=1  TO  Nl   :  P(I)=I  :  NE: 


N 
R=ABS(X (P(JS) ) /X (P(K) ) ) 


GOSUB 
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370       IF  R  <  .999  OR  R  >  1.001  THEN  JE=K-1  :  GOSUB  14000 

:  JS=K 

380  NEXT  K 

390  JE=N  :  GOSUB  14000 

400  TP=0:TN=0 

410  FOR  1=1  TO  N 

420      IF  S(I,1)  >  O  THEN  TP=TP+S ( I , 1 )  ELSE 

TN=TN+ABS(S(I, 1) ) 

430  NEXT  I 

440  IF  TP  >  TN  THEN  TS=TN  ELSE  TS=TP 

450  CLS:PRINT3>20,  "WILCOXSON  SIGNED  RANKS  TEST"  :  PRINT:  PRINT 

460  PRINT  "HYPOTHESIS  NULL 

ALTERNATIVE" 

470  PRINT  "     (1)  E(X)   =   E(Y)  E(X)  =/= 

E':Y)  " 

480  PRINT  "     (2)  E(X)  <=   E<Y)  E<X>    > 

E(Y)  " 

490  PRINT  "     (3)  E<X)  =>   E<Y)  ECX)   < 

E(Y)  " 

500  PRINT:  PRINT  "Where  X  and  Y  represent  Data  Sets  1  Z<    2 

Resp  ectively" 

510"  JP=640:!<S=3:  GOSUB  8000 

520  IF  10=1  V=TS  ELSE  IF  10=2  V=TP  ELSE  IF  10=3  V=TN 

530  CLS:PRINT320, "WILCOXSON  SIGNED  RANKS  TEST" : PRINT: PRINT 

540  PRINT  "Test  Statistic  =  ";V 

545  PRINT  "Sample  Size  =  ";N 

550  GOSUB  2000 

560  ON  I O  GOSUB  600 , 800 , 1 000 

570  PRINT  3960, "Press  SPACE  BAR  when  ready  to  return  to 

main  menu"; 

580  IF  INKEY$=  ""  GOTO  580  ELSE  RUN  "MENU/BAS" 

600  REM  t    t    t    NULL  HYPOTHESIS  #1  *  * 

610  PRINT  "P-ValLie  =  ";:  PRINT  USING  "#.####":  2  *  FOX 

620  PRINT: INPUT"Enter  Desired  Alpha  Level ";AL 

630  IF  AL<0  OR  AL>1  THEN  PRINT  "***  ERROR  ttt    ENTER  A  VLAUE 

BETWEEN  O  &  1 " : GOTO  490 

640  PRINT: PRINT  "Far  Alpha  =";AL;",  "; 

650  IF  AL<2  *  FOX  THEN  PRINT  "accept":  ELSE  PRINT  "reject"; 

660  PRINT  "  the  hypothesis  that  E(X)  =  E(Y)" 

670  PRINT  "-VS-  the  alternative  that  E  •:  X )  =/=  E(Y)." 

680  RETURN 

800  REM  *  *  *  NULL  HYPOTHESIS  #2  *  * 

810  PRINT  "P-Value  =  ";:PRINT  USING  "#.####"; FOX 

820  PRINT: INPUT"Enter  Desired  Alpha  Level ";AL 

830  IF  AL<0  OR  AL>1  THEN  PRINT  "tt*    ERROR  ***  ENTER  A  VLAUE 

BETWEEN  O  Zt     1  "  :  GOTO  620 

840  PRINT: PRINT  "For  Alpha  =";AL;",  "; 

850  IF  AL<FaX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 

860  PRINT  "  the  hypothesis  that  E(X)  <=  E<Y)" 

870  PRINT  "-VS-  the  alternative  that  E<X)  >  E(Y)." 

880  RETURN 
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iOOO  REM  «  «  *  NULL  HYPOTHESIS  #3  *  * 

1010  PRINT  "P-Value  =  ";:  PRINT  USING  "#.####'■:  FOX 

1020  PRINT: INPUT"Enter  Desired  Alpha  Level ■; AL 

1030  IF  AL<0  OR  AL>1  THEN  PRINT  ":♦**  ERROR  ***  ENTER  A 

VLAUE  BETWEEN  O  &  1 " :  GOTO  L'ZO 

1040  PRINT: PRINT  "For  Alpha  =";AL;",  "; 

1050  IF  AL<FOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 

1060  PKTNT  "  the  hypothesis  that  E(X)  =>  E(Y)" 

1070  PRINT  "-VS-  the  alternative  that  E(X)  <  E(Y)." 

1080  RETURN 

2000  REM  *  :«  «   WILCOX  SON  DISTRIBUTION  *  * 

2020  PX=V: IN=N:FLAG=0 

2030  IF  V>IN* (IN+l>/2  OR  V<0  OR  IN<1  THEN  PRINT  "  *  *  ERROR 

*  «":GOTO  100 

2040  I F  V > I N*  ( I N+ 1 )  /4  THEN  V=  I N*  ( I N+ 1 )  /2-v'- 1 :  FLAG=  1 

2050  IF  V>27  OR  <IN>9  AND  V>IN)  THEN  GOSUB  2600  ELSE  GOSUB 

2070 

2060  RETURN 

2070  IF  INT(V)=V  GOSUB  2110  ELSE  GOTO  2090 

2080  RETURN 

2090  V=INT(V) : GOSUB  21 10: TFOX=FOX: V=V+1 : GOSUB 

21 10: FOX= (FQX+TFOX ) /2: Q0X=1-F0X 

2100  RETURN 

2110  II (2)=0: II (3)=0: II (4)=0: II (5)=0: II (6)=0 

2120  C=l 

2130  N=INT(.5« (-1+SQR( l+a*V) ) ) 

2140  FOR  1=1  TO  N 

2150  I1(1)=0 

2160  II  (1)=I1  ■:  1)+1 

2170  IF  I<2  GOTO  2320 

2130  I1(2)=I1(1) 

2190  II (2)=I1 (2)+l 

2200  IF  I<3  GOTO  2320 

2210  II (3)=I1 (2) 

2220  II (3)=I1 (3)+l 

2230  IF  I<4  GOTO  2320 

2240  II (4) =11 (3) 

2250  II (4) =11 (4)+l 

2260  IF  I<5  GOTO  2320 

2270  II '5) =11 (4) 

2280  II (5) =11 (5)+l 

2290  IF  I<6  GOTO  2320 

2300  II (6>=I1 <5) 

23 1 0  1 1  ( 6 )  =  1 1  ( 6 )  + 1 

2320  S=I1(1)+I1(2)+I1(3)+I1(4)+I1(5)+I1(6) 

2340  IF  S>V  THEN  GOTO  2360 

2350  C=C+1:0N  I  GOTO  2550,2540,2530,2510,2510,2500 

2360  IF  I>1  GOTO  2380 

2370  I i ( 1 ) = I N : GOTO  2550 

2380  FOR  IL=1  TO  I-l 

2390  S=0:J=I-IL 
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2400  IF  J-1<=0  THEN  GOTO  2440 

2410  FOR  K=l  TO  J-1 

2420  S=S+I1 CK) 

2430  NEXT  K 

2440  S=S+(IL+1) 1 (J)+(IL+l)*(IL+2)/2 

2450  IF  S>V  THEN  GOTO  2480 

2460  I1(I+1-IL)=IN 

2470  ON  J  GOTO  2540,2530,2520,2510,2500 

2480  NEXT  IL 

2490  II (1)=IN:G0T0  2550 

2500  IF  I1(6)<IN  GOTO  2310 

2510  IF  I1(5)<IN  GOTO  2280 

2520  IF  Ii(4)<IN  GOTO  2250 

2530  IF  I1(3)<IN  GOTO  2220 

2540  IF  I1(2)<IN  GOTO  2190 

2550  IF  IldXIN  GOTO  2160 

2560  NEXT  I 

2570  F0X=C/2CIN:a0X=l-F0X 

2575  IF  FLAG=1  THEN  FOX=QOX : Q0X=l-FaX 

2580  RETURN 

2590  REM  *  «  »  WILCOX  SON  SIGNED  RANK  APPROX  D I  ST  I  BUT  I  ON  :* 

2600  MU=IN* (IN+1) /4 

2610  VAR=(2N+1)  :*  (IN+l)N/24 

2620  NX=V+0.5 

2630  GOSUB  7010 

2640  L4=(3NC2+3N-1) / (ION* (IN+1) * (2N+1) ) 

2650  W=(NX-hU) /SDR(VAR) : F= (EXP (-0. 5*WC2) ) /SQR (6. 28313) 

2660  F3=F* (WC3-3*W) 

2670  F0X=FDX+L4*F3 

2680  F5=-F*  (WC5-10*Wi:3+15*W) 

2690  F7=-F«(WC7-21*WC5H-105*WC3-105*W) 

2700 

L6=4* (3NC4+6NC3-3N+1) / (35* (IN* (IN+1) * (2N+1) ) C2) 

2710  F0X=F0X+L6«F5+-5*F7*L4C2 

2720  IF  FLAG=1  THEN  F0X=1-F0X 

2730  QOX=l-FaX 

2740  RETURN 

7000  REM   *  «  *   NORMAL  DISTRIBUTION   *  * 

7010  SD=SQR(VAR) 

7020  XN=(NX-MU)/SD:AX=ABS(XN) 

7030  X2=.39S9423*EXP( (-AX*AX)/2) 

7040  NT=1/ (1+.2316419*AX) 

7050 

Q0X  =  X2* ( ( . 3193815*NT) - ( . 3565638*NTC2)  +  ( 1 . 78147S*NTC3) - ( 1  -  321 

256*NTC4) + ( 1 . 330274*NTC5) ) 

7060  IF  XN<0  GOTO  7090 

7070  F0X=1-Q0X 

7080  RETURN 

7090  FDX=QOX 

7100  QOX=l-FaX 

7110  RETURN 
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8000  REM  «  *  «   OPTION  SELECTION  ROUTINE  * 

8010  PRINT  3JP+5, "NUMBER  OF  OPTION  DESIRED 

" ; LB* ; CHR* (1 76 )  ;  RB* 

8020  ZI=INKEY$:IF  ZI=""  THEN  S020 

8030  PRINT  3JP+31,ZI; : I0=VAL(2I) 

8040  IF  I0<1  OR  IO>KS  THEN  PRINT  3JP,"*:** 


OPTION 

10000  REM 

10010 

10020 

10030 

10040 

SPECIAL 

10050 

1 0060 

1 0070 

10080 

10090 

101 00 

10110 

10120 

10130 

10140 

101 50 

10160 

ELEMENT 

10170 

10180 

1 0 1 90 

10200 

LARGE 

10210 

10220 

1 0230 

10240 

10250 

10260 

LARGE 

1 0270 

10280 

1 0290 

11000 

11010 

11020 

11030 

11040 

11050 

1 1 060 

11070 

1 1 080 

11090 


»**"; 


NOT  A  VALID 
TO  lOOO: NEXT: PRINT  3JP,BL$:G0T0  8010 
SORT  «  :* 


FOR  1=1 

*  «  ARRAY 

L=l  :  B(L)=NH-1 

J=B(L)  'SET  END  OF  ARRAY  SEGMENT 
I=M-1   'SET  START  OF  ARRAY  SEGMENT 

IF  J-M  <  3  THEN  10230    'HANDLE  1  OR  2  ELEMENTS  IN 
CASE 

Mi  =  INT':  (I+J) /2)       'SET  COMPARE  ELEMENTS 
REM  FIND  A  LARGE  ELEMENT  AMONG  THE  SMALL  ONES 
1  =  1  +  1 
IF  I=J  THEN  10170 

IF  X<P(I))  <=  X(P(M1>)  THEN  10070 

REM  FIND  A  SMALL  ELEMENT  AMONG  THE  LARGE  ONES 

J=J-1 

IF  I=J  THEN  10170 

IF  X(P<J))  >=  X(P(M1))  THEN  lOllO 

REM  EXCHANGE  OUT  OF  PLACE  ELEMENTS 

H=P (I )  :  P (I ) =P  <  J )  :  P ( J ) =H  :  GOTO  1 0070 

REM  ADJUST  COMPARE  ELEMENT  TO  NEW  MIDDLE 

IF  I  >  Ml  THEN  1=1-1 

IF  J=M1  THEN  10210 

H=P(I)  :  P(I)=P(M1)  :  P<M1)=H 

REM  SAVE  STARTING  POINT  FOR  ARRAY  SEGMENT  OF 
ELEMENTS 

L=L+1  :  B(L)=I 

GOTO  10020 
IF  J-M  <  2  THEN  10270 

IF  X(P(M))  <=  X(P(M+1))  THEN  10260 
H=P(M)  :  P<M)=P(M+1)  :  P(M+1)=H 
REM  SET  BEGIN  AND  ENDPOINTS  FOR  ARRAY  SEGMENT 
ELEMENTS 
M=B(L)+1  :  L=L-1 
IF  L>0  THEN  GOTO  10020 
RETURN 

REM  t    t     t    READ  FILE  DIRECTORY  :*  « 

CLS:PRINT5)20,  "WILCOXSON  SIGNED  RANKS  TEST":PRINT 
F$="EDITDATA/DIR: "+DN$ 
OPEN  " I " ,  1 ,  F$ 
INPUT  #1,F0 
FOR  1=1  TO  FO 

INPUT  #1,F$(I) 
NEXT  I 
CLOSE  1 
RETURN 


OF 
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12000  REM  *  *  «  PRINT  FILE  DIRECTORY  «  t 

12010  CLS:PRINT.920,  "WILCOXSON  SIGNED  RANKS  TEST":  PRINT 

12020  FOR  1=1  TO  FO  STEP  4 

12030      FOR  J=l  TO  4 

12040  IF  I+J-1>F0  THEN  12090 

12050  IF  J=l  AND  LEN (F$ ( I+J-1 ) X 12  AND  I>1  THEN 

PRINT  "" 

12060  PRINT  TAB(16*(J-1) ) ; " (";STR$(I+J-1) ; ") 

";F$<I+J-1) ; 

12070      NEXT  J 

12080  NEXT  I 

12090  RETURN 

14000  REM  t    t    t       AVERAGE  RANKS   *  * 

14010  IF  JS=JE  THEN  S (P ( JS) , 1 ) =SGN (S (P ( JS) , 1 ) )  t    JS  : 

RETURN 

14020  AV=0  :  FOR  I=JS  TO  JE  :  AV=AV+I  :  NEXT  I  : 

AV=AV/ (JE-JS+1) 

14030  FOR  I=JS  TO  JE  :  S (P ( I ) , 1 ) =SGN (S (P ( I ) , 1 ) )  *  AV  :  NEXT 

I 

14040  RETURN 

19000  REM  *  «  «  ERROR  HANDLING  ROUTINE   *  * 

19010  IF  ERR=106  AND  (SL=2005  OR  SL=2410  OR  SL=2710  OR  3110 

OR  SL=3510  OR  SL=3910)  THEN  INPUT  "NO  FILES  ON  THIS  DRIVE. 

PRESS  ENTER  TO  CONTINUE" ; A$  :  SL=0  :  RESUME  llOSO 

19100  PRINT  "ERR=";ERR: "ERL=";ERL  :  ON  ERROR  GOTO  O 


AQ 


F^-ROGIR^r^   S,    Mi!=itNJ!N|/:B«=^S 


10  CLEAR  1000 

20  ON  ERROR  GOTO  19000 

30  DEFINT  PzDEFSTR  2 : DEFDBL  F,L 

40  DIM  F*(95) ,P(3000) ,X(3000) ,B<50) 

50  LB$=CHR$ (123) : RB$=CHR$ (125) : MB$=CHR$ ( 176) : BL$=CHR$ (252) 

60  REM  *  *  *   MANN-WHITNEY  TEST   *  * 

70  CLS:  PRINT  5)20,  "MANN-WHITNEY  TEST" 

75  PR I NT: PR I NT: PR I NT  "Disk  Drive  Containing  Data  Files 

(0-3)?" 

30  DN$=INKEY*: IF  DN$=""  THEN  80  ELSE  I=VAL(DN$) 

35  IF  KO  OR  I>3  THEN  GOTO  75 

90  SL=100  :  GOSUB  11000  :  IF  SL=0  THEN  70  ELSE  SOSUB  12000 

:  PRINT  "" 

100  PRINT: INPUT  "Enter  the  Two  Data  Files  (e.g.  3,5)":Fi,F: 

110  IF  F1>F0  OR  F2>F0  OR  FKl  OR  F2<  1  THEN  100 

120  OPEN  "I", 1,F$(F1)+": "+DN* 

130  INPUT  #1,N 

140  FOR  1=1  TO  N  :  INPUT  #1,X(I)  :  NEXT  I 

150  CLOSE  1 

160  OPEN  "I", 1,F$(F2)+": "+DN$ 

170  INPUT  #1,M 

ISO  FOR  I=N+1  TO  N+M  :  INPUT  #1,X(I)  :  NEXT  I 

190  CLOSE  1 

200  MT=M 

210  N1=N+M:M=1  :  FOR  1=1  TO  Nl   :  P(I)=I  :  NEXT  I  :  GOSUB 

1 0000 

220  JS=1 

230  FOR  K=2  TO  Nl 

240  IF  X ■P(JS) )<>X (P<K) )  THEN  JE=K-i:  GOSUB  14000:  J3=K 

250  NEXT  K 

260  JE=N1  :  GOSUB  14000 

270  TS=0 

230  FOR  1=1  TO  N 

290   TS=TS+X ( I ) 

300  NEXT  I 

305  UX=TS-N* (N+1) /2  :  UY=NT-UX 

310  IF  UX  <  UY  THEN  UU=UX  ELSE  UU=UY 

320  CLS:  PRINT.920,  "MANN-WHITNEY  TEST"  :  PRINT:  PRINT 

330  PRINT  "HYPOTHESIS  NULL 

ALTERNATIVE" 

340  PR  I  NT  "     ( 1 )  E  ( X )   =   E  ( Y )  E  ( X )  =/  = 

E(Y)  " 

350  PRINT  "     (2)  E(X)  <=   E(Y)  E(X)   > 

E(Y)  " 

360  PRINT  "     (3)  E(X)  =>   E(Y)  E(X) 

E(Y)  " 
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370  PRINT:  PRINT  "Where  X  and  Y  represent  Data  Sets  1  St  2 

Respect i  vel y " 

380  JP=640:KS=3:GOSUB  8000 

385  IF  10=1  THEN  U=UU  ELSE  IF  10=2  THEN  U=UX  ELSE  IF  10=3 

THEN  U=UY 

390  CLS : PR I NT320 , " MANN : -WH I TNEY  TEST " : PR I NT : PR I NT 

400  PRINT  "Test  Statistic  =  ";U 

405  PRINT  "Sample  Sizes  are    ";N;: PRINT  "  and  " ; MT 

410  GOSUB  2000 

420  ON  10  GOSUB  600,300,1000 

430  PRINT  5)960,  "Press  SPACE  BAR  when  ready  to  return  to 

main  menu"; 

440  IF  INKEY$=  ""  GOTO  440  ELSE  RUN  "MENU/BAS" 

600  REM  t    t    t    NULL  HYPOTHESIS  #1  «  ;« 

610  PRINT  "P-Value  =  ";:PRINT  USING  "#-####";2  *  FOX 

620  PRINT: INPUT"Enter  Desired  Alpha  Level "; AL 

630  IF  AL<0  OR  AL>1  THEN  PRINT  "*:*«  ERROR  ttt    ENTER  A  VLAUE 

BETWEEN  O  S.  1  "  :  GOTO  620 

640  PRINT:PRINT  "For  Alpha  =":AL;",  "; 

650  IF  AL<2  *  FOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 

660  PRINT  "  the  hypothesis  that  E(X)  =  E(Y)" 

670  PRINT  "-VS-  the  alternative  that  E(X)  =/=  E(Y)." 

680  RETURN 

800  REM  *  *  *  NULL  HYPOTHESIS  #2  *  « 

810  PRINT  "P-Value  =  ";:PRINT  USING  "#.####"; FOX 

820  PRINT: INPUT"Enter  Desired  Alpha  Level";AL 

830  IF  AL<0  OR  ALM  THEN  PRINT  "«*«  ERROR  ***  ENTER  A  VLAUE 

BETWEEN  O  ^    1 " : GOTO  620 

840  PRINT:PRINT  "For  Alpha  =";AL;",  "; 

850  IF  AL<FOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 

860  PRINT  "  the  hypothesis  that  E(X)  <=  E(Y)" 

370  PRINT  "-VS-  the  alternative  that  E(X)  >  E(Y)." 

880  RETURN 

lOOO  REM  *    t    t    NULL  HYPOTHESIS  #3  *  t 

1010  PRINT  "P-Value  =  ";:PRINT  USING  "#.####"; FOX 

1020  PRINT: INPUT"Entsr  Desired  Alpha  Level "; AL 

1030  IF  AL<0  OR  ALM  THEN  PRINT  "***  ERROR  ttt    ENTER  A 

VALUE  BETWEEN  O  &  1 " : GOTO  620 

1040  PRINT: PRINT  "For  Alpha  =";AL;",  "; 

1050  IF  AL<FOX  THEN  PRINT  "accept ";  ELSE  PRINT  "reject"; 

1060  PRINT  "  the  hypothesis  that  E(X)  =>  E(Y)" 

1070  PRINT  "-VS-  the  alternative  that  E':X)  <  E(Y)." 

1080  RETURN 

2000  REM   *  *  *   MANN-WHITNEY  DISTRIBUTION   *  t 

2010  M=MT 

2020  IF  M>N  THEN  M=N:N=MT 

2030  FLAG=0 

2040  IF  U>N/2-.5  THEN  U=N-U-1 : FLAG=1 

2050  I N=N+M : V=TS : I =N 

2060  IF  U>N  OR  U<0  OR  N< 1  OR  M< 1  THEN  PRINT  "  *  *  ERROR 

*":GOTO  100 
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2070  IF  V>27  OR  M>9  THEN  GOSUB  2600  ELSE  GOSUB  2090 

2080  RETURN 

2090  IF  INT(V)=V  THEN  GOSUB  2130  ELSE  GOTO  2110 

2100  RETURN 

2110  V=INT(V) : GOSUB  2130: TFOX=FOX: V=V+1 : GOSUB 

2130: FOX= (FOX+TFOX ) /2: Q0X=1-F0X 

2120  RETURN 

2130  II  (2)  =0:  II  (3)=0:  1 1  (4)  =0:  1 1  (5)  =0:  1 1  (6)  =0 

2140  C=0 

2150  I1(1)=0 

2160  I1(1)=I1(1)+1 

2170  IF  I<2  GOTO  2310 

2180  II (2)=I1 (1) 

2190  II  (2)=I1  (2)-<-l 

2200  IF  I<3  GOTO  2320 

2210  II (3) =11 (2) 

2220  II (3) =11 (3) +1 

2230  IF  I<4  GOTO  2320 

2240  II (4) =11 (3) 

2250  II (4)=I1 (4)+l 

2260  IF  I<5  GOTO  2320 

2270  1 1  ( 5 ) = I  1  <  4 ) 

2280  II (5) =11 (5)+l 

2290  IF  I<6  GOTO  2320 

2300  II (6) =11 (5) 

2310  II (6)=I1 <6)+l 

2320  S=I1(1)+I1(2)+I1(3)+I1(4)+I1(5)+I1(6) 

2330  IF  S>V  THEN  GOTO  2350 

2340  C=C+1:0N  I  GOTO  2540,2530,2520,2510,2500,2490 

2350  IF  I>1  GOTO  2370 

2360  I1(1)=IN:G0T0  2540 

2370  FOR  IL=1  TO  I-l 

2380  S=0:J=I-IL 

2390  IF  J-1<=0  THEN  GOTO  2430 

2400  FOR  K=l  TO  J-1 

2410  S=S+I1 (K) 

2420  NEXT  K 

2430  S=S+<IL+i) 1 (J)+(IL+1) * (IL+2) /2 

2440  IF  S>V  THEN  GOTO  2470 

2450  II (I+1-IL)=IN 

2460  ON  J  GOTO  2530,2520,2510,2500,2490 

2470  NEXT  IL 

2480  II (1>=IN:G0T0  2540 

2490  IF  I1(6)<IN  GOTO  2310 

2500  IF  11(5) 'IN  GOTO  2280 

2510  IF  I1(4)<IN  GOTO  2250 

2520  IF  I1(3)<IN  GOTO  2220 

2530  IF  I1(2)<IN  GOTO  2190 

2540  IF  1 1-  IX  IN  GOTO  2160 

2550  S=l 

2560  FOR  K=0  TO  N-1 : S=S* ( IN-K) / (N-K) : NEXT  K 


t  ^ 


2570  F0X=C/S:Q0X=1-F0X 

2580  IF  FLAG=1  THEN  FOX=QOX: Q0X=1-F0X 

2590  RETURN 

2600  REM  *  *  *  MANN-WHITNEY  APPROXIMATE  DISTRIBUTION  :*  % 

2610  MU=N*M/2: VAR=N*M* (N+M+1) /12:NX=U+.5 

2620  GOSUB  7010 

2630  L4=(MC2+NE2-?-IN+M*N)/(20*M*N*(IN+l)  ) 

2640 

L6=  (2*  (MC4+NC4)  +4*M*N«  (MC2+NC2)  +6*  (MC2)  *  (NC2)  +4:«  (MCS+NCS)  +7 

«M*N« (M+N) +MC2+NC2+2*M*N-IN) / (210* (MC2) * (NC2) * ( ( IN+1 >  C2) ) 

2650  W=  (NX-MU)  /SCR  (VAR)  :  F=  (EXP  (-0.  5*WC2)  )  /SQR  (6.  28318) 

2660  F3=F* (WC3-3*W) 

2670  F5=-F« (WC5-10*WC3+15*W) 

26S0  F7=-F*  (WC7-21.1tWC5+105*WJ:3-105*W) 

2690  F0X=F0X+L4  *  F3 

2700  F0X=F0X+L6  :*  F5 

2710  F0X=F0X+.5  *  F7  *  L4!:2 

2720  IF  FLAG=1  THEN  F0X=1-F0X 

2730  Q0X=1-F0X 

2740  RETURN 

7000  REM  %     %    %       NORMAL  DISTRIBUTION   *  * 

7010  SD=SQR(VAR) 

7020  XN= (NX-MU) /SD:AX=ABS(XN) 

7030  XZ=.3989423*EXP( (-AXtAX) /2) 

7040  NT=1/ (1+.2316419*AX) 

7050 

QOX=XZ*  (  ( .  3193S15.*NT)  -  ( .  356563S:*NTC2)  +  ( 1 .  7Si478:*NTC3)  -  ( 1 .  821 

256*NTi:4)  +  (1.330274*NTC5)  ) 

7060  IF  XN<0  GOTO  7090 

7070  F0X=1-Q0X 

7080  RETURN 

7090  FOX=QOX 

7100  D0X=1-F0X 

7110  RETURN 

3000  REM  %    t    %       OPTION  SELECTION  ROUTINE  *  « 

8010  PRINT  S>JP+5,  "NUMBER  OF  OPTION  DESIRED 

" ; LBS ; CHR$ ( 1 76 ) ; RB$ 

8020  ZI=INKEY$:IF  ZI=""  THEN  S020 

8030  PRINT  3JP+31,ZI; : IO=VAL(ZI) 

8040  IF  I0<1  OR  IO>KS  THEN  PRINT  3JP,":*:t*   NOT  A  VALID 

OPTION   **:*■';:  FOR  1  =  1  TO  1000:  NEXT:  PRINT  3JP,BL$:G0T0  8010 

10000  REM  %     t     %    ARRAY  SORT  t     % 

10010  L=l  :  B(L)=N1+1 

10020  J=B(L)  ^SET  END  OF  ARRAY  SEGMENT 

10030  I=M-1   'SET  START  OF  ARRAY  SEGMENT 

10040  IF  J-M  <  3  THEN  10230    'HANDLE  1  OR  2  ELEMENTS  IN 

SPECIAL  CASE 

10050    M1  =  INT(  (I+J) /2)       -'SET  COMPARE  ELEMENTS 

10060    REM  FIND  A  LARGE  ELEMENT  AMONG  THE  SMALL  ONES 

10070    1=1+1 

10080    IF  I=J  THEN  10170 


yy. 


10090 

10100 

10110 

10120 

10130 

10140 

10150 

10160 

ELEMENT 

10170 

lOlSO 

10190 

1 0200 

LARGE 

1 02 1 0 

10220 

10230 

1 0240 

1 0250 

10260 

LARGE 

10270 

10280 

10290 

11000 

11010 

11020 

11030 

11040 

11050 

11060 

1 1 070 

11080 

1 1 090 

12000 


IF  X(P(I))  <=  X(P(M1))  THEN  10070 

REM  FIND  A  SMALL  ELEMENT  AMONG  THE  LARGE  ONES 

J=J-1 

IF  I=J  THEN  10170 

IF  X(P(J))  >=  X(P(M1))  THEN  10110 
REM  EXCHANGE  OUT  OF  PLACE  ELEMENTS 
H=P(I)  :  P(I)=P(J)  :  P(J)=H  :  GOTO  10070 
REM  ADJUST  COMPARE  ELEMENT  TO  NEW  MIDDLE 


IF  I  >  Ml  THEN  1=1-1 
IF  J=M1  THEN  10210 

H=P(I)  :  P(I)=P(M1)  :  P(M1)=H 

REM  SAVE  STARTING  POINT  FOR  ARRAY  SEGMENT  OF 
ELEMENTS 

L=L+1  :  B(L)=I 
GOTO  10020 
IF  J-M  <  2  THEN  10270 

IF  X(P(M))  <=  X(P(M+1))  THEN  10260 
H=P(M)  :  P'M)=P(M+1>  :  P(M+1)=H 
REM  SET  BEGIN  AND  ENDPOINTS  FDR  ARRAY  SEGMENT  OF 
ELEMENTS 
M=B(L)+1  :  L=L-1 
IF  L>0  THEN  GOTO  10020 
RETURN 

REM  «  *  *  READ  FILE  DIRECTORY  *  * 
CLS:PRINTa20, "MANN-WHITNEY  TEST": PRINT 
F$="EDITDATA/DIR: "+DN$ 
OPEN  "  I  " , 1 , F* 
INPUT  #1,F0 
FOR  1=1  TO  FO 

INPUT  #l,F*':i) 
NEXT  I 
CLOSE  1 
RETURN 

REM  «  *  *   PRINT  FILE  DIRECTORY  *  t 
1 20 1 O  CLS : PR I NT320 , " MANN-WH I TNE Y  TEST " : PR I NT 
12020  FOR  1=1  TO  FO  STEP  4 
FOR  J=l  TO  4 

IF  I+J-1>F0  THEN  12090 

IF  J=l  AND  LENvF$(I+J-l)  X12  AND  I>1  THEN 


12030 

12040 

12050 

PRINT  "" 

12060 

";F$(I+J-1) 

12070 

12080 

12090 

14000 

14010 

14020 


PRINT  TAB(16*(J-1)) 


"  ("jSTR^d+J-D 


14030 
14040 


NEXT  J 
NEXT  I 
RETURN 

REM  t    t     t    AVERAGE  RANPxS  *  * 
IF  JS=JE  THEN  X(P(JS))=JS  :  RETURN 
AV=0:FOR  I=JS  TO  JE:  AV=AV+I:NEXT  I 
FOR  I=JS  TO  JE: X (P(I) )=AV:  NEXT  I 
RETURN 


AV=AV/ (JE-JS+1) 
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19000  =•  ERROR  HANDLING 

19010  IF  ERR=106  AND  (SL=2005  OR  SL=2410  OR  SL=2710  OR  3110 

OR  SL=3510  OR  SL=3910)  THEN  INPUT  "NO  FILES  ON  THIS  DRIVE. 

PRESS  ENTER  TO  CONTINUE" ;A$  :  SL=0  :  RESUME  11080 

19100  PRINT  "ERR=";ERR; "ERL=";ERL  :  ON  ERROR  GOTO  O 


R-F^OC3IRi=»r^   4£3.  ,    Sr^  X  F<SvJC3V/:Brf=^^ 


10  ON  ERROR  GOTO  19000 

20  DEFSTR  Z:  DEFINT  P 

30  DIM  F*(95) ,P(1500) ,X1 (1500) ,B(50) , X2(1500) , X (1500) 

40  LB$=CHR* ( 1 23 ) : RB$=CHR$ ( 1 25 ) : MB$=CHR$ (176): BL*=CHR$ ( 252 ) 

50  REM  *    t    t       SMIRNOV  TEST   *  » 

60  CLS:  PRINT  5)20,  "SMIRNOV  TEST" 

70  PRINT:PRINT:PRINT  "Disk  Drive  Containing  Data  Files 

(0-3)7" 

80  DN$=INKEY*: IF  DN$=""  THEN  SO  ELSE  I=VAL(DN*) 

90  IF  KO  OR  I>3  THEN  GOTO  70 

100  SL=li0  :  GOSUB  11000  :  IF  SL=0  THEN  60  ELSE  GOSUB  13000 

:  PRINT  "" 

110  PRINT: INPUT  "Enter  the  Two  Data  Files  (e.g.  3,5)";F1,F2 

120  IF  F1>F0  OR  F2::F0  OR  FKl  OR  F2<  1  THEN  110 

130  OPEN  "I", 1,F*(F1)+": "+DN* 

140  INPUT  #1,N 

150  FOR  1=1  TO  N  :  INPUT  #1,X(I)  :  NEXT  I 

160  CLOSE  1 

170  N1=N:M=1:F0R  1=1  TO  N: P ( I ) =1 : NEXT  I : GOSUB  10000 

ISO  FOR  1=1  TO  N: XI (I)=X (P(I) ) :NEXT  I 

190  OPEN  "I", 1,F*(F2)+": "+DN* 

200  INPUT  #1,MM 

210  FOR  1=1  TO  MM  :  INPUT  #1,X(I)  :  NEXT  I 

220  CLOSE  1 

230  N1=MM:M=1:FGR  1=1  TO  MM: P ( I ) =1 : NEXT  I : GOSUB  lOOOO 

240  FOR  1=1  TO  MM:X2(I)=X (P(I) ) zNEXT  I 

242  IF  MM=>N  THEN  GOTO  250 

244  FOR  1=1  TO  N:X(I)=X1(I)  :  NEXT  I 

245  FOR  1=1  TO  MM  :  X1(I)=X2(I)  :  NEXT  I 

246  FOR  1=1  TO  N  :  X2 ( I ) =X ( I )  :  NEXT  I 
248  I=N  :  N=MM  :  MM=I 

250  1=1: J=l: XCDF=0:YCDF=0 

260  TP=0:TM=0 

270  IF  I>N  AND  J >MM  THEN  GOTO  380 

280  IF  X1(I)=X2(J)  THEN  I  =  Ih-1  :  J=J  +  1 :  GOTO  330 

290  IF  XI (I)  <  X2(J)  THEN  GOTO  300  ELSE  GOTO  310 

300  IF  I  <=N  THEN  1=1+1: GOTO  330: ELSE  IF  J<=MM  THEN 

J=J+1:G0T0  330: ELSE  SOTO  380 

310  IF  XI (I)  >  X2(J)  THEN  GOTO  320  ELSE  GOTO  330 

320  IF  J  <=MM  THEN  J=J+1  ELSE  IF  I  <=N  THEN  1=1+1  ELSE  GOTO 

380 

330  YCDF= ( I - 1 ) /N  :  XCDF= ( J- 1 ) /MM 

340  D=YCDF-XCDF 

350  IF  D<TM  THEN  TM=D 

360  IF  D>TP  THEN  TP=D 

370  GOTO  270 

3S0  CLS: PRINT320, "SMIRNOV  TEST" : PRINT: PRINT 


/o 


390  PRINT  "HYPOTHESIS  NULL 

ALTERNATIVE" 

400  PRINT  "     (1>  F(X)   =   G(X)  F(X)  =/= 

G  <  X  )  " 

410  PRINT  "     (2)  F(X)  <=   G(X)  F(X)    > 

G  (  X  )  " 

420  PRINT  "     (3)  F(X)  =>   G(X)  F(X)   < 

G(X)  " 

430  PRINT: PRINT  "Where  F(X)  and  G(X)  represent  Data  Sets  1 

it    2  Respectively" 

440  JP=640:KS=3:G0SUB  SOOO 

450  CLS:PRINTa)20,  "SMIRNOV  TEST"  :  PRINT:  PRINT 

460  ON  10  GOTO  470,480,490 

470  IF  TP>ABS(TM)  THEN  D=TP  ELSE  D=ABS (TM) : GOTO  500 

480  D=TP:GOTO  500 

490  D=ABS(TM) 

500  PRINT  "Test  Statistic  =  ";D 

505  PRINT  "Sample  Sizes  are  ";N;: PRINT  "  and  " ; MM 

510  GOSUB  2000 

520  ON  10  GOSUB  600,800,1000 

530  PRINT  5)960,  "Press  SPACE  BAR  when  ready  tc  return  to 

main  menu"; 

540  IF  INKEY$=  ""  GOTO  540  ELSE  RUN  "MENU/3AS" 

600  REM  t     t    t    NULL  HYPOTHESIS  #1  *  * 

610  PRINT  "P-Value  =  ";:PRINT  USING  "#.####"; 2  t    QOX 

620  PRINT: INPUT"Enter  Desired  Alpha  Level ";AL 

630  IF  AL<0  OR  AL>1  THEN  PRINT  "*:**  ERROR  ***  ENTER  A  VLAUE 

BETWEEN  0  &  1 " : GOTO  620 

640  PRINT: PRINT  "For  Alpha  =";AL;",  "; 

650  IF  AL<2«Q0X  THEN  PRINT  "accept";  ELSE  PRINT  "reject": 

660  PRINT  "  the  hypothesis  that  F-X)  =  G(X)" 

670  PRINT  "-VS-  the  alternative  that  F(X)  =/=  G(X)." 

680  RETURN 

800  REM  *  *  *  NULL  HYPOTHESIS  #2  *  * 

810  PRINT  "P-Value  =  ";:PRINT  USING  "#.####";QOX 

820  PRINT: INPUT"Enter  Desired  Alpha  Level";AL 

330  IF  hL<0  OR  AL>1  THEN  PRINT  "«**  ERROR  *«:«  ENTER  A  VLAUE 

BETWEEN  O  S<  1  "  :  GOTO  620 

840  PRINT:PRINT  "For  Alpha  =";AL;",  "; 

850  IF  AL<QOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 

860  PRINT  "  the  hypothesis  that  F(X)  <=  G(X)" 

870  PRINT  "-VS-  the  alternative  that  F(X)  >  G<X)." 

880  RETURN 

1000  REM  t    t     t    NULL  HYPOTHESIS  #3  «  * 

1010  PRINT  "P-Value  =  ";: PRINT  USING  "#.####"; QOX 

1020  PRINT: INPUT "Enter  Desired  Alpha  Level "; AL 

1030  IF  AL<0  OR  AL>1  THEN  PRINT  "**«  ERROR  ttt    ENTER  A 

VLAUE  BETWEEN  O  &  1 " : GOTO  620 

1040  PRINT:PRINT  "For  Alpha  =";AL:",  "; 

1050  IF  AL<QOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 

1060  PRINT  "  the  hypothesis  that  F(X>  =>  G ( X ) " 

1070  PRINT  "-VS-  the  alternative  that  F(X)  <  G(X)." 
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*  *   SMIRNOV   DISTRIBUTION   *  * 


12+SQR(N)+. 11/SQR(N) ) 
*  3QR(N/ (N+M) ) 
S=D  t    SQR(N/ (N+M) ) 


1080  RETURN 
2000  REM  * 
2010  M=MM 

2020  IF  N<=.1  THEN  S= (D-1 / (2) ) * ( 
ELSE  IF  (M/N)=INT(M/N)  THEN  S=D 

+2/(3  t    SQR(M) )  ELSE 

+2/(5  «  SQR(M)) 
2030  QOX=(. 135) C (SC2) 
2040  IF  QOX<0  THEN  QOX=0 
2050  IF  Q0X>1  THEN  Q0X=1 
2060  F0X=1-Q0X 
2070  PX=D 
2080  RETURN 

8000  REM  t    t    t       OPTION  SELECTION  ROUTINE  *  * 
3010  PRINT  5) JP+5,  "NUMBER  OF  OPTION  DESIRED 
" ; LB$; CHR$ ( 176) ; RB$ 

8020  2I=INKEY$:IF  21=""  THEN  8020 
8030  PRINT  o)JP+31,ZI;  :  IO=VAL(ZI) 


THEN  PRINT  3JP,"*«* 
TO  1000: NEXT: PRINT 
SORT   *  * 


SJP,BL$i 


VAL I D 
GOTO  SO  10 


8040  IF  I0<1  OR  IO>KS  THEN  PRINT  9JP, "ttt       NOT  A 
OPTION   ***";: FOR  1=1 
10000  REM  ♦  *  t    ARRAY 
L=l  :  B(L)=N1+1 

J=B(L)  'SET  END  OF  ARRAY  SEGMENT 
I=M-1   'SET  START  OF  ARRAY  SEGMENT 

IF  J-M  <  3  THEN  10230    'HANDLE  1  OR  2  ELEMENTS  IN 
CASE 

M1=INT( (I+J) /2)       'SET  COMPARE  ELEMENTS 
REM  FIND  A  LARGE  ELEMENT  AMONG  THE  SMALL  ONES 
1  =  1  +  1 
IF  I=J  THEN  10170 


10010 

10020 

1 0030 

1  0040 

SPECIAL 

1 0050 

10060 

10070 

10080 

10090 

10100 

10110 

10120 

10130 

10140 

10150 

10160 

ELEMENT 

101 70 

10130 

10190 

1 0200 

LARGE 

10210 

10220 

10230 

10240 

10250 

1 0260 

LARGE 

10270 


IF  X(P(I))  <=  X(P(M1))  THEN  10070 

REM  FIND  A  SMALL  ELEMENT  AMONG  THE  LARGE  ONES 

J=J-1 

IF  I=J  THEN  10170 

IF  X(P(J))  >=  X(P(M1))  THEN  10110 
REM  EXCHANGE  OUT  OF  PLACE  ELEMENTS 
H=P(I)  :  P(I)=P(J)  :  P(J)=H  :  GOTO  10070 
REM  ADJUST  COMPARE  ELEMENT  TO  NEW  MIDDLE 


IF  I  >  Ml  THEN  1=1-1 
IF  J=M1  THEN  10210 

H=P(I)  :  P(I)=P(M1)  :  P(M1)=H 

REM  SAVE  STARTING  POINT  FOR  ARRAY  SEGMENT  OF 
ELEMENTS 

L=L+1  :  B(L)=I 
GOTO  10020 
IF  J-M  <  2  THEN  10270 

IF  X(P(M))  <=  X(P(M+1))  THEN  10260 
H=P(M)  :  P(M)=P(M+1)  :  P(M+1)=H 
REM  SET  BEGIN  AND  ENDPOINTS  FOR  ARRAY  SEGMENT 
ELEMENTS 
M=B<L)+1  :  L=L-1 
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102S0  IF  L>0  THEN  GOTO  10020 

10290  RETURN 

11000  REM   *  *  X   READ  DATA  FILE  DIRECTORY   *  * 

11010  CLS:PRINT320, "SMIRNOV  TEST": PR I NT 

11020  F*="EDITDATA/DIR: "+DN$ 

11030  OPEN  "I", 1,F$ 

11040  INPUT  #1,F0 

11050  FOR  1=1  TO  FO 

11060       INPUT  #1,F$(I) 

11070  NEXT  I 

11080  CLOSE  1 

11090  RETURN 

13000  REM   *  *  *   PRINT  DATA  FILE  DIRECTORY   :*  * 

1 30 1 O  CLS : PR I NT320 , " SM I RNO V  TEST " : PR I NT 

13020  FOR  1=1  TO  FO  STEP  4 

13030      FOR  J=l  TO  4 

13040  IF  I+J-1>F0  THEN  13090 

13050  IF  J=l  AND  LEN (F$ ( I+J-1 ) ) < 12  AND  I>1  THEN 

PR I NT  " " 

13060  PRINT  TAB(16*(J-1)  );"  <"|;STR$(  I+J-1)  ; '■) 

":F*( I+J-1) ; 

13070      NEXT  J 

13080  NEXT  I 

13090  RETURN 

19000  'ERROR  HANDLING 

19010  IF  ERR=106  AND  <SL=2005  OR  SL=2410  OR  SL=2710  OR  3110 

OR  SL=3510  OR  SL=3910)  THEN  INPUT  "NO  FILES  ON  THIS  DRIVE. 

PRESS  ENTER  TO  CONTINUE" ;A$  :  SL=0  :  RESUME  11080 

19100  PRINT  '■ERR=";ERR;  "ERL=";ERL  :  ON  ERROR  GOTO  O 


7Q 


F='F^ZOGF^^=\t^ 


|vJL_  I  l_l_/Ei^S 


N    G 


DATA 


GOSUB    1 OOOO 


5  CLEAR  1000 
10  ON  ERROR  GOTO  19000 
20  DEFBTR  Z:  DEFINT  P 

30  LB$=CHR$ (123) : RB$=CHR* ( 125) : MB*=CHR* ( 176) : BL*=CHR* (252) 
40  DIM  F$(95) ,P(1500) , XI (1500) ,B(50) ,X(1500) 

50  LB$=CHR$ ( 1 23 ) : RB*=CHR$ ( 1 25 ) : MB$=CHR$ (176): BL$=CHR$ ( 252 ) 
60  REM   «  «  «   LILLIEFOR  TEST  FOR  NORMALITY   *  * 
70  CLS:PRINT  .520,  "N  ORMALITY      TEST" 
SO  PRINT: PRINT: PRINT  "Disk  Drive  Containing  Data  Files 
(0-3)7  "; 

90  DN$=INKEY*: IF  DN$=""  THEN  90  ELSE  PRINT  DN$:  I=VAL(DN$) 
100  IF  KO  OR  I>3  THEN  GOTO  80 

110  SL=110  :  GOSUB  11000  :  IF  SL=0  THEN  70  ELSE  GOSUB  12000 
:  PRINT  "" 

120  PRINT: INPUT  "Enter  the  Data  File  #";F1 
130  IF  F1>F0  OR  FKl   THEN  120 
140  OPEN  "I", 1,F$(F1)+": "+DN$ 
150  INPUT  #1,N 

160  FOR  1=1  TO  N  :  INPUT  #1,X(I)  :  NEXT  I 
170  CLOSE  1 

ISO    CLS:  PRINT    5)458,     "t     t     t       ORDER 
*    «" 

190  N1=N:M=1:F0R  1=1  TO  N: P ( I ) =1 : NEXT 
200  SM=0 : SUM=0 : S2=0 
210  FOR  1=1  TO 

N:X1 (I)=X (P(I) ) :SM=SM+X1 (I) :DV=X1 ( I ) -SUM: SUM=SUM+DV/ I : S2=S2+ 
DV* ( X 1 ( I ) -SUM) : NEXT  I 
220  CLS: PRINT  "SELECT  OPTION  " 

230  PRINT:PRINT:PRINT  "(1)    MEAN  and  VARIANCE  Ar&    both 
unknown " 

240  PRINT: PRINT  "(2)    MEAN 
250  PR I NT : PR I NT  " ( 3 ) 
260  JP=S96:KS=3: GOSUB 
290  ON  10  GOTO  300,320,340 

300  MU=SM/N: VR=S2/ (N-1 ) : SD=SQR (VR) : XCDF=0: T=0 
310  GOTO  370 

320  CLS: INPUT  "Enter  the  MEAN";MU 
330  VR=0:FOR  1  =  1  TO  N:  VR=VR-h  (  X  1  ( I ) -MU)  112:  NEXT 
I:SD=SQR(VR/ (N-1) ) :GOTO  370 
340  CLS:  INPUT  "Enter  the  MEAN";MU 
350  INPUT  "Enter  the  VARIANCE" ; VR: SD=SDR (VR) 

360  ^F  VR<0  THEN  PRINT  "*  «  *   ERROR  *  *  *   VARIANCE  MUST 
BE  POSITIVE": GOTO  350 
370  FOR  1=  1  TO  N 

380  XN=(X1 (I)-MU) /SD: AX=ABS(XN) 
390  XZ=.3989423XP( (-AX»AX) /2) 
400  NT=1/ (1+.2316419*AX) 


is  known,  VARIANCE  is  unknown" 
MEAN  and  VARIANCE  a^rs    both  known" 
8000 


30 


7S1478«NTC3) - ( 1 . S21 


L  I  T  Y 
Statistic  = 


T  E 


S  T 


410 

Q0X=X2* ( ( . 3193S15*NT) - ( . 3565638*NTC2) + ( 1 

256«NTi:4)  +  (1.33274:*NTi:5)  ) 

420  IF  XN<0  THEN  Q0X=1-Q0X 

F0X=1-Q0X 

D=ABS(XCDF-FOX) 

IF  D>T  THEN  T=D 

XCDF=I/N 

D=ABS(XCDF-FQX) 

IF  D>T  THEN  T=D 

NEXT  I 

CLS:PRINT320, "N  O  R  M  A 

PR I NT : PR I NT :PK I NT  "Test 

PRINT  "Sai..ple  Size  =  " ;  N 

ON  ID  GOSUB  2000,3000,4000 

PKINT  "P-Value  =  ";: PRINT  USING 

PRINT: INPUT"Enter  Desired  Alpha 

IF  AL<0  OR  AL>1  THEN  PRINT  ":**:* 
BETWEEN  O  Zc    1  "  :  GOTO  590 
610  PRINT: PRINT  "For  Alpha  =";AL;",  "; 

IF  AL<QOX  THEN  PRINT  "accept";  ELSE  PRINT 

PRINT  "  the  hypothesis  that  the  sample  is 

PRINT  "a  normal  population." 

PRINT  3960, "Press  SPACE  BAR  when  ready  to 
main  menu"; 

310  IF  INKEY$=  ""  GOTO  310  ELSE  RUN  "MENU/BAS" 

2000  REM  :*  t  *  LILLIEFQR^S  DISTRIBUTION  -  MEAN  AND  VARIANCE 
UNKNOWN 
2020  DENCM=1 
2030  D=  T/DENOM 


430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
580 
590 
600 


620 
630 
640 
800 


"#.####" ;QOX 
Level ":AL 
ERROR  ;**:* 


ENTER  A  VLAUE 

"reject"; 
from" 

return  to 


2040 

S=D*  (--01-^SQR(N)+.S5/SQR(N)  ) 

2050 

TDENOM=DENOh 

2060 

GOSUB  2500 

2070 

IF  FOX  <=  .85  THEN  DENGM  =  1-.24  *  FOX: GOTO  2200 

2080 

IF  FOX  <=  .9  THEN  DENOM  =  1.357  -  .66*  FOX  :  GOTO  22' 

2090 

IF  FOX  <:>  .95  THEN  DENOM  =  1.339  -  .64*  FOX  :  GOTO 

2200 

2100 

IF  FOX  <=  .975  THEN  DENGM  =  1.795  -  1.12*  FOX  : 

GOTO 

2200 

2110 

DENOM  =  2.133  -  1.467*  FOX 

2200 

IF  ABS(TDENOM  -  DENOM)  >  .001  THEN  GOTO  2030 

22 1 0 

RETURN 

2500 

QOX  =  . 135C '8:2) 

25 1 0 

IF  QOX  >  1  THEN  QOX  =  1 

2520 

FOX  =  1  -  QOX 

2530 

RETURN 

2540 

IF  PR  =>  .70  THEN  XQ  =  0.743*PR  +  0.2799  :  GOTO 

2590 

2550 

IF  PR  =>  .50  THEN  XQ  =  0.546*PR  +  0.4187  :  GOTO 

2590 

2560 

IF  PR  =>  .30  THEN  XQ  =  0.464*PR  +  0.4588  :  GOTO 

2590 

2570 

IF  PR  =>  .20  THEN  XQ  =  0. 4S5*PR  +  0.4525  :  GOTO 

2590 

2580 

IF  PR  =>  .10  THEN  XQ  =  0.587*PR  +  0.4321  :  GOTO 

2590 

2585 

XQ  =  .846* PR  +  0.4062 

81 


2590  XQ=XQ/DENOM 

2600  XQ=XQ-»-.20/N 

2610  RETURN 

3000  REM  *  *  *  LILLIEFOR'S  DISTRIBUTION  -  MEAN  KNOWN  AND 

VARIANCE  UNKNOWN  :  %    t 

3010  DENOM  =  1  :  PX  =  D 

3020  DENOM  =  PX /DENOM 

3022  S=D* ( . 055+SQR (N) +. 2S3/SQR (N) ) 

3026  TDENOM  =  DENOM 

3030  GOSUB  2500 

3040  IF  FOX  <=  .9  THEN  DENOM  =  1.83  -.8*  FOX  :  GOTO  3200 

3050  IF  FOX  <=  .95  THEN  DENOM  =  1.47  -  .4*  FOX  :  GOTO  3200 

3060  IF  FOX  <=  .975  THEN  DENOM  =  1.85  -  . S*  FOX  :  GOTO  3200 

3070  DENOM  =  1.07 

3200  IF  ABS<TDENaM-DENOM)  >  .001  THEN  GOTO  3020 

3210  RETURN 

4000  REM  *  *  *   KOLMOGOROV  DISTRIBUTION  *  * 

4010  S=T:*  (O.  12+SQR(N)+.  11/SQR(N)  ) 

4020  QOX  =  .  135C  (3!:2) 

4030  F0X=1-Q0X 

4040  RETURN 

8000  REM  *  *  lc   OPTION  SELECTION  ROUTINE  %    % 

8010  PRINT  5) JP+5,  "NUMBER  OF  OPTION  DESIRED 

" ; LB$ ; CHR$ (176); RE* 

8020  ZI=INKEY$:IF  21=""  THEN  8020 

8030  PR  I  NT  a)JP+3 1 ,  Z  I ;  :  1 0=VAL  (  Z  I ) 

8040  IF  I0<1  OR  IO>KS  THEN  PRINT  3JP,"«**   NOT  A  VALID 

OPTION   ***";:FOR  1=1  TO  1000: NEXT: PRINT  3JP,BL$:G0T0  30 lO 

8050  RETURN 

lOOOO  REM   *  *  *   ARRAY  SORT  %    % 

10010  L=l  :  B(L)=Nl-!-l 

10020  J=B(L)  'SET  END  OF  ARRAY  SEGMENT 

10030  I=M-1   'SET  START  OF  ARRAY  SEGMENT 

10040  IF  J-M  <  3  THEN  10230    'HANDLE  1  OR  2  ELEMENTS  IN 

SPECIAL  CASE 

1 0050    M 1  =  I  NT  ■:  ( I  +  J  )  /2 )       '  SET  COMPARE  ELEMENTS 

10060    REM  FIND  A  LARGE  ELEMENT  AMONG  THE  SMALL  ONES 

10070    1=1+1 

10080    IF  I=J  THEN  10170 

10090        IF  X<P<I))  <=  X(P(M1))  THEN  10070 

10100        REM  FIND  A  SMALL  ELEMENT  AMONG  THE  LARGE  ONES 

10110       J=J-1 

10120        IF  I=J  THEN  10170 

10130  IF  X(P(J))  >=  X':P(M1))  THEN  10110 

1014O  REM  EXCHANGE  OUT  OF  PLACE  ELEMENTS 

10150  H=P(I)  :  P(I)=P(J)  :  P<J)=H  :  GOTO  10070 

10160  REM  ADJUST  COMPARE  ELEMENT  TO  NEW  MIDDLE 

ELEMENT 

10170     IF  I  >  Ml  THEN  1=1-1 

10130     IF  J=M1  THEN  10210 

10190         H=P(I)  :  P(I)=P':M1)  :  P(M1)=H 

10200         REM  SAVE  STARTING  POINT  FOR  ARRAY  SEGMENT  OF 


32 


LARGE 
10210 
10220 
10230 
10240 
10250 
10260 
LARGE 
10270 
10280 
10290 
1 1 000 
11010 
11020 
1 1 030 
1 1 040 
11050 
1 1 060 
11070 
1 1 080 
11090 
12000 
12010 
12020 
12030 
12040 
12050 
PRINT 
1 2060 
";F$( 
12070 
1 2080 
12090 
1 9000 
19010 
OR  SL 
PRESS 
19100 


ELEMEhJTS 

L=L-»-l  :  B'L)=I 

GGTD  10020 
IF  J-M  <  2  THEN  10270 

IF  X(P(M))  <=  X(P(M-e-l))  THEN  10260 
H=P(M)  :  P(M)=P(M+1)  :  P(M+1)=H 
REM  SET  BEGIN  AND  ENDPOINTS  FOR  ARRAY  SEGMENT 
ELEMENTS 
M=B(L)+1  :  L=L-1 
IF  L>0  THEN  GOTO  10020 
RETURN 
REM  t     t    *    READ  DATA  FILE  DIRECTORY   *  * 


OF 


CLS:PRINTS)20,  "N  O  R  M  A 
F$="EDITDATA/DIR:  "-!-DN$ 
OPEN  " I " , 1 , F* 
INPUT  #1,F0 
FDR  1=1  TO  FO 

INPUT  #1,F$(I) 
NEXT  I 
CLOSE  1 
RETURN 


L  I  T  Y 


I 


E  S  T": PRINT 


REM 

CLS: 

FOR 


t     t    t       LIST  DATA 
PRINT320, "NORM 


FILE  DIREC 
A  L  I  T  Y 


ORY 


* 


T": PRINT 


1=1  TO  FO  STEP  4 

FOR  J=l  TO  4 

IF  I-!-J-l>FO  THEN  12090 

IF  J=l  AND  LEN(F$(I+J-1) )<12 


AND  I>1  THEN 


PRINT  TAB(16:*(J-1)  )  ;  "  (  "  :  STR$  ( I+J-1 )  ;  ") 
I+J-1); 

NEXT  J 

NEXT  I 

RETURN 

' ERROR  HANDLING 

IF  ERR=106  AND  (SL=2005  OR  SL=2410  OR  SL=2710 
=3510  OR  SL=3910)  THEN  INPUT  "NO  FILES  ON  THIS 

ENTER  TO  C0NTINUE";A5  :  SL=0  :  RESUME  1 1080 

PRINT  "ERR=";ERR; "ERL=";ERL  :  ON  ERROR  GOTO  O 


OR  3110 
DRIVE. 


F=-IRO<3iRrf=%l^ 


e:l_I!_!_/b*=»s 


5  CLEAR  1000 
10  ON  ERROR  GOTO  19000 
20  DEFSTR  Z:  DEFINT  P 

30  LB$=CHR$ (123): RB*=CHR$ ( 1 25 )  : MB$=CHR$  < 1 76 ) : 
40  DIM  F$(95) ,P(1500) ,X1 <1500) ,B<50) ,X (1500) 
50  LB$=CHR* ( 1 23 ) : RB*=CHR* ( 1 25 ) : MB$=CHR$ ( 1 76 ) : 
60  REM   *  :«  *   EXPONEhiTIAL  DISTRIBUTION  TEST 
CLS: PRINT  ^20, "E  XPONENTIAL     T 


"Disk  Drive  Containing  Data  Files 


DN$:  I=VAL(DN$) 


70 

SO  PRINT: PRINT: PRINT 
(0-3)?  "; 

90  DN$=INKEY$: IF  DN*=""  THEN  90  ELSE  PRINT 
100  IF  KO  OR  I>3  THEN  GOTO  80 
110  SL=110  :  GOSUB  11000  :  IF  SL=0  THEN  70 
:  PRINT  "" 
120  PRINT: INPUT  "Enter  the  Data  File  #";F1 

IF  F1>F0  OR  FKl   THEN  120 

OPEN  "I", 1,F$(F1) +": "+DN$ 

INPUT  #1,N 

FOR  1=1  TO  N  :  INPUT  #1,X(I)  :  NEXT  I 

CLOSE  1 

CLS: PRINT  3453,  "*  *  *   O  R  D  E  R  I 


BL$=CHR$(252) 

BL$=CHR$(252) 

t     t 
EST" 


ELSE  GOSUB  12000 


130 

140 

150 

160 

170 

ISO 

*  *" 

190 

200 

210 


N  G 


DATA 


N1=N:M=1:F0R  1=1  TO  N: P ( I ) =1 : NEXT  I : GOSUB  10000 

SM=0 

FOR  1=1  TO  N: XI (I) =X (P(I) ) :SM=SM+X1 (I) :NEXT  I:LA=N/SM 


220  CLS: PRINT  "SELECT  OPTION  " 

230  PRINT: PRINT: PRINT  "(1)     LAMBDA  is  unknown 

240  PR I NT: PR I NT  "(2)     LAMBDA  is  known" 

250  •JP=S96:KS=2:  GOSUB  3000 


260 
270 
'7  on 


ON  10  GOTO  290,270 


CLS: INPUT  "Enter  LAMBDA": LA 
IF  LA  <=  O  THEN  PRINT  "*  *  * 

BE  POSITIVE": GOTO  270 

290  T=0:XCDF=0 

300  FOR  1=1  TO  N 

310  FOX  =  1 -EXP (-LAI (I) ) 

320  D=ABS ( XCDF-FOX ) 

330  IF  D  >  T  THEN  T=D 

340  XCDF  =  I/N 

350  D=ABS( XCDF-FOX) 

360  IF  D  >  T  THEN  T=D 

370  NEXT  I 

380  CLS:PRINT5)20,  "E  X 


ERROR   *  *  *   LAMBDA  MUST 


590  PRINT: PRINT: PRINT 


400  PRINT  "Sample  Size  = 
410  ON  10  GOSUB  2000,3000 


PONENTIAL 
"Test    Statistic    = 
"  ;  N 


TEST 


";T 


420    PRINT    "P-Value    = 


PRINT    USING    "#.####" :QOX 


84 


430  PRINT: INPUT"Enter  Desired  Alpha  Level ";AL 

440  IF  AL<0  OR  AL>1  THEN  PRINT  "ttt    ERROR  *:**  ENTER  A  VLAUE 

BETWEEN  O  &  1 " : GOTO  430 

450  PRINT: PRINT  "For  Alpha  =";AL:",  "; 

460  IF  AL<QOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 

470  PRINT  "  the  hypothesis  that  the  sample  is  from" 

480  PRINT  "a  normal  population." 

490  PRINT  3960, "Press  SPACE  BAR  when  ready  to  return  to 

main  menu"; 

500  IF  INKEY$=  ""  GOTO  500  ELSE  RUN  "MENU/BAS" 

2000  REM  *  *  «  LILLIEFOR'S  DISTRIBUTION  -  MEAN  AND  VARIANCE 

UNKNOWN 

2020  DENOM=SQR (N) +. 26+. 5/SQR (N) 

2040  IF  T  >  1.3442/DENOM  +  . 2/N  THEN  FOX  =  1  :  GOTO  2220 

2060  F0X=.5:G0SUB  22S0 

2080  IF  T>XQ  THEN  F1=F0X : F2=l : GOTO  2120 

2100  F1=0  :  F2=F0X 

2120  FOX  =  (Fl+F2)/2 

2140  GOSUB  2280 

2160  IF  ABS'T-XQX.OOOl  THEN  GOTO  2220 

2180  IF  T::XQ  THEN  F1=F0X  ELSE  F2=F0X 

2200  GOTO  2120 

2220  Q0X=1-F0X 

2230  RETURN 

2280  IF  FOX  =>  .95  THEN  XQ  =  5.125*  FOX  -  3. 7808  :  GOTO 

2460 

2300  IF  FOX  =>  .90  THEN  XQ  =  2-026*  FOX  -  0-8367  :  GOTO 

2460 

2320  IF  FOX  =>  .80  THEN  XQ  =  1.124*  FOX  -  0.0249  :  GOTO 

2460 

2340  IF  FOX  =>  .70  THEN  XQ  =  0.743*  FOX  +  0.2799  :  GOTO 

2460 

2360  IF  FOX  =>  .50  THEN  XQ  =  0.546*  FOX  +  0.4137  :  GOTO 

2460 

2380  IF  FOX  =>  .30  THEN  XQ  =  0.464*  FOX  +  0.4588  :  GOTO 

2460 

2400  IF  FOX  =>  .20  THEN  XQ  =  0.485*  FOX  +  0.4525  :  GOTO 

2460 

2420  IF  FOX  =>  .10  THEN  XQ  =  0.587*  FOX  +  0.4321  :  GOTO 

2460 

2440  XQ  =  .846*  FOX  +  0.4062 

2460  XQ=XQ/DENOM  +  .2/N 

2470  RETURN 

3000  REM  *  *  *  LILLIEFOR'S  DISTRIBUTION  -  LAMEA  KNOWN  :  * 

3010  S=T* (O. 12+SQR(N)+. 11/SQR<N) ) 

3020  QOX=. 135n (SC2) 

3030  IF  Q0X>1  THEN  Q0X=1 

3040  F0X=1-Q0X 

3050  PX=T 

3060  RETURN 

3000  REM  *  «  *   OPTION  SELECTION  ROUTINE  *  * 


8! 


VALID 
GOTO  SO 10 


IN 


SOlO  PRINT  3 JP+5, "NUMBER  OF  OPTION  DESIRED 
" ; LB* ; CHR$ (176); RB* 

S020  ZI=INKEY$:IF  ZI=""'  THEN  8020 
8030  PRINT  3JP+31,2I; : IO=VAL(ZI> 

8040  IF  I0<1  OR  IO>KS  THEN  PRINT  5)JP,"«*X  NOT  A 
OPTION  ***";:FOR  1=1  TO  1000: NEXT: PRINT  3JP,BL$ 
8050  RETURN 

lOOOO  REM  t    t    t       ARRAY  SORT   «  * 
L=l  :  B(L)=N1+1 

J=B(L)  'SET  END  OF  ARRAY  SEGMENT 
I=M-1   'SET  START  OF  ARRAY  SEGMENT 
IF  J-M  <  3  THEN  10230    ="  HANDLE  1  OR  2  ELEMENTS 
CASE 

M1  =  INT(  (H-J)/2)       'SET  COMPARE  ELEMENTS 
REM  FIND  A  LARGE  ELEMENT  AMONG  THE  SMALL  ONES 
1  =  1  +  1 
IF  I=J  THEN  10170 

IF  XCPd))  <=  X<P<M1))  THEN  10070 

REM  FIND  A  SMALL  ELEMENT  AMONG  THE  LARGE  ONES 

J=J-1 

IF  I=J  THEN  10170 

IF  X(P(J))  >=  X'PvMD)  THEN  10110 
REM  EXCHANGE  OUT  OF  PLACE  ELEMENTS 
H=P(I)  :  P(I)=P(J)  :  P(J)=H  :  GOTO  10070 
REM  ADJUST  COMPARE  ELEMENT  TO  NEW  MIDDLE 


10010 

10020 

1 0030 

1 0040 

SPECIAL 

1 0050 

10060 

1 0070 

10080 

1 0090 

10100 

10110 

10120 

10130 

10140 

10150 

10160 

ELEMENT 

10170 

1 0 1 30 

10190 

1 0200 

LARGE 

10210 

1 0220 

1 0230 

1 0240 

10250 

10260 

LARGE 

1 0270 

1 0280 

1 0290 

11000 

11010 

11020 

11030 

1 1 040 

1 1 050 

1 1 060 

11070 

11080 

11090 

12000 


IF  I  >  Ml  THEN  1=1-1 

IF  J=M1  THEN  10210 

H=P(I)  :  P(I)=P(M1)  :  P(M1)=H 

REM  SAVE  STARTING  POINT  FOR  ARRAY  SEGMENT  OF 
ELEMENTS 

L=LH-1  :  B<L)=I 

GOTO  1 0020 
IF  J-M  <  2  THEN  10270 

IF  X(P(M))  <=  X(P(M-+-1))  THEN  10260 
l-i=P(M)  -  p<M)=P(M+l)  :  P(M+1)=H 
REM  SET  BEGIN  AND  ENDPOINTS  FOR  ARRAY  SEGMENT 
ELEMENTS 
M=B(L)-t-l  :  L=L-1 
IF  L>0  THEN  GOTO  10020 
RETURN 

REM   *  *  *   READ  DATA  FILE  DIRECTORY   .*  t 
CLS:PRINT320, "E  XPONENTIAL     TE 
F*="EDITDATA/DIR: "+DN$ 
OPEN  " I " , 1 , F$ 
INPUT  #1,F0 
FOR  1=1  TO  FO 

INPUT  #1,F$(I) 
NEXT  I 
CLOSE  1 
RETURN 
REM  t    t     t       LIST  DATA  FILE  DIRECTORY  t     % 


OF 


S  T": PRINT 
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12010  CLS:PRINT320, "E    XPONENTIAL  TES    T":PRINT 

12020  FOR  1=1  TO  FO  STEP  4 

12030      FDR  J=l  TO  4 

12040  IF  I+J-1>F0  THEN  12090 

12050  IF  J=l  AND  LEN (F$ < I+J-1 ) X 12  AND  I>1  THEN 

PRINT  "" 

12060  PRINT  TAB(16*(J-1) ) ; " (";STR^(I+J-1) ; ") 

";F*(I+J-1) ; 

12070      NEXT  J 

12080  NEXT  I 

12090  RETURN 

19000  'ERROR  HANDLING 

19010  IF  ERR=106  AND  <SL=2005  OR  SL=2410  OR  SL=2710  OR  3110 

OR  SL=3510  OR  SL=3910)  THEN  INPUT  "NO  FILES  ON  THIS  DRIVE. 

PRESS  ENTER  TO  CONTINUE";A$  :  SL=0  :  RESUME  11080 

19100  PRINT  "ERR=";ERR; "ERL=";ERL  :  ON  ERROR  GOTO  O 


87 


f=-f;:ogf£<^m  *^  ,  c::h  i  sgjlj^o*^^ 


5  CLEAR  1000 
10  ON  ERROR  GOTO  19000 
20  DEFSTR  2:  DEFINT  P 

30  LB$=CHR$ ( 1 23 ) : RB*=CHR$  < 1 25 ) : MB*=CHR * (176): BL$=CHR$ ( 252 ) 
40  DIM  F$(95) 

50  LB$=CHR$ ( 123) : RB*=CHR$ ( 125) : MB*=CHR$ ( 176) : 3L*=CHR$ (252) 
60  REM  *  *  :*   CHI -SQUARE  TEST  *  * 

70  CLS: PRINT  320, "C  HI-SQUARE      TEST" 
30  PRINT:PRINT:PRINT  "Disk  Drive  Containing  Data  Files 
(0-3)?  "; 

90  DN$=INKEY*: IF  DN$=""  THEN  90  ELSE  PRINT  DN$:  I=VAL(DN$) 
100  IF  KO  OR  I>3  THEN  GOTO  SO 

110  SL=110  :  GQSUB  11000  :  IF  SL=0  THEN  70  ELSE  GOSUB  12000 
:  PRINT  "" 

120  PRINT: INPUT  "Enter  the  Data  File  #";F1 
130  IF  F1>F0  OR  FKl   THEN  120 
140  OPEN  "I", 1,F$(F1)+": "+DN$ 
150  INPUT  #1,RT: INPUT  #i,CT: INPUT  #1,BT 
160  DIM  X (RT,CT) ,R(RT) ,C(CT) 
170  FOR  Rl=l  TO  RT 
ISO  FOR  Cl=l  TO  CT 
190  FOR  Bl=l  TO  BT 
200  INPUT  #1,X(R1,C1) 
210  NEXT  Bl 
220  NEXT  CI 
230  NEXT  Rl 

240  N=0 : DF= ( RT- 1 ) * ( CT- 1 ) 
250  FOR  1=1  TO  RT 
260  R ( I ) =0 
270  FOR  J=l  TO  CT 
2S0  R ( I ) =R ( I ) +X ( I , J )     ^ 
290  NEXT  J  _ 

300  N=N-!-R  (  I )  ""o/ 

310  NEXT  I 
320  FOR  J=l  TO  CT 
330  C(J)=0 
340  FOR  1=1  TO  RT 
350  C(J)=C(J)+X (I, J) 
360  NEXT  I 
370  NEXT  J 
3S0  T=0:C2=0:CF=0 
390  FOR  1=1  TO  RT 
400  FOR  J=l  TO  CT 

4 1 0  E=R  ( I )  t  C  ( J  )  /N  -  ^'^ ^'"■^  °  ■'^^    ^^  -'"'■' "  ' '  •  -^  -^  CUJ22-. 
420  T=T+( (X (I, J)-E) C2) /E 
430  IF  E  <  5  THEN  CF=CF+1 
440  IF  E  <  1  THEN  CZ=CZ+1 
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450  NEXT  J 

460  NEXT  I 

465  IF  RT=2  AND  CT=2  THEN  T= ( ( ABS ( X (1 , 1 ) *  X (2, 2) -X ( 1 , 2) 

X (2, 1) )-0.5«N) C2) *  N/(R(1)*  R(2)*  C ( 1 ) «  C(2)) 

470    CLSrPRINT    320, "C    HI-SQUARE  TEST" 

430  PRINTi PRINT: PRINT  "Test  Statistic  =  " ; T 

490  PRINT  "Degrees  o-f  Freedom  =";DF 

500  PRINT  "Cells  with  expected  -Frequency  <  5  is  ";100 

CF/(RT*  CT);:  PRINT  "7." 

510    PRINT    "Number    c-f    cells    with    expected    -frequency    <     1     is 

7    -"- 

520  BOSUB  2000 

530  PRINT  "P-Value  =  ";:PRINT  USING  "#.####" ;QOX 

540  PRINT: INPUT"Enter  Desired  Alpha  Level";AL 

550  IF  AL<0  OR  AL>1  THEN  PRINT  "%%%    ERROR  %%%    ENTER  A  VLAUE 

BETWEEN  O  S<  1  "  :  GOTO  540 

560  PR I NT: PR I NT  "For  Alpha  =";AL;",  "; 

570  IF  AL<QOX  THEN  PRINT  "accept";  ELSE  PRINT  "reject"; 

5S0  PRINT  "  the  hypothesis  that  the  column  cells" 

590  PRINT  "have  equal  probabilities." 

600  PRINT  3960, "Press  SPACE  BAR  when  ready  to  return  to 

main  menu"; 

610  IF  INK"EY$=  ""  GOTO  610  ELSE  RUN  "MENU/BAS" 

2000  REM  %    %    t       CHI -SQUARE  DISTRIBUTION  *  * 

2020  CX=T 

204O  IF  DF>30  GOTO  2460 

2060  DP=1 

2080  FOR  I=INT(DF)  TO  2  STEP  -2 

2100  DP=DP 

2120  NEXT  I 

2140  NP=CXE  ( INT  (  (DF-i-1 )  /2)  )  *EXP  (-CX/2)  /DF^ 

2160  IF  INT(DF/2)--=DF/2  GOTO  2220^^ 

2130  L1=SQR  (2/CX/3.  141593)  '/^^      "^ 

2200  GOTO  2240  '^^     '^'=^,, 

2220  L 1  =  1  "*" 

2240  LC=1 

22±>0  MC=1 

22S0  D=DF 

2300  D=D+2 

2320  MC=MC»CX/D 

2340  IF  MC<.000001  GOTO  2400 

2360  LC=LC-<-MC 

23S0  GOTO  2300 

2400  F0X=L1*NP:*  LC 

2420  Q0X=1-F0X 

2440  RETURN 

2460  NX=( (CX/DF) E (1/3)  -  (l-(2/(9*  DF) ) ) ) /SQR (2/ (9*  DF) ) 

24S0  MU=0  :  VAR=1  :  12=1 

3000  REM    *:«*    NORMAL    DISTRIBUTION 

*  * 

3010  SD=3CR(VAR) 

3020  XN=(NX-MU) /SD: AX=AES ( XN) 
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3030  X2=. 3989423*  EXP ( (-AXtAX) /2) 

3040  NT=1/ (1+. 2316419*  AX) 

3050 

Q0X=X2* ( ( . 3193ai5*NT) - ( . 3565638*NTC2) + ( 1 . 7S147S*NTC3) - ( 1 . 821 

256*NTC4) + ( 1 . 330274*NTn5) ) 

3060  IF  XN<0  GOTO  3090 

3070  FQX=1-Q0X 

3080  RETURN 

3090  FOX=QOX 

3100  D0X=1-F0X 

3110  RETURN 

11000  REM   *  «  *   READ  DATA  FILE  DIRECTORY   *  * 

11010    CL3:PRINT    320, "C    HI-SQUARE  TES    T":PRINT 

11020  F*="CTABLE/DIR: "+DN$ 

11030  OPEN  "I", 1,F$ 

11040  INPUT  #1,F0 

11050  FOR  1=1  TO  FO 

11060       INPUT  #1,F$(I) 

11070  NEXT  I 

11080  CLOSE  1 

11090  RETURN 

12000  REM    *  *  *   PRINT  OUT  DIRECTORY  CF  DATA  FILES   *  * 

12010    CLS: PRINT    320, "C    HI-SQUARE  TES    T": PRINT 

12020  FOR  1=1  TO  FO  STEP  4 

12030      FOR  J=l  TO  4 

12040  IF  I+J-1>F0  THEN  12090 

12050  IF  J=l  AND  LEN ■ F$ ( I+J-l ) )< 12  AND  I>1  THEN 

PR I NT  "  " 

12060  PRINT  TAB(16*  (J-1)  )  :  "  C'jSTR^d-^J-i)  ;  ") 

".:F$<I+J-1)  ; 

12070      NEXT  J 

120S0  NEXT  I 

12090  RETURN 

19000  'ERROR  HANDLING 

19010  IF  ERR=106  AND  (SL=2005  OR  SL=2410  OR  SL=2710  OR  3110 

OR  SL=3510  OR  SL=3910)  THEN  INPUT  "NO  FILES  ON  THIS  DRIVE. 

PRESS  ENTER  TO  CONTINUE"; AS  :  SL=0  :  RESUME  11080 

19100  PRINT  "ERR=":ERR; "ERL=";ERL  :  ON  ERROR  GOTO  O 


Qn 


R-FRCDGi=?*^!^ 


1  O  , 


|NJF=-^-r*^"r/Brf=%S 


RB$=CHR*(12; 


MS$=CHR«(176) 


10  CLEAR  350 

20  DEFINT  I-K  :  DEFSTR  Z 

30  DEFDBL  F,L 

40  LB$=CHR$(123) 

BL$=CHR*(252) 

50  CLS  :  PRINT  TAB  (20)  ;  "  «**   MENU   :***" 

60  PRINT" (1)  Probability  Distributions 

70  PRINT" (2)  Inverse  Probability  Distributions 

SO  PRINT  "<3)  Return  to  Program  Selection  Menu" 

90  KP=704  :  KS=3   :  GOSUB  130 

100  DN  10  GGSUB  500,1000,120 

110  GOTO  50 

120  RUN  "MENU/BAS" 

130  REN 

140  PRINT  5)KP+5,  "OPTION  DESIRED  "  ;  LB$;  NB$;  R3$; 

150  ZI=INKEY$  :  IF  ZI=""  THEN  150EL3E  IQ=VAL(ZI) 

160  PRINT  ;5KP+21,  ZI;:  FOR  1  =  1  TO  100  :  NEXT 

170  IF  I0<1  OR  IO>KS  THEN  PRINT  S)KP,  BL$;  :  PRINT 


3KP 


1 000 


NEXT 


PRINT 


NOT  A  VALID  OPTION    :*:*:«";:  FOR  1  =  1  TO 

QKP,     BL$;  :  GOTO  140 

ISO  RETURN 

190  PRINT  a:'KP+5,L3$;  "  ";L*;"  or  ";R$;" 

"; LB$:MB$;RB$; 

200  ZQ=INKEY$  :  IF  ZQ=""  THEN  200ELSE  PRINT  :DKP+13,  Z3",  :  FDR 

1=1  TO  100  :  NEXT 


" ; RB$ ; " 


210  IF  ZQ<>L$  AND  ZQOCHR*  ( ASC  (L$ 

ZQ<>CHR$(A3C(R$)-32)  THEN  190 

220  RETURN 

500  REM 

510  CLS  :  PRINT 

520  PR I NT : PR I NT 

PR I NT "MENU: " 

PRINT"  (1)        Nor.Tial 

PRINT" (2) 


530 
540 
550 
560 
570 


TAB (13)" PROBAB I L I T Y 

T  7'=0 


Di  stri  buti  on 


'2)  AND  ZQOR*  AND 


DISTRIBUTIONS 


Wilcoxson  Signed— Ranks  Distribution" 
Mann— Whitney  Distribution" 
Kologmorov  Distribution" 
Sfnirnov  Distribution" 

Lilli-for's  Normality  Test  Distribution" 
Li  1 1  i  e-F or  "■■  5  Exponent  i  al  i  ty  Test 


PRINT  "(3 

PRINT  "  (4 
5S0  PRINT  " (5 
590  PRINT  " (6 
600  PRINT  " <7 
Di  stri but ion" 

610  PRINT  "(S)   Chi-Square  Distribution" 
620  KP=S96  :  KS=a  :  GOSUB  130 

630  ON  10  GOSUB  2000 , 3000 , 5000 , 6000 , 7000 , SOOO , 9000 , 1 0000 
640  IF  FLAG=1  THEN  FLAG=0: F0X=1-F0X : Q0X=1-QGX 

650  PRINT  :  PRINT  :  PRINT  "Prob.   (  X  <=";PX;")  =  ";  :  PRINT 
USING"#. ####"; FOX 
660  PRINT  "Prob.   (  X  >  ":PX;")  =  ";  :  PRINT 
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USING"#.  ####'* ;  QOX 

670    PRINT    3960, "       Run    Again    or    Quit?";:     KT-97S    :     L$="r"     : 

R$="q"     :     GOSUB    190 

6e0  CLS  :  I 2=0  :  IF  ZQ="r"  OR  ZQ="R"  GOTO  630ELSE  RETURN 

690  GOTO  50 

1000  CLS  :  PRINT  TAB < 15) " INVERSE  PROBABILITY  DISTRIBUTIONS 

1010  PRINT  5)192,  "MENU:  ":  PRINT"  (1)   Inverse  Normal 

1020  PRINT" (2)   Inverse  Wi 1  cox son-Si gned  Ranks" 

1030  PRINT  "(3)   Inverse  Mann-Whitney  Distribution" 

1040  PRINT  "(4)   Inverse  Kolmogorov  Distribution" 

1050  PRINT  "(5)   Inverse  Smirnov  Distribution" 

1060  PRINT  "(6)   Inverse  Lillie-for  Normality  Test 

Distribution" 

1070  PRINT  "(7)   Inverse  Lillie-for  Exponential  Test 

Distribution" 

lOSO  PRINT  "(8)   Inverse  Chi-Square  Distribution" 

1090  KP=a96:KS=S: GOSUB  130 

1 1 00  I 2=0 

1110  ON  10  GOSUB 

13000, 14000, 15000, 16000, 17000, ISOOO, 19000, 19220 

1120  PRINT  :  PRINT  "For   F(k)  =  Prob .   (  X  <=  a     )  =";PR 

1130  IF  XQ=9E36  THEN  PRINT  ">i    =   infinity"  ELSE  IF  XQ=9E-36 

PRINT  "m  =   -  infinity" 

1140  IF  XQ=9E36  OR  XQ=9E-36  THEN  GOTO  1160 

11 50  PRINT  ";:  =  ";  :  IF  ABS(XQ)>1000  THEN  PRINT  INT(XQ) 

ELSE  IF  ABS(XQ)>100  PRINT  USING"####. #" : XQ  ELSE  IF 

ABS(XQ)>10  PRINT  US I NG "###.##"; XQ  ELSE  PRINT 

USING"##. ####" : XQ 

1160  PRINT  3960,"   Run  Again  or  Cuit"::  L$="r"  :  R$="q"  : 

KP=976  :  GOSUB  190 

1170  IF  ZQ="R"  OR  Za="r"  GOTO  llOOELSE  RETURN 

2000  REM  t     t    t       NORMAL  DISTRIBUTION   *  * 

2010  IF  12=1  GOTO  2070 

2020  CLS: PRINT  320, "NORMAL  DISTRIBUTION" 

2030  PRINT:PRINT: INPUT"Enter  Mean " : MU 

2040  INPUT"Enter  Var i ance" : VAR 

2050  IF  VAR<=0  PRINT"*:*  ERROR  **    ";VAR"IS  NOT  A  VALID 

VARIANCE": GOTO  2040 

2060  INPUT"Enter  Normal  Statistic   x  " ; NX  :  IF  12=0  THEN 

PX=NX 

2070  SD=SQR(VAR) 

20S0  XN= (NX-MU) /SD: AX=ABS ( XN) 

2090  X  Z  = . 3989  423  X  P (  ( - A  X  *  A  X ) / 2 ) 

2100  NT=1/  (H-.2316419*AX) 

2110 

QOX=XZ« ( < . 3193815*NT) - < . 356563S*NTC2) + ( 1 . 78147S*NT:3) - ( 1 . S21 

256:*NTC4)-i-(l.  330274 «NTC5)  ) 

2120  IF  XN<0  GOTO  2150 

2130  F0X=1-Q0X 

2140  RETURN 

215C  FQX=QOX 

2160  QOX=l-FaX 
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2170  RETURN 

3000  REM  *  4:  *   WILCOXSON  SIGNED-RANKS  DISTRIBUTION   *  * 

3010  CLS:PRINT5>10, "WILCOXSON  SIGNED-RANKS  DISTRIBUTION" 

3020  PR I NT: PR I NT 

3030  INPUT  "ENTER  Wilcox sen  Stat i sti c " ; V: PX=V 

3040  IP  V  <  O  THEN  PRINT  "*  *  *  ERROR  *  *  *   STATISTIC  MUST 

BE  POSITIVE": GOTO  3030 

3050  INPUT  "ENTER  Sample  Size"; IN 

3060  IF  IN  <  1  THEN  PRINT  " *  t    %    ERROR  *  *  *   SAMPLE  SIZE 

MUST  BE  >  l":GaTO  3050 

3070  IF  V>IN:«  (IN+l)/2  THEN  PRINT  "  %    %    ERROR  %     t        " :  GOTO 

3030 

3080  IF  V>IN*  (IN+1) /4  THEN  V=IN:*  ( IN+l  ) /2-V-l  :  FLAG=1 

3090  IF  V>27  OR  (IN>9  AND  V>IN)  THEN  GCSUB  3620  ELSE  G0SU3 

3110 

3100  RETURN 

3110  IF  INT(y)=V  GGSUB  3150  ELSE  GOTO  3130 

3120  RETURN 

3130  V=INT(V) : GGSUB  3150: TFOX=FOX : V=V+1 : 60SUB 

3150: FOX= (FOX+TFOX ) /2: Q0X=1-F0X 

3140  RETURN 

3150  II (2)=0: II <3)=0: II (4)=0: II (5)=0: II (6)=0 

3160  C=l 

3 1 70  N=  I  NT  ( .  5*  ( - 1 +SQR  <  1  +Bt^J )  )  ) 

3100  FOR  1=1  TO  N 

3190  I1(1)=0 

3200  Il<l)=Il(l)-!-l 

3210  IF  I<2  GOTO  3360 

3220  1 1  <  2 ) = 1 1  ( 1 ) 

3230  Il(2)=Ii(2)+l 

3240  IF  I<3  GOTO  3360 

3250  II (3) =11 (2) 

3260  I  1  (3)  =11  ■:3>+l 

3270  IF  I<4  GOTO  3360 

3280  1 1 <  4 ) = 1 1  ( 3 ) 

3290  II  (4)  =11  (4) +1 

3300  IF  I<5  GOTO  3360 

3310  II (5) =11 (4) 

3320  II  (5)  =11  ':5)+l 

3330  IF  I<6  GOTO  3360 

3340  II  (6) =11  (5) 

3350  II  (6)  =11  (6)-^l 

3360  S=I1  (l)-t-Il  (2) +11  <3)+Il  (4)+I  1  •5)-!-Il  (6) 

3370  IF  S>V  THEN  GOTO  3390 

3380  C=C+1:0N  I  GOTO  3580,3570,3560,3540,3540,3530 

3390  IF  I>1  GOTO  3410 

3400  1 1  ( 1) = I N : GOTO  3580 

3410  FOR  IL=1  TO  I-l 

3420  S=0:J=I-IL 

3430  IF  J-1<=0  THEN  GOTO  3470 

3440  FDR  K=l  TO  J-1 

3450  8=5+ 1 1 ( K ) 


3460  NEXT  K 

3470  S=S+(IL+1)  1  ':J)  +  (IL+1)  *  (IL-J-2)  /2 

3480  IF  S>V  THEN  GOTO  3510 

3490  I1(I+1-IL)=IN 

3500  ON  J    GOTO  3570,3560,3550,3540,3530 

3510  NEXT  IL 

3520  I1(1)=IN:G0T0  3580 

3530  IF  I1(6)<IN  GOTO  3350 

3540  IF  I1(5)<IN  GOTO  3320 

3550  IF  I1(4)<IN  GOTO  3290 

3560  IF  I1(3)<IN  GOTO  3260 

3570  IF  I1(2)<IN  GOTO  3230 

3580    IF    IldXIN    GOTO    3200 

3590  NEXT  I 

3600  F0X=C/2!:  IN:  QOX=1-FOX 

3610  RETURN 

3620  REM  t    t    t    WILCOX SON  SIGNED  RANK  APPRO X  DISTIBUTICN  * 

3630  MU=IN« (IN+1) /4 

3640  V AR= ( 2N+ 1 ) * ( I N+ 1 ) N/ 24 

3650  NX=V+0.5 

3660  GOSUB  2070 

3670  L4=  ( 3N C  2+3N- 1  >  /  (  1  ON*  ( I N+ 1 )  :*  ( 2N-i- 1  )  ) 

3680  W=<NX-MU)  /SQR(VAR)  :  F=  <EXP  (-0.  5*UJ[:2)  )  /SQR  <6.  28313) 

3690  F3=F* (WE3-3*W) 

3700  F0X=F0X+L4*  F3 

37 1 0  F5=-F  *  ( W  E  5- 1 0  * W  C  3+ 1 5  :* W ) 

3720  F7=-F:*  ( WC7-21  ;*WC5-i-105:*WC3-105*W) 

3730 

L6=4* (3NC4+6NC3-3N+1 ) / (35* (IN* (IN+1 ) * (2N+i) ) C2) 

3740  F0X=F0X+L6*  F5+.5*F7:*  L4l2 

3750  aGX=l-FOX 

3760  RETURN 

5000  REM   *  *  *   MANN-WHITNEY  DISTRIBUTION   ^c  * 

50 1 0  CLS :  PR  I  NT:D  1 0 ,  "  MANN-WH I TNEY  D I STR I  BUT  I GN  " 

5020  PR I NT: PR I NT 

5030  INPUT  "ENTER  Mann-Whitney  Stat i sti c" ; U: PX=U 

5040  IF  U  <  O  THEN  PRINT  "*  *  *   ERROR   *  *  *   STATISTIC 

MUST  BE  POSITIVE  " : GOTO  5030 

5050  INPUT  "ENTER  Smaller  Sample  Size";N 

5060  IF  N  <  1  THEN  PRINT  "«  *  *  ERROR  *  *  *  SAMPLE  SIZE 

MUST  BE  >  1":G0T0  5050 

5070  INPUT  "ENTER  Larger  Sample  Size":M 

5080  IF  M  <  N  THEN  PRINT  "  *  *  *   ERROR  *  *  *  YOU  ENTERED 

THE  SMALLER  SAMPLE  SIZE  LAST  " : GOTO  5050 

5090  FLAG=0 

5100  IF  U>N*M/2-.5  THEN  U=N*M-U-1 : FLAG=1 

51 10  IN=N+M:  V=U-(-N*  (N+1 )  /2:  I=N 

5120  IF  U>N*M   THEN  PRINT  "  *  *  ERROR  *  *   THE  STATISTIC 

MUST  BE  >  THE  PRODUCT  OF  THE  SAMPLE  SIZES": GOTO  5030 

5130  IF  V>27  OR  M>9  THEN  GOSUB  5650  ELSE  GOSUB  5150 

5140  RETURN 

5150  IF  INT(V)=V  THEN  GOSUB  5190  ELSE  GOTO  5170 


94 


5160  RETURN 

5170  V=INT(V) :GOSUB  5190: TFOX=FOX : V=V+1 : GOSUB 

5190:FGX=(F0X+TF0X)/2:Q0X=1-FQX 

5180  RETURN 

5190  II (2)=0: II (3)=0: II (4)=0: II (5)=0: II (6)=0 

5200  C=0 

5210  I1(1)=0 

5220  II  (1)=I1  <l)-«-l 

5230  IF  I<2  GOTO  5370 

5240  II (2)=I1 (1) 

5250  II (2)=I1 (2)+l 

5260  IF  I<3  GOTO  53SO 

5270  II (3) =11 (2) 

5280  II (3)=I1 (3)+l 

5290  IF  I<4  GOTO  5380 

5300  11 (4) =11 (3) 

5310  II (4)=I1 (4)+l 

5320  IF  I<5  GOTO  5380 

5330  I1':5)=I1<4) 

5340  II (5) =11 (5)+l 

5350  IF  I<6  GOTO  5380 

5360  II (6) =11 (5) 

5370  11 (6) =11 (6) +1 

5380  8=11  (D+Il  (2) +11  (3) +11  (4) +11  (5) +11  (6) 

5390  IF  S>V  THEN  GOTO  5410 

5400  C=C+1:0N  I  GOTO  5600,5590,5530,5570,5550,5550 

5410  IF  I>1  GOTO  5430 

5420  1 1  ( 1 ) = I N : GOTO  5600 

5430  FOR  IL=1  TO  I-l 

5440  S=0:J=I-IL 

5450  IF  J-1<=0  THEN  GOTO  5490 

5460  FOR  K=l  TO  J-1 

5470  8=8+11 <K) 

5480  NEXT  K 

5490  S=S+(IL+1) 1 (J) +( IL+1 ) * (IL+2) /2 

5500  IF  S>V  THEN  GOTO  5530 

5510  II (I+1-IL)=IN 

5520  ON  J  GOTO  5590,5580,5570,5560,5550 

5530  NEXT  IL 

5540  II (1)=IN:G0T0  5600 

5550  IF  I1-:6)<IN  GOTO  5370 

5560  IF  I1(5)<IN  GOTO  5340 

5570  IF  I1<4)<IN  GOTO  5310 

5580  IF  I1':3)<IN  GOTO  5280 

5590  IF  I1':2)<IN  GOTO  5250 

5600  IF  lidXIN  GOTO  5220 

5610  8=1 

5620  FOR  K=0  TO  N-1 : 8=8* (IN-K) / CN-K) : NEXT  K 

5630  F0X=C/8:Q0X=1-FQX 

5640  RETURN 

5650  REM  *  :*  *  MANN-WHITNEY  APPROXIMATE  DISTRIBUTION  *  * 

5660  MU=N/2: VAR=N« (N+M+1 > / 12: NX=U+. 5 


5670  GGSU3  2070 

5680  L4=(MC2+NE2+IN-!-M*N)  /  (20:*N:*<INh-1)  ) 

5690 

L6=(2*  (MC4-i-NC4)-i-4*N:*  (ME2-!-Ni:2) +6*  •M!:2)  :*  <Nl2)+4*  (MCS+NCS) 

*N*  (M-!-N)  +MC2+NC2-H2:*N-IN)  /  (210*  (MC2)  *  ':Nr2)  *  <  { IN+1  )  112)  ) 

5700  W= ( N X -MU ) / SDR ( VAR ) : F= ( E  X  P ( -0 . 5  E  2 ) ) / SQR ( 6 . 283 1 S ) 

57 1 0  F3=F  * ( W  C  3-3 ) 

5720  F5=-F*  (WC5-10!:3+15) 

5730  F7=-F  « ( W  C  7-2 1 C  5+ 1 05  C  3- 1 05 ) 

5740  F0X=F0X+L4«  F3 

5750  F0X=FDX+L6;*  F5 

5760  F0X=F0X+.5*  F7*  L4C2 

5770  RETURN 

6000  REIi   *  *  *   KOLMOGROV  DISTRIBUTION   :*  :* 

6010  CLS:PRINT310,  "KGLr-tOGOROV  DISTRIBUTION" 

6020  PR I NT: PR I NT 

6030  INPUT  "ENTER  SKclmcgcrcv  Statistic" ;D 

6040  IF  D<0  OR  D>1  THEN  PRINT  "*  *  ERROR  *  * " : GOTO  6030 

6050  INPUT  "ENTER  Sample  Sizs";N 

6060  IF  N<1  THEN  PRINT  "  :*  *  ERROR  *  * " :  GOTO  6050 

6070  S=D  :*  ( 0 .  1 2+SQR  ( N )  +0 .  1  1  /  SQR  ( N )  ) 

6080  GOSUB  6110 

6090  PX=D 

6100  RETURN 

6110  aCX=(.  135)  L  (Si:2) 

6120  IF  QOX<0  THEN  QOX=0 

6130  IF  Q0X>1  THEN  Q0X=1 

6140  FGX=1-Q0X 

6150  RETURN 

7000  REM  t    t     t       SMIRNOV  DISTRIBUTION   *  * 

7010  CLS:PRINT310,  "SMIRNOV  DISTRIBUTION" 

7020  PR I NT : PR I NT 

7030  INPUT  "ENTER  Smirnov  Statistic";D 

7040  IF  D<0  OR  DM  THEN  PRINT  "*  *  ERROR  *  * "  :  GOTO  7030 

7050  INPUT  "ENTER  Smaller  Sample  Size";N 

7060  IF  N<1  THEN  PRINT  "*  t    ERROR  t     t       SAMPLE  SIZE  MUST 

BE  >  1":G0T0  7050 

7070  INPUT  "ENTER  Larger  Sample  Size":M 

7080  IF  M<N  THEN  PRINT  "t     *  ERROR  t     *   YOU  ENTERED  THE 

SMALLER  SAMPLE  SIZE  LAST": GOTO  7050 

7090  IF  N<  =  .  1  THEN  S=  (D-1  /  (2)  )  *  ( .  12+3QR  (N)  ^-.  1 1 /SQR  <M)  ) 

ELSE  IF  (M/N)=INT(M/N)  THEN 

3=D*  SQR'N/ (N-!-M)  )+2/ <33c  SQR(M)) 

ELSE  S=D*  SQR  (N/  (N+M)  )  -t-2/  <5*  SQR  (M)  ) 
7100  GOSUB  6110 
7110  PX=D 
7120  RETURN 

8000  REM   *  *  -*   L  ILL  IFOR'S  TEST  FOR  NORMALITY   *  i 
8010  J=IO 

8020  CLS: PRINT  "SELECT  OPTION" : PRINT: PRINT 
8030  PRINT  "<:i)    MEAN  and  VARIANCE  Br&    unknown ":  PRINT 
8040  PRINT  "(2)    MEAN  is  known.  VARIANCE  is  unknown" 
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e050  KF-S96:KS=2:GDSU3  140 

3060  CLS: PRINTS 10,  "LILLIFOR'S  TEST  FOR  NORMALITY" 

8070  PR I NT: PR I NT 

8080  INPUT  "ENTER  SKcIjt.dq or ov  Statist!  c";D 

3090  IF  D<0  OR  D>1  THEN  PRINT  "t    t    ERROR  *  * " : GOTO  SOSO 

8100  INPUT  "ENTER  Sainple  Si2&";N 

8110  IF  N<1  THEN  PRINT  " «  *  ERROR  *  »   SAMPLE  SIZE  MUST  BE 

POSITIVE": GOTO  8100 

8120  DEN0M=1.0:PX=D 

8130  D=PX/DENOM 

3140  IF  10=1  THEN  8=0* (-. Ol+SQR (N) +0. S5/SQR (N) )  ELSE 

S=D:*  ( O.  055+SQR  <N)  +.  2S3/SQR  (N)  ) 

8150  TDEN0M=DEN0M 

SI 60  GOSUB  6110 

8170  IF  I0>1  THEN  GOTO  8260 

8180  IF  FOX  <=  .85  THEN  DEN0M=i--24«  FOX: GOTO  3230 

8190  IF  FOX  <=  .90  THEN  DENOM= 1.357-  .66*  FOX: GOTO  3230 

8200  IF  FOX  <=  .95  THEN  DENOM  =  1.339  -  .64:*  FOX  :  GOTO 

8230 

8210  IF  FOX  <=  .975  THEN  DENOM  =  1.795  -  1.12*  FOX  :  GOTO 

8230 

8220  DENOM  =  2.133-1.467*  FOX 

8230  IF  ABS(TDENOM-DENOM)  >  .001  THEN  GOTO  3130 

3240  10= J 

3250  RETURN 

8260  IF  FOX  <=  -9  THEN  DEN0M=1 . 83-. 3*  FOX  : GOTO  3230 

8270  IF  FOX  <=  .95  THEN  DENOM  =1.47  -  -4*  FOX  :  GOTO  3230 

8230  IF  FOX  <=  -975  THEN  DENOM  =  1.35  -  .3*  FOX  :  GOTO  3230 

3290  DENOM  =  1 . 07  :  GOTO  8230 

9000  REM   *  *  *   L ILL IFOR'S  TEST  FOR  EXPONENT I AL I TY   *  * 

9010  CLS: PRINTS 10,  "L ILL  IFOR'S  TEST  FOR  EXPONENT  I AL I TY" 

9020  PR I NT : PR I NT 

9030  INPUT  "ENTER  SKclmcgcrov  Stati5tic";D 

9040  IF  D<0  OR  D>1  THEN  PRINT  "*  *  ERROR  *  * " : GOTO  SOSO 

9050  INPUT  "ENTER  Sample  Sizs":N 

9060  IF  N<1  THEN  PRINT  "*  *  ERROR  *  *   SAMPLE  SIZE  MUST  BE 

POS I T I VE " : GOT  0  31 00 

9070  IF  D  >  1  .3442/ <SQR  (N) +.26+.5/SQR':N)  ) +.2/N  THEN  FR=1  : 

GOTO  9160 

9030  PR=.5:G0SUB  19090 

9090  IF  D>XQ  THEN  P1=PR: P2=l : GOTO  9110 

9100  P1=0:P2=PR 

9110  PR= ( P 1 +P2 ) / 2 

9120  GOSUB  19090 

9130  IF  ABS(D-XQ)<.000i  THEN  GOTO  9160 

9140  IF  D::XQ  THEN  P1=PR  ELSE  P2=PR 

9150  GOTO  9110 

9160  FOX=PR 

9170  Q0X=1-F0X 

9130  PX  =  D 

9190  RETURN 
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10010  IF  12=1  GOTO  lOOSO 

10020    CLS:     PRIlviT    5)20,  "CHI -SQUARE    DISTRIBUTION" 

10030    PRINT    3192,  "";  :  INPUT"Enter    Chi-Sqciare    Statistic       ;■; 

" :  CX 

10040  IF  CX<=0  PRINT"**  ERROR  **    CHI -SQUARE  STATISTIC 

MUST  BE  >  0":SOTO  10030 

10050  IF  12=0  THEN  PX=CX 

10060  INPUT"Entsr  Degrees  o-f  Freedom"  ;DF 

10070  IF  DF<=0  OR  DF<>INT(DF)  THEN  PRINT"  *t  ERROR  *:* 

ENTER  AN  INTEGER  >  0":GOTO  10060 

10080  IF  DF>30  GOTO  10290 

10090  DP=1 

10100  FOR  I=INT(DF)  TO  2  STEP  —2 

101 10  DP=DP 

10120  NEXT  I 

10130  NP=CXE (INT( (DF+1) /2) ) *EXP<-CX/2) /DP 

10140  IF  INT(DF/2)=DF/2  GOTO  10170 

10150  PC=SQR(2/CX/3. 141593) 

10160  GOTO  lOlBO 

10170  PC=1 

10180  LC=1 

10190  MC=1 

10200  D=DF 

10210  D=D+2 

10220  MC=MC«  CX/D 

10230  IF  MC<.0O0OOl  GOTO  10260 

10240  LC=LC+MC 

10250  GOTO  10210 

10260  FOX=PC*NP*  LC 

10270  Q0X=1-F0X 

10280  RETURN 

10290  NX=(  (CX/DF)  C  (1/3)  -  <  1- (2/ (9:SDF)  )  )  ) /SQR  <2/ (9=cDF)  ) 

10300  MU=0  :  VAR=1  :  12=1 

10310  GDSUB  2070 

10320  RETURN 

13000  REM  t    M    t  INVERSE  NORMAL  DISTRIBUTION   4c  * 

13010  IF  12=1  GOTO  13100 

13020  CLS  :  PRINT  320,  "INVERSE  NORMAL" 

13030  PRINT  :  PRINT  :  INPUT "Enter  Mean "; MO 

13040  INPUT  "Enter  Variance" ;VR 

13050  IF  VR<=0  PRINT" *t  ERROR  tt  ":VR:"IS  NOT  A  VALID 

VAR I ANCE " : GOTO  1 3040 

13060  INPUT"Enter  Probability  (i.e.   F(>;)  =   Prob.   <  X  <=  ,< 

)  =  ";PR 

13070  IF  PR=0  THEN  XQ=9E-36  ELSE  IF  PR=i  XQ=9E+36 

13080  IF  PR=0  OR  PR=1  RETURN 

13090  IF  PR<0  CR  PR>1  THEN  PRINT"**  ERROR  **    ";PR;"IS  NOT 

A  VALID  PROBABILITY": GOTO  13060 

131 00  I F  PR  > . 5  THEN  RL=SQR ( LOG ( 1 / (  ( 1 -PR )  C  2 )  )  )   ELSE 

RL=SQR ( LOG (  1 / ( PR  E  2 )  )  ) 

13110 

XQ=RL-( (2. 515517+ (.302S53*RL) + ( . 01032S*RLC2) ) / (1*(1. 4327SS*R 
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L)  +  (.  189269:*RLC2)  +  (.00130S*RLL3)  )  ) 

13120  IF  PR>=.5  THEN  XQ=MO+ (SQR < VR) *XQ)  ELSE 

XQ=MD-<SQR(VR)  :*XQ) 

13130  RETURN 

14000  REM  t    t    t        INVERSE  WILCOX  SON  SIGNED-RANKS  :*  * 

14010  CLS: PRINT  "INVERSE  WILCOXSON  SIGNED-RANKS 

DISTRIBUTION" 

14020  PRINT: PRINT 

14030  INPUT  "ENTER  Sample  Size"; IN 

14040  IF  IN  <  1  THEN  PRINT  "«  *  «   ERROR  *  *  *  SAMPLE  SIZE 

MUST  BE  POSITIVE": GOTO  14030 

14050  M0=IN:*  (  IN+1) /4-.5 

1 4060  VR= ( 2N+ 1 ) t  < I N+ 1) N/24 

14070  GOSUB  13060 

14071  V=INT(XQ> 

14072  GOSUB  3090 

14073  TFOX=FOX 

14074  V=V+1 

14075  GOSUB  3090 

14076  XQ=INT (XQ) + (CDBL (PR) -TFOX) / (FOX-CDBL (TFOX ) ) 
14080  MO=XQ-INT(XQ) 

14090  IF  M0<.26  THEN  XD=INT ( XD) : RETURN 

14100  IF  M0<.75  THEN  XQ=INT ( XD) +. 5  ELSE  Xa=INT(XQ)+l 

14110  RETURN 

15000  REM  :«  *  t        INVERSE  MANN-WHITNEY   *  * 

15010  CLS: PRINT  "INVERSE  MANN- WHITNEY  DISTRIBUTION" 

15020  PR I NT : P R I NT 

15030  INPUT  "ENTER  Smaller  Sample  Size":N 

15040  IF  N<1  THEN  PRINT  "*  *  t    ERROR  *  *  *  SAMPLE  SIZE  MUST 

BE  PCS I T I VE " : GOTO  1 5030 

15050  INPUT  "Enter  Larger  Sample  Size";M 

15060  IF  M<1  THEN  PRINT  "*  *  t    ERROR  *  *  *  SAMPLE  SIZE  MUST 

BE  POSITIVE": GOTO  15050 

15070  M0=N*M/2 

1 5080  VR=N:* M *  ( N+M+ 1  )  /  1 2 

15090  GOSUB  13060 

15091  U=INT':XQ) 

15092  GOSUB  5130 

15093  TFOX=FaX 

15094  U=U+1 

15096  GOSUB  5130 

1 5097  XQ=  I  NT  (  XQ )  -»-  ( CDBL  ( PR )  -CDBL  ( TFOX  )  )  /  ( FOX-CDBL  ( TFOX  >  ) 
151 00  MC=  X Q- 1 NT ( X  3 ) 

15110  IF  MG< -  26  THEN  XQ= I NT ( XQ )  : RETURN 

15120  IF  M0<.75  THEN  XQ=INT  (  XQ) +.  5  ELSE  XQ=INT  ( XQ) -i-1 

15130  RETURN 

16000  REM  t    t    t        INVERSE  KOLMOGDROV  *  :* 

16010  CLS  :  PRINT  5)20,  "INVERSE  KQLMGGQRQV" 

16020  PR I NT: PR I NT 

16030  INPUT "ENTER  Sample  Size";N 

16040  IF  N<1  THEN  PRINT  "*  *  ERROR  :lc  t " :  GOTO  16030 

16050  INPUT"'Enter  Prcbabilitv  (i.e.   F(;;)  =   Prob.   (  X  <=  x 
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=  ";PR 

6060  IF  PR=1  THEN  XQ=  0.0  ELSE  IF  PR=0  XQ=9E+36 

6070  IF  PR=0  OR  PR=1  RETURN 

60e0  IF  PR<0  GR  PR>1  THEN  PRINT"**  ERROR  Mt  ";PR;"IS  NOT 

A  VALID  PROBABILITY": GOTO  16050 

6090  S=SQR (LOG  < 1-PR) /LOG ( .  135)  ) 

6100  D=S/  (.  12+SQR(N)-f-.  11/SQR(N)  ) 

6110  XD=D 

6120  RETURN 

7000  REM  *  *  *   INVERSE  SMIRNOV  *  * 

7010  CLS  :  PRINT  5)20,  "INVERSE  SMIRNOV" 

7020  PR I NT: PR I NT 

7030  INPUT "ENTER  Smaller  Sample  Size";N 

7040  IF  N<1  THEN  PRINT  "*  «  ERROR  *  « " : GOTO  17030 

7050  INPUT "ENTER  Larger  Sample  Size";M 

7060  IF  M<N  THEN  PRINT  "*  *  ERROR  *  « " : GOTO  17030 

7070  INPUT "Enter  Probability  (i.e.   F(k)  =   Prob-   (  X  <=  x 
=    ";PR 

7080  IF  PR=1  THEN  XQ=  0.0   ELSE  IF  PR=0  THEN  XQ=9E+36 

7090  IF  PR=0  OR  PR=1  RETURN 

7100  IF  PR<0  OR  PR>1  THEN  PRINT"  ic*  ERROR  tt  ":PR;"I3  NOT 

A  VALID  PROBABILITY": GOTO  17060 

7110  S=SQR ( LOG ( 1 -PR ) / LOG ( . 1 35 ) ) 
17120  IF  N<=0.1*M  THEN  D=S/ ( .  12+SQR  (N) +.  1  1 /SDR  (N)  ) -:-l  /  (2*rl) 
ELSE  IF  •M/N)=INT(M/N)  THEN 

D= (S-2/ (3*  SCR (M) ) ) /SQR (N*M/ ( N+M) ) 
ELSE  D= (S-2/ (5*  SQR (M) ) ) /SQR (N*M/ (N+M) ) 
17130  XQ=D 
17140  RETURN 

ISOOO  REM  *  *  *   INVERSE  LILLIEFOR"*  S  NORMALALITY   *  * 
13010  J=ia: CLS: PRINT  "SELECT  OPTION" : PRINT: PRINT 
1S020  PRINT  "(1)    MEAN  and  VARIANCE  Ars    unknown  ': PRINT 
1S030  PRINT  "(2)    MEAN  is  kncv^n,  VARIANCE  is  unknown" 
1 8040  KP=896 : KS=2 : GOSUB  1 40 

13050  CLS:  PRINT  5)13,  "INVERSE  LILLIEFOR'S  NORMALITY  TEST 
DISTRIBUTION" 
1 3060  PR I NT : PR I NT 

18070  INPUT "ENTER  Sample  Size":N 

180S0  IF  N<1  THEN  PRINT  "*  *  ERROR  *  * " : GOTO  16030 
1S090  INPUT  "Enter  Probability  (i.e.   F(>;)  =   Prob.   (  X  <>  >i 
)  =  ";PR 

iSlOO  IF  PR=0  THEN  Xa=9E-36  ELSE  IF  PR= 1  XS=9E+36 
131 10  IF  PR=0  OR  PR=1  THEN  I a=-J:  RETURN 

1S120  IF  PR<0  OR  PR::1  THEN  PRINT"**  ERROR  **    ";PR:"IS  NOT 
A  VALID  PROBABILITY": GOTO  13060 
1 3 I 30  S=SQR ( LOG ( ( 1 -PR ) ) /LOG ( . 1 35 ) ) 

18140  IF  10=1  THEN  D=S/ (-.  01+SQR  <N) -h.  S5/SQR  (N)  )  ELSE 
D=S/  ( .  055+ .  283 /SQR  (N)  -i-SQR  (N)  ) 
18150  IF  I0>1  THEN  GOTO  18240 

18160  IF  PR  <=  .85  THEN  DEN0M=i-.24*  FOX : GOTO  1S210 
18170  IF  PR  <=  .90  THEN  DENOM= 1.357-  .66*PR:GQT0  18210 
18130  IF  PR  <=  .95  THEN  DENOM  =  1.339  -  0.64* PR  :  GOTO 
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1S210 
1  3190 
18210 
13200 
18210 
18220 
1 8230 
18240 
18250 
18260 
18270 
19000 
1 90 1 0 
DI3TR 
1  9020 
1 9030 
19040 
1 9050 
)  =  " 
19060 
1 9070 
1 9080 

A  VAL 
1  9090 

19100 
19110 
19120 
19130 
19140 
19150 
19160 
19170 
19180 
191 90 
1 9200 
19210 
i  9220 
19230 
20000 
200 1 0 
20020 
ENTER 
20030 
)  =  " 
20O40 
20050 
20060 
A  VAL 
20070 
20080 
20090 
20100 


IF  PR  <=  .975  THEN  DENOM  =  1.795 


DENOM  =  2. 133-1.467*PR 
XQ=D*DENQM 
10  =  J 
RETURN 
IF  PR  <=  , 
IF  PR  <=  , 
IF  PR  <=  , 
DENOM  =  1, 
REM  *  t    t 
CLS : PR I NT 
I BUT ION " 
PR I NT: PR I NT 

INPUT "ENTER  Sample  Size":N 
IF  N<i  THEN  PRINT  "*  t    ERROR 


1.12*PR  :  GOTO 


9  THEN  DENDM  =  1.83  -  .8*PR  : 
95  THEN  DENOM  =  1.47  -  . 4*PR 
975  THEN  DENOM  =  1.85  -  .  8:*PR 
07: GOTO  1S210 

INVERSE  LILLIEFOR  EXPONENTIAL  TEST  * 
:D13,  "INVERSE  LILLIEFOR' 3  EXPONENTIAL 


GOTO  18210 

GOTO  13210 
:  GOTO  18210 


* 

TEST 


INPUT"Enter  Probability  (i.e. 

;PR 
IF  PR=0  THEN  XQ=9E-36  ELSE 
IF  PR=0  OR  PR=1  RETURN 
IF  PR<0  OR  PR>1  THEN  PRINT 

ID  PROBABILITY": GOTO  19050 
DENOM=SQR (N) +. 26+. 5/SQR (N) 


*":GOTO  16030 
F(x)  =   Pr ob. 


A   -^  = 


IF  PR=1  XQ=9E+36 


:**  ERROR  ** 


";PR: "IS  NOT 


IF 
IF 
IF 
IF 


PR 
PR 
PR 
PR 


95  THEN  XQ  =  5 


90  THEN 
SO  THEN 


XQ  = 
XQ  = 


IF  PR 

IF  PR 


IF 
IF 
XQ 


PR 
PR 


7 

5 


1 


THEN 
THEN 
THEN 
THEN 
THEN 


XQ  = 
XQ  = 
XQ  = 
XQ  = 
XQ  = 
4062 


. i25*PR  -  3.7S0S 
2.026*PR  -  .8367 
1. 124* PR  -  .0249 
743*PR  +  .2799  : 
*-  .4137  ; 
. 45SS  : 


.546.*PR 
.464:icpR 
.485* PR 
.587:*PR 


.4525 


:  GOT 

i-> 

19190 

:  GOT 

0 

191 90 

:  GOT 

0 

191 90 

GOTO 

1 

9 1 90 

GOTO 

1 

19190 

GOTO 

± 

9 1 90 

GOTO 

19190 

GOTO 

1 

9190 

.846;*PR  -* 
XQ=XQ/ DENOM 
XQ  =  XQ  +  .  2/N 
RETURN 
REM 

IF  12=1  GOTO  20070 

CLS  :  PRINT  5)20.  "  INVERSE  CHI-SQUARE" 
PRINT  :  INPUT"Enter  Decrees  or  Freedom' 


r.jr 


IF  DF<=0  OR  DF<>INT<DF>  THEN 
AN  INTEGER  >  0":GOTO  20010 

'Enter  Probability    i, 


PRINT"*;*  ERROR  ** 


INPUT 

PR 

IF  PR=0 

IF  PR=1 


=  Prcb. 


Y  ;■  = 


THEN  XQ=0  :  RETURN 
THEN  XQ=9E36  :  RETURN 

IF  PR<0  OR  PR>1  THEN  PRINT"** 
ID  PROBABILITY": GOTO  20030 

M0=0  :  VR=1  :  12=1  :  P2=PR 

IF  DF>1  GOTO  20130 

PR=. 5* ( i-PZ) 

GOSUB  13100 


ERROR  ** 


" : PR ; " I 3 


NO" 
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20110  XQ=XQC2  :  PR=PZ 

20120  RETURN 

20130  IF  DF>2  GOTO  20160 

20140  XQ=-2*  LOG(l-PZ) 

20150  RETURN 

20160  GOSUB  13100 

20170  DQ=XQ 

20180  IF  DF>(2+INT(4.1cABS(DQ)  )  )  GOTO  202S0 

20190  X 9= 1.0000886  -  -  2237368 /DF  -  . 01513904/DFC2 

20200  X9=X9+(DFC-.5  %    DQ  *  (.4713941  +  . 026070S3/DF  - 

.00S986007/DFC2) ) 

20210  X9=X9  +  (1/DF  *  DQE2  *  (.000134S02S  +  .0112S1S6/DF  + 

. 02277679 /DFC2> ) 

20220  X9=X9  +  (DFC-1.5  *  DQC3  *  (-.008553069  -  .01 153761 /DF 

-  .01323293/DFi:2)  ) 

20230  X9=X9  +  (DFC-2  *  DQC4  *  (.00312558  +  . 005169654/DF  - 

.  006950356/ DF  1:2)  ) 

20240  X9=X9  +  (DFC-2.5  *  DQE5  :*  (-.0008426812  + 

. 0025300 1 / DF  +  . 00 1 060438 / DF  C  2 )  ) 

20250  X9=X9  +  (DFC-3  *  D0r6  t     (.00009780499  -  . 0014501 17/DF 

+  .001565326/DFC2) ) 

20260  XQ=DF*X9C3 

20270  RETURN 

20280  X9=l  -  2/(9*DF)  +  (4*DQC4  +  16*DQC2  -28) / ( 1215*DFE2) 

20290  X9=X9  +  (8*DQC6  +  720*DQC4  +  3126*DQC2  + 

2904)  /  ( 229635  *DF  1:3) 

20300  XS=DQ/3  +  (-DQC3  +  3*DQ) / ( 1624DF) 

20310  XS=X8  -  (3*DQC5  +    40*DQC3  +  45*DQ)  /  <5S32*DFi:2) 

20320  XS=X3  +  (301*DQC7  -  1517:*DQC5  -32769*DQC3  - 

79349:*DQ)  /  (7873200*DFC3) 

20330  Xa=  SQR(2/DF)  *  X8 

20340  X  Q= DF  *  ( ( X  8+  X  9 ) C  3 ) 

20350  RETURN 


lOJ. 
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